source: git/Singular/links/ssiLink.cc @ cfe30a

fieker-DuValspielwiese
Last change on this file since cfe30a was cfe30a, checked in by Andreas Steenpass <steenpass@…>, 11 years ago
chg: do another waitpid() after nanosleep() succeeded (cherry picked from commit 9d9a7faa6bb7522c1db8caabc4c47c0aa0840caf) Signed-off-by: Andreas Steenpass <steenpass@mathematik.uni-kl.de>
  • Property mode set to 100644
File size: 50.8 KB
Line 
1/****************************************
2 * Computer Algebra System SINGULAR     *
3 ****************************************/
4/***************************************************************
5 * File:    ssiLink.h
6 *  Purpose: declaration of sl_link routines for ssi
7 ***************************************************************/
8#include <stdio.h>
9#include <fcntl.h>
10#include <errno.h>
11#include <unistd.h>
12#include <sys/types.h>
13#include <signal.h>
14#include <sys/types.h>          /* for portability */
15#include <sys/select.h>
16#include <sys/socket.h>
17#include <netinet/in.h>
18#include <ctype.h>   /*for isdigit*/
19#include <netdb.h>
20#include <sys/wait.h>
21#include <time.h>
22
23#include <kernel/mod2.h>
24#include <Singular/si_signals.h>
25// #include "mod2.h"
26
27#include <Singular/tok.h>
28#include <Singular/ipid.h>
29#include <Singular/ipshell.h>
30#include <omalloc/omalloc.h>
31#include <libpolys/polys/monomials/ring.h>
32#include <libpolys/polys/matpol.h>
33#include <libpolys/polys/simpleideals.h>
34#include <libpolys/polys/monomials/p_polys.h>
35#include <libpolys/coeffs/longrat.h>
36#include <libpolys/misc/intvec.h>
37#include <libpolys/coeffs/bigintmat.h>
38#include <libpolys/misc/options.h>
39#include <kernel/timer.h>
40#include <Singular/rlimit.h>
41#include <Singular/subexpr.h>
42#include <Singular/links/silink.h>
43#include <Singular/cntrlc.h>
44#include <Singular/lists.h>
45#include <Singular/blackbox.h>
46#include <Singular/links/s_buff.h>
47#include <Singular/links/ssiLink.h>
48
49#ifdef HAVE_SIMPLEIPC
50#include <Singular/links/simpleipc.h>
51#endif
52
53#include <Singular/si_signals.h>
54
55#define SSI_VERSION 5
56
57// 64 bit version:
58//#if SIZEOF_LONG == 8
59#if 0
60#define MAX_NUM_SIZE 60
61#define POW_2_28 (1L<<60)
62#define LONG long
63#else
64// 32 bit version:
65#define MAX_NUM_SIZE 28
66#define POW_2_28 (1L<<28)
67#define LONG int
68#endif
69typedef struct
70{
71  s_buff f_read;
72  FILE *f_write;
73  ring r;
74  pid_t pid; /* only valid for fork/tcp mode*/
75  int fd_read,fd_write; /* only valid for fork/tcp mode*/
76  char level;
77  char send_quit_at_exit;
78  char quit_sent;
79
80} ssiInfo;
81
82link_list ssiToBeClosed=NULL;
83volatile BOOLEAN ssiToBeClosed_inactive=TRUE;
84
85// the helper functions:
86void ssiSetCurrRing(const ring r)
87{
88  if (!rEqual(r,currRing,1))
89  {
90    char name[20];
91    int nr=0;
92    do
93    { sprintf(name,"ssiRing%d",nr); nr++; }
94    while(IDROOT->get(name, 0)!=NULL);
95    idhdl h=enterid(omStrDup(name),0,RING_CMD,&IDROOT,FALSE);
96    IDRING(h)=r;
97    r->ref++;
98    rSetHdl(h);
99  }
100}
101// the implementation of the functions:
102void ssiWriteInt(ssiInfo *d,const int i)
103{
104  fprintf(d->f_write,"%d ",i);
105  //if (d->f_debug!=NULL) fprintf(d->f_debug,"int: %d ",i);
106}
107
108void ssiWriteString(ssiInfo *d,const char *s)
109{
110  fprintf(d->f_write,"%d %s ",(int)strlen(s),s);
111  //if (d->f_debug!=NULL) fprintf(d->f_debug,"stringi: %d \"%s\" ",strlen(s),s);
112}
113
114
115void ssiWriteBigInt(const ssiInfo *d, const number n)
116{
117  // syntax is as follows:
118  // case 2 Q:     3 4 <int>
119  //        or     3 3 <mpz_t nominator>
120  if(SR_HDL(n) & SR_INT)
121  {
122    fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
123    //if (d->f_debug!=NULL) fprintf(d->f_debug,"bigint: short \"%ld\" ",SR_TO_INT(n));
124  }
125  else if (n->s==3)
126  {
127    fputs("3 ",d->f_write);
128    mpz_out_str(d->f_write,10,n->z);
129    fputc(' ',d->f_write);
130    //gmp_fprintf(d->f_write,"3 %Zd ",n->z);
131    //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"bigint: gmp \"%Zd\" ",n->z);
132  }
133  else WerrorS("illiegal bigint");
134}
135
136void ssiWriteNumber(const ssiInfo *d, const number n)
137{
138  // syntax is as follows:
139  // case 1 Z/p:   3 <int>
140  // case 2 Q:     3 4 <int>
141  //        or     3 0 <mpz_t nominator> <mpz_t denominator>
142  //        or     3 1  dto.
143  //        or     3 3 <mpz_t nominator>
144  //        or     3 5 <mpz_t raw nom.> <mpz_t raw denom.>
145  //        or     3 6 <mpz_t raw nom.> <mpz_t raw denom.>
146  //        or     3 7 <mpz_t raw nom.>
147  if(rField_is_Zp(d->r))
148  {
149    fprintf(d->f_write,"%d ",(int)(long)n);
150    //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: \"%ld\" ",(int)(long)n);
151  }
152  else if (rField_is_Q(d->r))
153  {
154    if(SR_HDL(n) & SR_INT)
155    {
156      #if SIZEOF_LONG == 4
157      fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
158      #else
159      long nn=SR_TO_INT(n);
160      if ((nn<POW_2_28)||(nn>= -POW_2_28))
161        fprintf(d->f_write,"4 %ld ",nn);
162      else
163      {
164        mpz_t tmp;
165        mpz_init_set_si(tmp,nn);
166        fputs("8 ",d->f_write);
167        mpz_out_str (d->f_write,32, tmp);
168        fputc(' ',d->f_write);
169        mpz_clear(tmp);
170      }
171      #endif
172      //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: short \"%ld\" ",SR_TO_INT(n));
173    }
174    else if (n->s<2)
175    {
176      //gmp_fprintf(d->f_write,"%d %Zd %Zd ",n->s,n->z,n->n);
177      fprintf(d->f_write,"%d ",n->s+5);
178      mpz_out_str (d->f_write,32, n->z);
179      fputc(' ',d->f_write);
180      mpz_out_str (d->f_write,32, n->n);
181      fputc(' ',d->f_write);
182
183      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: s=%d gmp/gmp \"%Zd %Zd\" ",n->s,n->z,n->n);
184    }
185    else /*n->s==3*/
186    {
187      //gmp_fprintf(d->f_write,"3 %Zd ",n->z);
188      fputs("8 ",d->f_write);
189      mpz_out_str (d->f_write,32, n->z);
190      fputc(' ',d->f_write);
191
192      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: gmp \"%Zd\" ",n->z);
193    }
194  }
195  else WerrorS("coeff field not implemented");
196}
197
198void ssiWriteRing(ssiInfo *d,const ring r)
199{
200  /* 5 <ch> <N> <l1> <v1> ...<lN> <vN> <number of orderings> <ord1> <block0_1> <block1_1> .... */
201  if (d->r!=NULL) rKill(d->r);
202  d->r=r;
203  if (r!=NULL)
204  {
205    /*d->*/r->ref++;
206    fprintf(d->f_write,"%d %d ",n_GetChar(r->cf),r->N);
207
208    int i;
209    for(i=0;i<r->N;i++)
210    {
211      fprintf(d->f_write,"%d %s ",(int)strlen(r->names[i]),r->names[i]);
212    }
213    /* number of orderings:*/
214    i=0;
215    // remember dummy ring: everything 0:
216    if (r->order!=NULL) while (r->order[i]!=0) i++;
217    fprintf(d->f_write,"%d ",i);
218    /* each ordering block: */
219    i=0;
220    if (r->order!=NULL) while(r->order[i]!=0)
221    {
222      fprintf(d->f_write,"%d %d %d ",r->order[i],r->block0[i], r->block1[i]);
223      switch(r->order[i])
224      {
225        case ringorder_a:
226        case ringorder_wp:
227        case ringorder_Wp:
228        case ringorder_ws:
229        case ringorder_Ws:
230        case ringorder_aa:
231        {
232          int ii;
233          for(ii=r->block0[i];ii<=r->block1[i];ii++)
234            fprintf(d->f_write,"%d ",r->wvhdl[i][ii-r->block0[i]]);
235        }
236        break;
237
238        case ringorder_a64:
239        case ringorder_M:
240        case ringorder_L:
241        case ringorder_IS:
242          Werror("ring oder not implemented for ssi:%d",r->order[i]);
243          break;
244
245        default: break;
246      }
247      i++;
248    }
249  }
250  else /* dummy ring r==NULL*/
251  {
252    fprintf(d->f_write,"0 0 0 "/*,r->ch,r->N, blocks*/);
253  }
254}
255
256void ssiWritePoly(ssiInfo *d, int typ, poly p)
257{
258  fprintf(d->f_write,"%d ",pLength(p));//number of terms
259  int i;
260
261  while(p!=NULL)
262  {
263    ssiWriteNumber(d,pGetCoeff(p));
264    //nWrite(fich,pGetCoeff(p));
265    fprintf(d->f_write,"%ld ",p_GetComp(p,d->r));//component
266
267    for(int j=1;j<=rVar(d->r);j++)
268    {
269      fprintf(d->f_write,"%ld ",p_GetExp(p,j,d->r ));//x^j
270    }
271    pIter(p);
272  }
273}
274
275void ssiWriteIdeal(ssiInfo *d, int typ,ideal I)
276{
277   // syntax: 7 # of elements <poly 1> <poly2>.....
278   // syntax: 8 <rows> <cols> <poly 1> <poly2>.....
279   matrix M=(matrix)I;
280   int mn;
281   if (typ==MATRIX_CMD)
282   {
283     mn=MATROWS(M)*MATCOLS(M);
284     fprintf(d->f_write,"%d %d ", MATROWS(M),MATCOLS(M));
285   }
286   else
287   {
288     mn=IDELEMS(I);
289     fprintf(d->f_write,"%d ",IDELEMS(I));
290   }
291
292   int i;
293   int tt;
294   if (typ==MODUL_CMD) tt=VECTOR_CMD;
295   else                tt=POLY_CMD;
296
297   for(i=0;i<mn;i++)
298   {
299     ssiWritePoly(d,tt,I->m[i]);
300   }
301}
302
303void ssiWriteCommand(si_link l, command D)
304{
305  ssiInfo *d=(ssiInfo*)l->data;
306  // syntax: <num ops> <operation> <op1> <op2> ....
307  fprintf(d->f_write,"%d %d ",D->argc,D->op);
308  if (D->argc >0) ssiWrite(l, &(D->arg1));
309  if (D->argc < 4)
310  {
311    if (D->argc >1) ssiWrite(l, &(D->arg2));
312    if (D->argc >2) ssiWrite(l, &(D->arg3));
313  }
314}
315
316void ssiWriteProc(ssiInfo *d,procinfov p)
317{
318  if (p->data.s.body==NULL)
319    iiGetLibProcBuffer(p);
320  if (p->data.s.body!=NULL)
321    ssiWriteString(d,p->data.s.body);
322  else
323    ssiWriteString(d,"");
324}
325
326void ssiWriteList(si_link l,lists dd)
327{
328  ssiInfo *d=(ssiInfo*)l->data;
329  int Ll=lSize(dd);
330  fprintf(d->f_write,"%d ",Ll+1);
331  int i;
332  for(i=0;i<=Ll;i++)
333  {
334    ssiWrite(l,&(dd->m[i]));
335  }
336}
337void ssiWriteIntvec(ssiInfo *d,intvec * v)
338{
339  fprintf(d->f_write,"%d ",v->length());
340  int i;
341  for(i=0;i<v->length();i++)
342  {
343    fprintf(d->f_write,"%d ",(*v)[i]);
344  }
345}
346void ssiWriteIntmat(ssiInfo *d,intvec * v)
347{
348  fprintf(d->f_write,"%d %d ",v->rows(),v->cols());
349  int i;
350  for(i=0;i<v->length();i++)
351  {
352    fprintf(d->f_write,"%d ",(*v)[i]);
353  }
354}
355
356void ssiWriteBigintmat(ssiInfo *d,bigintmat * v)
357{
358  fprintf(d->f_write,"%d %d ",v->rows(),v->cols());
359  int i;
360  for(i=0;i<v->length();i++)
361  {
362    ssiWriteBigInt(d,(*v)[i]);
363  }
364}
365
366char *ssiReadString(ssiInfo *d)
367{
368  char *buf;
369  int l;
370  l=s_readint(d->f_read);
371  buf=(char*)omAlloc0(l+1);
372  int c =s_getc(d->f_read); /* skip ' '*/
373  int ll=s_readbytes(buf,l,d->f_read);
374  //if (ll!=l) printf("want %d, got %d bytes\n",l,ll);
375  buf[l]='\0';
376  return buf;
377}
378
379int ssiReadInt(s_buff fich)
380{
381  return s_readint(fich);
382}
383
384number ssiReadBigInt(ssiInfo *d)
385{
386  int sub_type=-1;
387  sub_type=s_readint(d->f_read);
388  switch(sub_type)
389  {
390   case 3:
391     {// read int or mpz_t or mpz_t, mpz_t
392       number n=nlRInit(0);
393       s_readmpz(d->f_read,n->z);
394       n->s=sub_type;
395       return n;
396     }
397   case 4:
398     {
399       int dd;
400       dd=s_readint(d->f_read);
401       return INT_TO_SR(dd);
402     }
403   default:
404       Werror("error in reading bigint: invalid subtype %d",sub_type);
405       return NULL;
406   }
407}
408
409static number ssiReadQNumber(ssiInfo *d)
410{
411  int sub_type=-1;
412  sub_type=s_readint(d->f_read);
413  switch(sub_type)
414  {
415     case 0:
416     case 1:
417       {// read mpz_t, mpz_t
418         number n=nlRInit(0);
419         mpz_init(n->n);
420         s_readmpz(d->f_read,n->z);
421         s_readmpz(d->f_read,n->n);
422         n->s=sub_type;
423         return n;
424       }
425
426     case 3:
427       {// read mpz_t
428         number n=nlRInit(0);
429         s_readmpz(d->f_read,n->z);
430         n->s=3; /*sub_type*/
431         return n;
432       }
433     case 4:
434       {
435         LONG dd=s_readlong(d->f_read);
436         //#if SIZEOF_LONG == 8
437         return INT_TO_SR(dd);
438         //#else
439         //return nlInit(dd,NULL);
440         //#endif
441       }
442     case 5:
443     case 6:
444       {// read raw mpz_t, mpz_t
445         number n=nlRInit(0);
446         mpz_init(n->n);
447         s_readmpz_base (d->f_read,n->z, 32);
448         s_readmpz_base (d->f_read,n->n, 32);
449         n->s=sub_type-5;
450         return n;
451       }
452     case 8:
453       {// read raw mpz_t
454         number n=nlRInit(0);
455         s_readmpz_base (d->f_read,n->z, 32);
456         n->s=sub_type=3; /*subtype-5*/
457         return n;
458       }
459
460     default: Werror("error in reading number: invalid subtype %d",sub_type);
461              return NULL;
462  }
463  return NULL;
464}
465number ssiReadNumber(ssiInfo *d)
466{
467  if (rField_is_Q(d->r))
468  {
469     return ssiReadQNumber(d);
470  }
471  else if (rField_is_Zp(d->r))
472  {
473    // read int
474    int dd;
475    dd=s_readint(d->f_read);
476    return (number)(long)dd;
477  }
478  else Werror("coeffs not implemented");
479  return NULL;
480}
481
482ring ssiReadRing(ssiInfo *d)
483{
484/* syntax is <ch> <N> <l1> <v1> ...<lN> <vN> <number of orderings> <ord1> <block0_1> <block1_1> .... */
485  int ch, N,i,l;
486  char **names;
487  ch=s_readint(d->f_read);
488  N=s_readint(d->f_read);
489  if (N!=0)
490  {
491    names=(char**)omAlloc(N*sizeof(char*));
492    for(i=0;i<N;i++)
493    {
494      names[i]=ssiReadString(d);
495    }
496  }
497  // read the orderings:
498  int num_ord; // number of orderings
499  num_ord=s_readint(d->f_read);
500  int *ord=(int *)omAlloc0((num_ord+1)*sizeof(int));
501  int *block0=(int *)omAlloc0((num_ord+1)*sizeof(int));
502  int *block1=(int *)omAlloc0((num_ord+1)*sizeof(int));
503  int **wvhdl=(int**)omAlloc0((num_ord+1)*sizeof(int*));
504  for(i=0;i<num_ord;i++)
505  {
506    ord[i]=s_readint(d->f_read);
507    block0[i]=s_readint(d->f_read);
508    block1[i]=s_readint(d->f_read);
509    switch(ord[i])
510    {
511      case ringorder_a:
512      case ringorder_wp:
513      case ringorder_Wp:
514      case ringorder_ws:
515      case ringorder_Ws:
516      case ringorder_aa:
517      {
518        wvhdl[i]=(int*)omAlloc((block1[i]-block0[i]+1)*sizeof(int));
519        int ii;
520        for(ii=block0[i];ii<=block1[i];ii++)
521          wvhdl[i][ii-block0[i]]=s_readint(d->f_read);
522      }
523      break;
524
525      case ringorder_a64:
526      case ringorder_M:
527      case ringorder_L:
528      case ringorder_IS:
529        Werror("ring oder not implemented for ssi:%d",ord[i]);
530        break;
531
532      default: break;
533    }
534  }
535  if (N==0)
536  {
537    omFree(ord);
538    omFree(block0);
539    omFree(block1);
540    omFree(wvhdl);
541    return NULL;
542  }
543  else return rDefault(ch,N,names,num_ord,ord,block0,block1,wvhdl);
544
545  return rDefault(ch,N,names,num_ord,ord,block0,block1,wvhdl);
546}
547
548poly ssiReadPoly(ssiInfo *D)
549{
550// < # of terms> < term1> < .....
551  int n,i,l;
552  n=ssiReadInt(D->f_read);
553  //Print("poly: terms:%d\n",n);
554  poly p;
555  int j;
556  j=0;
557  poly ret=NULL;
558  poly prev=NULL;
559  for(l=0;l<n;l++) // read n terms
560  {
561// coef,comp.exp1,..exp N
562    p=p_Init(D->r);
563    pGetCoeff(p)=ssiReadNumber(D);
564    int d;
565    d=s_readint(D->f_read);
566    p_SetComp(p,d,D->r);
567    for(i=1;i<=rVar(D->r);i++)
568    {
569      d=s_readint(D->f_read);
570      p_SetExp(p,i,d,D->r);
571    }
572    p_Setm(p,D->r);
573    p_Test(p,D->r);
574    if (ret==NULL) ret=p;
575    else           pNext(prev)=p;
576    prev=p;
577 }
578 return ret;
579}
580
581ideal ssiReadIdeal(ssiInfo *d)
582{
583  int n,i;
584  ideal I;
585  n=s_readint(d->f_read);
586  I=idInit(n,1);
587  for(i=0;i<IDELEMS(I);i++) // read n terms
588  {
589    I->m [i]=ssiReadPoly(d);
590  }
591  return I;
592}
593
594matrix ssiReadMatrix(ssiInfo *d)
595{
596  int n,m,i,j;
597  m=s_readint(d->f_read);
598  n=s_readint(d->f_read);
599  matrix M=mpNew(m,n);
600  poly p;
601  for(int i=1;i<=MATROWS(M);i++)
602    for(int j=1;j<=MATCOLS(M);j++)
603    {
604      p=ssiReadPoly(d);
605      MATELEM(M,i,j)=p;
606    }
607  return M;
608}
609
610command ssiReadCommand(si_link l)
611{
612  ssiInfo *d=(ssiInfo*)l->data;
613  // syntax: <num ops> <operation> <op1> <op2> ....
614  command D=(command)omAlloc0(sizeof(*D));
615  int argc,op;
616  argc=s_readint(d->f_read);
617  op=s_readint(d->f_read);
618  D->argc=argc; D->op=op;
619  leftv v;
620  if (argc >0)
621  {
622    v=ssiRead1(l);
623    memcpy(&(D->arg1),v,sizeof(*v));
624    omFreeBin(v,sleftv_bin);
625  }
626  if (argc <4)
627  {
628    if (D->argc >1)
629    {
630      v=ssiRead1(l);
631      memcpy(&(D->arg2),v,sizeof(*v));
632      omFreeBin(v,sleftv_bin);
633    }
634    if (D->argc >2)
635    {
636      v=ssiRead1(l);
637      memcpy(&(D->arg3),v,sizeof(*v));
638      omFreeBin(v,sleftv_bin);
639    }
640  }
641  else
642  {
643    leftv prev=&(D->arg1);
644    argc--;
645    while(argc >0)
646    {
647      v=ssiRead1(l);
648      prev->next=v;
649      prev=v;
650      argc--;
651    }
652  }
653  return D;
654}
655
656procinfov ssiReadProc(ssiInfo *d)
657{
658  char *s=ssiReadString(d);
659  procinfov p=(procinfov)omAlloc0Bin(procinfo_bin);
660  p->language=LANG_SINGULAR;
661  p->libname=omStrDup("");
662  p->procname=omStrDup("");
663  p->data.s.body=s;
664  return p;
665}
666lists ssiReadList(si_link l)
667{
668  ssiInfo *d=(ssiInfo*)l->data;
669  int nr;
670  nr=s_readint(d->f_read);
671  lists L=(lists)omAlloc(sizeof(*L));
672  L->Init(nr);
673
674  int i;
675  leftv v;
676  for(i=0;i<nr;i++)
677  {
678    v=ssiRead1(l);
679    memcpy(&(L->m[i]),v,sizeof(*v));
680    omFreeBin(v,sleftv_bin);
681  }
682  return L;
683}
684intvec* ssiReadIntvec(ssiInfo *d)
685{
686  int nr;
687  nr=s_readint(d->f_read);
688  intvec *v=new intvec(nr);
689  for(int i=0;i<nr;i++)
690  {
691    (*v)[i]=s_readint(d->f_read);
692  }
693  return v;
694}
695intvec* ssiReadIntmat(ssiInfo *d)
696{
697  int r,c;
698  r=s_readint(d->f_read);
699  c=s_readint(d->f_read);
700  intvec *v=new intvec(r,c,0);
701  for(int i=0;i<r*c;i++)
702  {
703    (*v)[i]=s_readint(d->f_read);
704  }
705  return v;
706}
707bigintmat* ssiReadBigintmat(ssiInfo *d)
708{
709  int r,c;
710  r=s_readint(d->f_read);
711  c=s_readint(d->f_read);
712  bigintmat *v=new bigintmat(r,c,coeffs_BIGINT);
713  for(int i=0;i<r*c;i++)
714  {
715    (*v)[i]=ssiReadBigInt(d);
716  }
717  return v;
718}
719
720void ssiReadBlackbox(leftv res, si_link l)
721{
722  ssiInfo *d=(ssiInfo*)l->data;
723  int throwaway;
724  throwaway=s_readint(d->f_read);
725  char *name=ssiReadString(d);
726  int tok;
727  blackboxIsCmd(name,tok);
728  if (tok>MAX_TOK)
729  {
730    blackbox *b=getBlackboxStuff(tok);
731    res->rtyp=tok;
732    b->blackbox_deserialize(&b,&(res->data),l);
733  }
734  else
735  {
736    Werror("blackbox %s not found",name);
737  }
738}
739
740//**************************************************************************/
741
742BOOLEAN ssiOpen(si_link l, short flag, leftv u)
743{
744  if (l!=NULL)
745  {
746    const char *mode;
747    ssiInfo *d=(ssiInfo*)omAlloc0(sizeof(ssiInfo));
748    if (flag & SI_LINK_OPEN)
749    {
750      if (l->mode[0] != '\0' && (strcmp(l->mode, "r") == 0))
751        flag = SI_LINK_READ;
752      else flag = SI_LINK_WRITE;
753    }
754
755    if (flag == SI_LINK_READ) mode = "r";
756    else if (strcmp(l->mode, "w") == 0) mode = "w";
757    else if (strcmp(l->mode, "fork") == 0) mode = "fork";
758    else if (strcmp(l->mode, "tcp") == 0) mode = "tcp";
759    else if (strcmp(l->mode, "connect") == 0) mode = "connect";
760    else mode = "a";
761
762
763    SI_LINK_SET_OPEN_P(l, flag);
764    l->data=d;
765    omFree(l->mode);
766    l->mode = omStrDup(mode);
767
768    if (l->name[0] == '\0')
769    {
770      if (strcmp(mode,"fork")==0)
771      {
772        link_list n=(link_list)omAlloc(sizeof(link_struct));
773        n->u=u;
774        n->l=l;
775        n->next=(void *)ssiToBeClosed;
776        ssiToBeClosed=n;
777
778        int pc[2];
779        int cp[2];
780        pipe(pc);
781        pipe(cp);
782        pid_t pid = fork();
783        if (pid == -1 && errno == EAGAIN)   // RLIMIT_NPROC too low?
784        {
785          raise_rlimit_nproc();
786          pid = fork();
787        }
788        if (pid == -1)
789        {
790          WerrorS("could not fork");
791        }
792        if (pid==0) /*fork: child*/
793        {
794          /* block SIGINT */
795          sigset_t sigint;
796          sigemptyset(&sigint);
797          sigaddset(&sigint, SIGINT);
798          sigprocmask(SIG_BLOCK, &sigint, NULL);
799
800          link_list hh=(link_list)ssiToBeClosed->next;
801          /* we know: l is the first entry in ssiToBeClosed-list */
802          while(hh!=NULL)
803          {
804            SI_LINK_SET_CLOSE_P(hh->l);
805            ssiInfo *dd=(ssiInfo*)hh->l->data;
806            s_close(dd->f_read);
807            s_free(dd->f_read);
808            fclose(dd->f_write);
809            if (dd->r!=NULL) rKill(dd->r);
810            omFreeSize((ADDRESS)dd,(sizeof *dd));
811            hh->l->data=NULL;
812            link_list nn=(link_list)hh->next;
813            omFree(hh);
814            hh=nn;
815          }
816          ssiToBeClosed->next=NULL;
817#ifdef HAVE_SIMPLEIPC
818          memset(sem_acquired, 0, SIPC_MAX_SEMAPHORES*sizeof(sem_acquired[0]));
819#endif   // HAVE_SIMPLEIPC
820          si_close(pc[1]); si_close(cp[0]);
821          d->f_write=fdopen(cp[1],"w");
822          d->f_read=s_open(pc[0]);
823          d->fd_read=pc[0];
824          d->fd_write=cp[1];
825          l->data=d;
826          omFree(l->mode);
827          l->mode = omStrDup(mode);
828          singular_in_batchmode=TRUE;
829          SI_LINK_SET_RW_OPEN_P(l);
830          //myynest=0;
831          fe_fgets_stdin=fe_fgets_dummy;
832          if ((u!=NULL)&&(u->rtyp==IDHDL))
833          {
834            idhdl h=(idhdl)u->data;
835            h->lev=0;
836          }
837          loop
838          {
839            leftv h=ssiRead1(l); /*contains an exit.... */
840            if (feErrors != NULL && *feErrors != '\0')
841            {
842              // handle errors:
843              PrintS(feErrors); /* currently quite simple */
844              *feErrors = '\0';
845            }
846            ssiWrite(l,h);
847            h->CleanUp();
848            omFreeBin(h, sleftv_bin);
849          }
850          /* never reached*/
851        }
852        else if (pid>0) /*fork: parent*/
853        {
854          d->pid=pid;
855          si_close(pc[0]); si_close(cp[1]);
856          d->f_write=fdopen(pc[1],"w");
857          d->f_read=s_open(cp[0]);
858          d->fd_read=cp[0];
859          d->fd_write=pc[1];
860          SI_LINK_SET_RW_OPEN_P(l);
861          d->send_quit_at_exit=1;
862        }
863        else
864        {
865          Werror("fork failed (%d)",errno);
866          l->data=NULL;
867          omFree(d);
868          return TRUE;
869        }
870      }
871      // ---------------------------------------------------------------------
872      else if (strcmp(mode,"tcp")==0)
873      {
874        int sockfd, newsockfd, portno, clilen;
875        struct sockaddr_in serv_addr, cli_addr;
876        int n;
877        sockfd = socket(AF_INET, SOCK_STREAM, 0);
878        if(sockfd < 0)
879        {
880          WerrorS("ERROR opening socket");
881          l->data=NULL;
882          omFree(d);
883          return TRUE;
884        }
885        memset((char *) &serv_addr,0, sizeof(serv_addr));
886        portno = 1025;
887        serv_addr.sin_family = AF_INET;
888        serv_addr.sin_addr.s_addr = INADDR_ANY;
889        do
890        {
891          portno++;
892          serv_addr.sin_port = htons(portno);
893          if(portno > 50000)
894          {
895            WerrorS("ERROR on binding (no free port available?)");
896            l->data=NULL;
897            omFree(d);
898            return TRUE;
899          }
900        }
901        while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
902        Print("waiting on port %d\n", portno);mflush();
903        listen(sockfd,1);
904        newsockfd = si_accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
905        if(newsockfd < 0)
906        {
907          WerrorS("ERROR on accept");
908          l->data=NULL;
909          omFree(d);
910          return TRUE;
911        }
912        PrintS("client accepted\n");
913        d->fd_read = newsockfd;
914        d->fd_write = newsockfd;
915        d->f_read = s_open(newsockfd);
916        d->f_write = fdopen(newsockfd, "w");
917        SI_LINK_SET_RW_OPEN_P(l);
918        si_close(sockfd);
919      }
920      // no ssi-Link on stdin or stdout
921      else
922      {
923        Werror("invalid mode >>%s<< for ssi",mode);
924        l->data=NULL;
925        omFree(d);
926        return TRUE;
927      }
928    }
929    // =========================================================================
930    else /*l->name=NULL*/
931    {
932      // tcp mode
933      if(strcmp(mode,"tcp")==0)
934      {
935        int sockfd, newsockfd, portno, clilen;
936        struct sockaddr_in serv_addr, cli_addr;
937        int n;
938        sockfd = socket(AF_INET, SOCK_STREAM, 0);
939        if(sockfd < 0)
940        {
941          WerrorS("ERROR opening socket");
942          l->data=NULL;
943          omFree(d);
944          return TRUE;
945        }
946        memset((char *) &serv_addr,0, sizeof(serv_addr));
947        portno = 1025;
948        serv_addr.sin_family = AF_INET;
949        serv_addr.sin_addr.s_addr = INADDR_ANY;
950        do
951        {
952          portno++;
953          serv_addr.sin_port = htons(portno);
954          if(portno > 50000)
955          {
956            WerrorS("ERROR on binding (no free port available?)");
957            l->data=NULL;
958            return TRUE;
959          }
960        }
961        while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
962        //Print("waiting on port %d\n", portno);mflush();
963        listen(sockfd,1);
964        char* cli_host = (char*)omAlloc(256);
965        char* path = (char*)omAlloc(1024);
966        int r = si_sscanf(l->name,"%255[^:]:%s",cli_host,path);
967        if(r == 0)
968        {
969          WerrorS("ERROR: no host specified");
970          l->data=NULL;
971          omFree(d);
972          omFree(path);
973          omFree(cli_host);
974          return TRUE;
975        }
976        else if(r == 1)
977        {
978          WarnS("program not specified, using /usr/local/bin/Singular");
979          strcpy(path,"/usr/local/bin/Singular");
980        }
981        char* ssh_command = (char*)omAlloc(256);
982        char* ser_host = (char*)omAlloc(64);
983        gethostname(ser_host,64);
984        sprintf(ssh_command,"ssh %s %s -q --batch --link=ssi --MPhost=%s --MPport=%d &",cli_host,path,ser_host,portno);
985        //Print("client on %s started:%s\n",cli_host,path);
986        omFree(path);
987        omFree(cli_host);
988        if (TEST_OPT_PROT) { Print("running >>%s<<\n",ssh_command); }
989        system(ssh_command);
990        omFree(ssh_command);
991        omFree(ser_host);
992        clilen = sizeof(cli_addr);
993        newsockfd = si_accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
994        if(newsockfd < 0)
995        {
996          WerrorS("ERROR on accept");
997          l->data=NULL;
998          omFree(d);
999          return TRUE;
1000        }
1001        //PrintS("client accepted\n");
1002        d->fd_read = newsockfd;
1003        d->fd_write = newsockfd;
1004        d->f_read = s_open(newsockfd);
1005        d->f_write = fdopen(newsockfd, "w");
1006        si_close(sockfd);
1007        SI_LINK_SET_RW_OPEN_P(l);
1008        d->send_quit_at_exit=1;
1009        link_list newlink=(link_list)omAlloc(sizeof(link_struct));
1010        newlink->u=u;
1011        newlink->l=l;
1012        newlink->next=(void *)ssiToBeClosed;
1013        ssiToBeClosed=newlink;
1014        fprintf(d->f_write,"98 %d %d %u %u\n",SSI_VERSION,MAX_TOK,si_opt_1,si_opt_2);
1015      }
1016      // ----------------------------------------------------------------------
1017      else if(strcmp(mode,"connect")==0)
1018      {
1019        char* host = (char*)omAlloc(256);
1020        int sockfd, portno, n;
1021        struct sockaddr_in serv_addr;
1022        struct hostent *server;
1023
1024        si_sscanf(l->name,"%255[^:]:%d",host,&portno);
1025        //Print("connect to host %s, port %d\n",host,portno);mflush();
1026        if (portno!=0)
1027        {
1028          sockfd = socket(AF_INET, SOCK_STREAM, 0);
1029          if (sockfd < 0) { WerrorS("ERROR opening socket"); return TRUE; }
1030          server = gethostbyname(host);
1031          if (server == NULL) {  WerrorS("ERROR, no such host");  return TRUE; }
1032          memset((char *) &serv_addr, 0, sizeof(serv_addr));
1033          serv_addr.sin_family = AF_INET;
1034          memcpy((char *)&serv_addr.sin_addr.s_addr,
1035                (char *)server->h_addr,
1036                server->h_length);
1037          serv_addr.sin_port = htons(portno);
1038          if (si_connect(sockfd,(sockaddr*)&serv_addr,sizeof(serv_addr)) < 0)
1039          { Werror("ERROR connecting(errno=%d)",errno); return TRUE; }
1040          //PrintS("connected\n");mflush();
1041          d->f_read=s_open(sockfd);
1042          d->fd_read=sockfd;
1043          d->f_write=fdopen(sockfd,"w");
1044          d->fd_write=sockfd;
1045          SI_LINK_SET_RW_OPEN_P(l);
1046          omFree(host);
1047        }
1048        else
1049        {
1050          l->data=NULL;
1051          omFree(d);
1052          return TRUE;
1053        }
1054      }
1055      // ======================================================================
1056      else
1057      {
1058        // normal link to a file
1059        FILE *outfile;
1060        char *filename=l->name;
1061
1062        if(filename[0]=='>')
1063        {
1064          if (filename[1]=='>')
1065          {
1066            filename+=2;
1067            mode = "a";
1068          }
1069          else
1070          {
1071            filename++;
1072            mode="w";
1073          }
1074        }
1075        outfile=myfopen(filename,mode);
1076        if (outfile!=NULL)
1077        {
1078          if (strcmp(l->mode,"r")==0)
1079          {
1080            fclose(outfile);
1081            d->f_read=s_open_by_name(filename);
1082          }
1083          else
1084          {
1085            d->f_write = outfile;
1086            fprintf(d->f_write,"98 %d %d %u %u\n",SSI_VERSION,MAX_TOK,si_opt_1,si_opt_2);
1087          }
1088        }
1089        else
1090        {
1091          omFree(d);
1092          l->data=NULL;
1093          return TRUE;
1094        }
1095      }
1096    }
1097  }
1098
1099  return FALSE;
1100}
1101
1102//**************************************************************************/
1103BOOLEAN ssiPrepClose(si_link l)
1104{
1105  if (l!=NULL)
1106  {
1107    ssiInfo *d = (ssiInfo *)l->data;
1108    if (d!=NULL)
1109    {
1110      if (d->send_quit_at_exit)
1111      {
1112        fputs("99\n",d->f_write);
1113        fflush(d->f_write);
1114      }
1115      d->quit_sent=1;
1116    }
1117  }
1118  return FALSE;
1119}
1120
1121BOOLEAN ssiClose(si_link l)
1122{
1123  if (l!=NULL)
1124  {
1125    SI_LINK_SET_CLOSE_P(l);
1126    ssiInfo *d = (ssiInfo *)l->data;
1127    if (d!=NULL)
1128    {
1129      if ((d->send_quit_at_exit)
1130      && (d->quit_sent==0))
1131      {
1132        fputs("99\n",d->f_write);
1133        fflush(d->f_write);
1134      }
1135      if (d->r!=NULL) rKill(d->r);
1136      if ((d->pid!=0)
1137      && (si_waitpid(d->pid,NULL,WNOHANG)==0))
1138      {
1139        struct timespec t;
1140        t.tv_sec=0;
1141        t.tv_nsec=100000000; // <=100 ms
1142        struct timespec rem;
1143        int r;
1144        do
1145        {
1146          r = nanosleep(&t, &rem);
1147          t = rem;
1148        } while ((r < 0) && (errno == EINTR)
1149            && (si_waitpid(d->pid,NULL,WNOHANG) == 0));
1150        if ((r == 0) && (si_waitpid(d->pid,NULL,WNOHANG) == 0))
1151        {
1152          kill(d->pid,15);
1153          t.tv_sec=5; // <=5s
1154          t.tv_nsec=0;
1155          do
1156          {
1157            r = nanosleep(&t, &rem);
1158            t = rem;
1159          } while ((r < 0) && (errno == EINTR)
1160              && (si_waitpid(d->pid,NULL,WNOHANG) == 0));
1161          if ((r == 0) && (si_waitpid(d->pid,NULL,WNOHANG) == 0))
1162          {
1163            kill(d->pid,9); // just to be sure
1164            si_waitpid(d->pid,NULL,0);
1165          }
1166        }
1167      }
1168      if (d->f_read!=NULL) s_close(d->f_read);
1169      if (d->f_read!=NULL) s_free(d->f_read);
1170      if (d->f_write!=NULL) fclose(d->f_write);
1171      if ((strcmp(l->mode,"tcp")==0)
1172      || (strcmp(l->mode,"fork")==0))
1173      {
1174        link_list hh=ssiToBeClosed;
1175        if (hh!=NULL)
1176        {
1177          if (hh->l==l)
1178          {
1179             ssiToBeClosed=(link_list)hh->next;
1180             omFreeSize(hh,sizeof(link_struct));
1181          }
1182          else while(hh->next!=NULL)
1183          {
1184            link_list hhh=(link_list)hh->next;
1185            if (hhh->l==l)
1186            {
1187              hh->next=hhh->next;
1188              omFreeSize(hhh,sizeof(link_struct));
1189              break;
1190            }
1191            else
1192              hh=(link_list)hh->next;
1193          }
1194        }
1195      }
1196      omFreeSize((ADDRESS)d,(sizeof *d));
1197    }
1198    l->data=NULL;
1199  }
1200  return FALSE;
1201}
1202
1203//**************************************************************************/
1204leftv ssiRead1(si_link l)
1205{
1206  ssiInfo *d = (ssiInfo *)l->data;
1207  leftv res=(leftv)omAlloc0(sizeof(sleftv));
1208  int t=0;
1209  t=s_readint(d->f_read);
1210  //Print("got type %d\n",t);
1211  switch(t)
1212  {
1213    case 1:res->rtyp=INT_CMD;
1214           res->data=(char *)(long)ssiReadInt(d->f_read);
1215           break;
1216    case 2:res->rtyp=STRING_CMD;
1217           res->data=(char *)ssiReadString(d);
1218           break;
1219    case 3:res->rtyp=NUMBER_CMD;
1220           res->data=(char *)ssiReadNumber(d);
1221           break;
1222    case 4:res->rtyp=BIGINT_CMD;
1223           res->data=(char *)ssiReadBigInt(d);
1224           break;
1225    case 15:
1226    case 5:{
1227             d->r=ssiReadRing(d);
1228             res->rtyp=RING_CMD;
1229             res->data=(char*)d->r;
1230             // we are in the top-level, so set the basering to d->r:
1231             if (d->r!=NULL)
1232             {
1233               d->r->ref++;
1234               ssiSetCurrRing(d->r);
1235             }
1236             if (t==15) return ssiRead1(l);
1237           }
1238           break;
1239    case 6:res->rtyp=POLY_CMD;
1240           if (d->r==NULL) goto no_ring;
1241           res->data=(char*)ssiReadPoly(d);
1242           break;
1243    case 7:res->rtyp=IDEAL_CMD;
1244           if (d->r==NULL) goto no_ring;
1245           res->data=(char*)ssiReadIdeal(d);
1246           break;
1247    case 8:res->rtyp=MATRIX_CMD;
1248           if (d->r==NULL) goto no_ring;
1249           res->data=(char*)ssiReadMatrix(d);
1250           break;
1251    case 9:res->rtyp=VECTOR_CMD;
1252           if (d->r==NULL) goto no_ring;
1253           res->data=(char*)ssiReadPoly(d);
1254           break;
1255    case 10:res->rtyp=MODUL_CMD;
1256           if (d->r==NULL) goto no_ring;
1257           res->data=(char*)ssiReadIdeal(d);
1258           break;
1259    case 11:
1260           {
1261             res->rtyp=COMMAND;
1262             res->data=ssiReadCommand(l);
1263             int nok=res->Eval();
1264             if (nok) WerrorS("error in eval");
1265             break;
1266           }
1267    case 12: /*DEF_CMD*/
1268           {
1269             res->rtyp=0;
1270             res->name=(char *)ssiReadString(d);
1271             int nok=res->Eval();
1272             if (nok) WerrorS("error in name lookup");
1273             break;
1274           }
1275    case 13: res->rtyp=PROC_CMD;
1276             res->data=ssiReadProc(d);
1277             break;
1278    case 14: res->rtyp=LIST_CMD;
1279             res->data=ssiReadList(l);
1280             break;
1281    case 16: res->rtyp=NONE; res->data=NULL;
1282             break;
1283    case 17: res->rtyp=INTVEC_CMD;
1284             res->data=ssiReadIntvec(d);
1285             break;
1286    case 18: res->rtyp=INTMAT_CMD;
1287             res->data=ssiReadIntmat(d);
1288             break;
1289    case 19: res->rtyp=BIGINTMAT_CMD;
1290             res->data=ssiReadBigintmat(d);
1291             break;
1292    case 20: ssiReadBlackbox(res,l);
1293             break;
1294    // ------------
1295    case 98: // version
1296             {
1297                int n98_v,n98_m;
1298                BITSET n98_o1,n98_o2;
1299                n98_v=s_readint(d->f_read);
1300                n98_m=s_readint(d->f_read);
1301                n98_o1=s_readint(d->f_read);
1302                n98_o2=s_readint(d->f_read);
1303                if ((n98_v!=SSI_VERSION) ||(n98_m!=MAX_TOK))
1304                {
1305                  Print("incompatible versions of ssi: %d/%d vs %d/%d",
1306                                  SSI_VERSION,MAX_TOK,n98_v,n98_m);
1307                }
1308                #ifndef NDEBUG
1309                if (TEST_OPT_DEBUG)
1310                  Print("// opening ssi-%d, MAX_TOK=%d\n",n98_v,n98_m);
1311                #endif
1312                si_opt_1=n98_o1;
1313                si_opt_2=n98_o2;
1314                return ssiRead1(l);
1315             }
1316    case 99: ssiClose(l); m2_end(0);
1317    case 0: if (s_iseof(d->f_read))
1318            {
1319              ssiClose(l);
1320              res->rtyp=DEF_CMD;
1321              break;
1322            }
1323    default: Werror("not implemented (t:%d)",t);
1324             omFreeSize(res,sizeof(sleftv));
1325             res=NULL;
1326             break;
1327  }
1328  return res;
1329no_ring: WerrorS("no ring");
1330  omFreeSize(res,sizeof(sleftv));
1331  return NULL;
1332}
1333//**************************************************************************/
1334BOOLEAN ssiWrite(si_link l, leftv data)
1335{
1336  if(SI_LINK_W_OPEN_P(l)==0)
1337     if (slOpen(l,SI_LINK_OPEN|SI_LINK_WRITE,NULL)) return TRUE;
1338  ssiInfo *d = (ssiInfo *)l->data;
1339  d->level++;
1340  //FILE *fich=d->f;
1341  while (data!=NULL)
1342  {
1343    int tt=data->Typ();
1344    void *dd=data->Data();
1345    if ((dd==NULL) && (data->name!=NULL) && (tt==0)) tt=DEF_CMD;
1346      // return pure undefined names as def
1347
1348    switch(tt /*data->Typ()*/)
1349    {
1350          case NONE/* nothing*/:fputs("16 ",d->f_write);
1351                          break;
1352          case STRING_CMD: fputs("2 ",d->f_write);
1353                           ssiWriteString(d,(char *)dd);
1354                           break;
1355          case INT_CMD: fputs("1 ",d->f_write);
1356                        ssiWriteInt(d,(int)(long)dd);
1357                        break;
1358          case BIGINT_CMD:fputs("4 ",d->f_write);
1359                        ssiWriteBigInt(d,(number)dd);
1360                        break;
1361          case NUMBER_CMD:
1362                          if (d->r!=currRing)
1363                          {
1364                            fputs("15 ",d->f_write);
1365                            ssiWriteRing(d,currRing);
1366                            if (d->level<=1) fputc('\n',d->f_write);
1367                          }
1368                          fputs("3 ",d->f_write);
1369                          ssiWriteNumber(d,(number)dd);
1370                        break;
1371          case RING_CMD:fputs("5 ",d->f_write);
1372                        ssiWriteRing(d,(ring)dd);
1373                        break;
1374          case POLY_CMD:
1375          case VECTOR_CMD:
1376                        if (d->r!=currRing)
1377                        {
1378                          fputs("15 ",d->f_write);
1379                          ssiWriteRing(d,currRing);
1380                          if (d->level<=1) fputc('\n',d->f_write);
1381                        }
1382                        if(tt==POLY_CMD) fputs("6 ",d->f_write);
1383                        else             fputs("9 ",d->f_write);
1384                        ssiWritePoly(d,tt,(poly)dd);
1385                        break;
1386          case IDEAL_CMD:
1387          case MODUL_CMD:
1388          case MATRIX_CMD:
1389                        if (d->r!=currRing)
1390                        {
1391                          fputs("15 ",d->f_write);
1392                          ssiWriteRing(d,currRing);
1393                          if (d->level<=1) fputc('\n',d->f_write);
1394                        }
1395                        if(tt==IDEAL_CMD)       fputs("7 ",d->f_write);
1396                        else if(tt==MATRIX_CMD) fputs("8 ",d->f_write);
1397                        else                    fputs("10 ",d->f_write);
1398                        ssiWriteIdeal(d,tt,(ideal)dd);
1399                        break;
1400          case COMMAND:
1401                   fputs("11 ",d->f_write);
1402                   ssiWriteCommand(l,(command)dd);
1403                   break;
1404          case DEF_CMD: /* not evaluated stuff in quotes */
1405                   fputs("12 ",d->f_write);
1406                   ssiWriteString(d,data->Name());
1407                   break;
1408          case PROC_CMD:
1409                   fputs("13 ",d->f_write);
1410                   ssiWriteProc(d,(procinfov)dd);
1411                   break;
1412          case LIST_CMD:
1413                   fputs("14 ",d->f_write);
1414                   ssiWriteList(l,(lists)dd);
1415                   break;
1416          case INTVEC_CMD:
1417                   fputs("17 ",d->f_write);
1418                   ssiWriteIntvec(d,(intvec *)dd);
1419                   break;
1420          case INTMAT_CMD:
1421                   fputs("18 ",d->f_write);
1422                   ssiWriteIntmat(d,(intvec *)dd);
1423                   break;
1424          case BIGINTMAT_CMD:
1425                   fputs("19 ",d->f_write);
1426                   ssiWriteBigintmat(d,(bigintmat *)dd);
1427                   break;
1428          default:
1429            if (tt>MAX_TOK)
1430            {
1431              blackbox *b=getBlackboxStuff(tt);
1432              fputs("20 ",d->f_write);
1433              b->blackbox_serialize(b,dd,l);
1434            }
1435            else
1436            {
1437              Werror("not implemented (t:%d, rtyp:%d)",tt, data->rtyp);
1438              d->level=0;
1439              return TRUE;
1440            }
1441            break;
1442    }
1443    if (d->level<=1) { fputc('\n',d->f_write); fflush(d->f_write); }
1444    data=data->next;
1445  }
1446  d->level--;
1447  return FALSE;
1448}
1449
1450BOOLEAN ssiGetDump(si_link l);
1451BOOLEAN ssiDump(si_link l);
1452
1453si_link_extension slInitSsiExtension(si_link_extension s)
1454{
1455  s->Open=ssiOpen;
1456  s->Close=ssiClose;
1457  s->Kill=ssiClose;
1458  s->Read=ssiRead1;
1459  s->Read2=(slRead2Proc)NULL;
1460  s->Write=ssiWrite;
1461  s->Dump=ssiDump;
1462  s->GetDump=ssiGetDump;
1463
1464  s->Status=slStatusSsi;
1465  s->type="ssi";
1466  return s;
1467}
1468
1469const char* slStatusSsi(si_link l, const char* request)
1470{
1471  ssiInfo *d=(ssiInfo*)l->data;
1472  if (d==NULL) return "not open";
1473  if (((strcmp(l->mode,"fork")==0)
1474  ||(strcmp(l->mode,"tcp")==0)
1475  ||(strcmp(l->mode,"connect")==0))
1476  && (strcmp(request, "read") == 0))
1477  {
1478    fd_set  mask, fdmask;
1479    struct timeval wt;
1480    if (s_isready(d->f_read)) return "ready";
1481    loop
1482    {
1483      /* Don't block. Return socket status immediately. */
1484      wt.tv_sec  = 0;
1485      wt.tv_usec = 0;
1486
1487      FD_ZERO(&mask);
1488      FD_SET(d->fd_read, &mask);
1489      //Print("test fd %d\n",d->fd_read);
1490    /* check with select: chars waiting: no -> not ready */
1491      switch (si_select(d->fd_read+1, &mask, NULL, NULL, &wt))
1492      {
1493        case 0: /* not ready */ return "not ready";
1494        case -1: /*error*/      return "error";
1495        case 1: /*ready ? */    break;
1496      }
1497    /* yes: read 1 char*/
1498    /* if \n, check again with select else ungetc(c), ready*/
1499      int c=s_getc(d->f_read);
1500      //Print("try c=%d\n",c);
1501      if (c== -1) return "eof"; /* eof or error */
1502      else if (isdigit(c))
1503      { s_ungetc(c,d->f_read); return "ready"; }
1504      else if (c>' ')
1505      {
1506        Werror("unknown char in ssiLink(%d)",c);
1507        return "error";
1508      }
1509      /* else: next char */
1510    }
1511  }
1512  else if (strcmp(request, "read") == 0)
1513  {
1514    if (SI_LINK_R_OPEN_P(l) && (!s_iseof(d->f_read)) && (s_isready(d->f_read))) return "ready";
1515    else return "not ready";
1516  }
1517  else if (strcmp(request, "write") == 0)
1518  {
1519    if (SI_LINK_W_OPEN_P(l)) return "ready";
1520    else return "not ready";
1521  }
1522  else return "unknown status request";
1523}
1524
1525int slStatusSsiL(lists L, int timeout)
1526{
1527// input: L: a list with links of type
1528//           ssi-connect, ssi-fork, ssi-tcp, MPtcp-fork or MPtcp-launch.
1529//           Note: Not every entry in L must be set.
1530//        timeout: timeout for select in micro-seconds
1531//           or -1 for infinity
1532//           or 0 for polling
1533// returns: ERROR (via Werror): L has wrong elements or link not open
1534//           -2: select returns an error
1535//           -1: the read state of all links is eof
1536//           0:  timeout (or polling): none ready,
1537//           i>0: (at least) L[i] is ready
1538  si_link l;
1539  ssiInfo *d;
1540  int d_fd;
1541  fd_set  mask, fdmask;
1542  FD_ZERO(&fdmask);
1543  FD_ZERO(&mask);
1544  int max_fd=0; /* 1 + max fd in fd_set */
1545
1546  /* timeout */
1547  struct timeval wt;
1548  struct timeval *wt_ptr=&wt;
1549  int startingtime = getRTimer()/TIMER_RESOLUTION;  // in seconds
1550  if (timeout== -1)
1551  {
1552    wt_ptr=NULL;
1553  }
1554  else
1555  {
1556    wt.tv_sec  = timeout / 1000000;
1557    wt.tv_usec = timeout % 1000000;
1558  }
1559
1560  /* auxiliary variables */
1561  int i;
1562  int j;
1563  int k;
1564  int s;
1565  char fdmaskempty;
1566
1567  /* check the links and fill in fdmask */
1568  /* check ssi links for ungetc_buf */
1569  for(i=L->nr; i>=0; i--)
1570  {
1571    if (L->m[i].Typ()!=DEF_CMD)
1572    {
1573      if (L->m[i].Typ()!=LINK_CMD)
1574      { WerrorS("all elements must be of type link"); return -2;}
1575      l=(si_link)L->m[i].Data();
1576      if(SI_LINK_OPEN_P(l)==0)
1577      { WerrorS("all links must be open"); return -2;}
1578      if (((strcmp(l->m->type,"ssi")!=0) && (strcmp(l->m->type,"MPtcp")!=0))
1579      || ((strcmp(l->mode,"fork")!=0) && (strcmp(l->mode,"tcp")!=0)
1580        && (strcmp(l->mode,"launch")!=0) && (strcmp(l->mode,"connect")!=0)))
1581      {
1582        WerrorS("all links must be of type ssi:fork, ssi:tcp, ssi:connect");
1583        return -2;
1584      }
1585      if (strcmp(l->m->type,"ssi")==0)
1586      {
1587        d=(ssiInfo*)l->data;
1588        d_fd=d->fd_read;
1589        if (!s_isready(d->f_read))
1590        {
1591          FD_SET(d_fd, &fdmask);
1592          if (d_fd > max_fd) max_fd=d_fd;
1593        }
1594        else
1595          return i+1;
1596      }
1597      else
1598      {
1599        Werror("wrong link type >>%s<<",l->m->type);
1600        return -2;
1601      }
1602    }
1603  }
1604  max_fd++;
1605
1606do_select:
1607  /* copy fdmask to mask */
1608  FD_ZERO(&mask);
1609  for(k = 0; k < max_fd; k++)
1610  {
1611    if(FD_ISSET(k, &fdmask))
1612    {
1613      FD_SET(k, &mask);
1614    }
1615  }
1616
1617  /* check with select: chars waiting: no -> not ready */
1618  s = si_select(max_fd, &mask, NULL, NULL, wt_ptr);
1619  if (s==-1)
1620  {
1621    WerrorS("error in select call");
1622    return -2; /*error*/
1623  }
1624  if (s==0)
1625  {
1626    return 0; /*poll: not ready */
1627  }
1628  else /* s>0, at least one ready  (the number of fd which are ready is s)*/
1629  {
1630    j=0;
1631    while (j<=max_fd) { if (FD_ISSET(j,&mask)) break; j++; }
1632    for(i=L->nr; i>=0; i--)
1633    {
1634      if (L->m[i].rtyp==LINK_CMD)
1635      {
1636        l=(si_link)L->m[i].Data();
1637        if (strcmp(l->m->type,"ssi")==0)
1638        {
1639          d=(ssiInfo*)l->data;
1640          d_fd=d->fd_read;
1641          if(j==d_fd) break;
1642        }
1643        else
1644        {
1645          Werror("wrong link type >>%s<<",l->m->type);
1646          return -2;
1647        }
1648      }
1649    }
1650    // only ssi links:
1651    loop
1652    {
1653      /* yes: read 1 char*/
1654      /* if \n, check again with select else ungetc(c), ready*/
1655      /* setting: d: current ssiInfo, j current fd, i current entry in L*/
1656      int c=s_getc(d->f_read);
1657      //Print("try c=%d\n",c);
1658      if (c== -1) /* eof */
1659      {
1660        FD_CLR(j,&fdmask);
1661        fdmaskempty = 1;
1662        for(k = 0; k < max_fd; k++)
1663        {
1664          if(FD_ISSET(k, &fdmask))
1665          {
1666            fdmaskempty = 0;
1667            break;
1668          }
1669        }
1670        if(fdmaskempty)
1671        {
1672          return -1;
1673        }
1674        if(timeout != -1)
1675        {
1676          timeout = si_max(0,
1677             timeout - 1000000*(getRTimer()/TIMER_RESOLUTION - startingtime));
1678          wt.tv_sec  = timeout / 1000000;
1679          wt.tv_usec = (timeout % 1000000);
1680        }
1681        goto do_select;
1682      }
1683
1684      else if (isdigit(c))
1685      { s_ungetc(c,d->f_read); return i+1; }
1686      else if (c>' ')
1687      {
1688        Werror("unknown char in ssiLink(%d)",c);
1689        return -2;
1690      }
1691      /* else: next char */
1692    }
1693  }
1694}
1695
1696int ssiBatch(const char *host, const char * port)
1697/* return 0 on success, >0 else*/
1698{
1699  si_link l=(si_link) omAlloc0Bin(sip_link_bin);
1700  char *buf=(char*)omAlloc(256);
1701  sprintf(buf,"ssi:connect %s:%s",host,port);
1702  slInit(l, buf);
1703  if (slOpen(l,SI_LINK_OPEN,NULL)) return 1;
1704  SI_LINK_SET_RW_OPEN_P(l);
1705
1706  idhdl id = enterid(omStrDup("link_ll"), 0, LINK_CMD, &IDROOT, FALSE);
1707  IDLINK(id) = l;
1708
1709  loop
1710  {
1711    leftv h=ssiRead1(l); /*contains an exit.... */
1712    if (feErrors != NULL && *feErrors != '\0')
1713    {
1714      // handle errors:
1715      PrintS(feErrors); /* currently quite simple */
1716      *feErrors = '\0';
1717    }
1718    ssiWrite(l,h);
1719    h->CleanUp();
1720    omFreeBin(h, sleftv_bin);
1721  }
1722  /* never reached*/
1723  exit(0);
1724}
1725
1726static int ssiReserved_P=0;
1727static int ssiReserved_sockfd;
1728static  struct sockaddr_in ssiResverd_serv_addr;
1729static int  ssiReserved_Clients;
1730int ssiReservePort(int clients)
1731{
1732  if (ssiReserved_P!=0)
1733  {
1734    WerrorS("ERROR already a reverved port requested");
1735    return 0;
1736  }
1737  int portno;
1738  int n;
1739  ssiReserved_sockfd = socket(AF_INET, SOCK_STREAM, 0);
1740  if(ssiReserved_sockfd < 0)
1741  {
1742    WerrorS("ERROR opening socket");
1743    return 0;
1744  }
1745  memset((char *) &ssiResverd_serv_addr,0, sizeof(ssiResverd_serv_addr));
1746  portno = 1025;
1747  ssiResverd_serv_addr.sin_family = AF_INET;
1748  ssiResverd_serv_addr.sin_addr.s_addr = INADDR_ANY;
1749  do
1750  {
1751    portno++;
1752    ssiResverd_serv_addr.sin_port = htons(portno);
1753    if(portno > 50000)
1754    {
1755      WerrorS("ERROR on binding (no free port available?)");
1756      return 0;
1757    }
1758  }
1759  while(bind(ssiReserved_sockfd, (struct sockaddr *) &ssiResverd_serv_addr, sizeof(ssiResverd_serv_addr)) < 0);
1760  ssiReserved_P=portno;
1761  listen(ssiReserved_sockfd,clients);
1762  ssiReserved_Clients=clients;
1763  return portno;
1764}
1765
1766extern si_link_extension si_link_root;
1767si_link ssiCommandLink()
1768{
1769  if (ssiReserved_P==0)
1770  {
1771    WerrorS("ERROR no reverved port requested");
1772    return NULL;
1773  }
1774  struct sockaddr_in cli_addr;
1775  int clilen = sizeof(cli_addr);
1776  int newsockfd = si_accept(ssiReserved_sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
1777  if(newsockfd < 0)
1778  {
1779    Werror("ERROR on accept (errno=%d)",errno);
1780    return NULL;
1781  }
1782  si_link l=(si_link) omAlloc0Bin(sip_link_bin);
1783  si_link_extension s = si_link_root;
1784  si_link_extension prev = s;
1785  while (strcmp(s->type, "ssi") != 0)
1786  {
1787    if (s->next == NULL)
1788    {
1789      prev = s;
1790      s = NULL;
1791      break;
1792    }
1793    else
1794    {
1795      s = s->next;
1796    }
1797  }
1798  if (s != NULL)
1799    l->m = s;
1800  else
1801  {
1802    si_link_extension ns = (si_link_extension)omAlloc0Bin(s_si_link_extension_bin);
1803    prev->next=slInitSsiExtension(ns);
1804    l->m = prev->next;
1805  }
1806  l->name=omStrDup("");
1807  l->mode=omStrDup("tcp");
1808  l->ref=1;
1809  ssiInfo *d=(ssiInfo*)omAlloc0(sizeof(ssiInfo));
1810  l->data=d;
1811  d->fd_read = newsockfd;
1812  d->fd_write = newsockfd;
1813  d->f_read = s_open(newsockfd);
1814  d->f_write = fdopen(newsockfd, "w");
1815  SI_LINK_SET_RW_OPEN_P(l);
1816  ssiReserved_Clients--;
1817  if (ssiReserved_Clients<=0)
1818  {
1819    ssiReserved_P=0;
1820    si_close(ssiReserved_sockfd);
1821  }
1822  return l;
1823}
1824/*---------------------------------------------------------------------*/
1825/**
1826 * @brief additional default signal handler
1827
1828  // some newer Linux version cannot have SIG_IGN for SIGCHLD,
1829  // so use this nice routine here:
1830  //  SuSe 9.x reports -1 always
1831  //  Redhat 9.x/FC x reports sometimes -1
1832  // see also: hpux_system
1833  // also needed by getrusage (timer etc.)
1834
1835 @param[in] sig
1836**/
1837/*---------------------------------------------------------------------*/
1838void sig_chld_hdl(int sig)
1839{
1840  pid_t kidpid;
1841  int status;
1842
1843  loop
1844  {
1845    kidpid = si_waitpid(-1, &status, WNOHANG);
1846    if (kidpid==-1)
1847    {
1848      /* continue on interruption (EINTR): */
1849      if (errno == EINTR) continue;
1850      /* break on anything else (EINVAL or ECHILD according to manpage): */
1851      break;
1852    }
1853    else if (kidpid==0) break; /* no more children to process, so break */
1854
1855    //printf("Child %ld terminated\n", kidpid);
1856    link_list hh=ssiToBeClosed;
1857    while((hh!=NULL)&&(ssiToBeClosed_inactive))
1858    {
1859      if((hh->l!=NULL) && (hh->l->m->Open==ssiOpen))
1860      {
1861        ssiInfo *d = (ssiInfo *)hh->l->data;
1862        if(d->pid==kidpid)
1863        {
1864          if(ssiToBeClosed_inactive)
1865          {
1866            ssiToBeClosed_inactive=FALSE;
1867            slClose(hh->l);
1868            ssiToBeClosed_inactive=TRUE;
1869            break;
1870          }
1871          else break;
1872        }
1873        else hh=(link_list)hh->next;
1874      }
1875      else hh=(link_list)hh->next;
1876    }
1877  }
1878}
1879
1880static BOOLEAN DumpSsiIdhdl(si_link l, idhdl h)
1881{
1882  int type_id = IDTYP(h);
1883
1884  // C-proc not to be dumped, also LIB-proc not
1885  if (type_id == PROC_CMD)
1886  {
1887    if (IDPROC(h)->language == LANG_C) return FALSE;
1888    if (IDPROC(h)->libname != NULL) return FALSE;
1889  }
1890  // do not dump links
1891  if (type_id == LINK_CMD) return FALSE;
1892
1893  // do not dump ssi internal rings: ssiRing*
1894  if ((type_id == RING_CMD) && (strncmp(IDID(h),"ssiRing",7)==0))
1895    return FALSE;
1896
1897  command D=(command)omAlloc0(sizeof(*D));
1898  sleftv tmp;
1899  memset(&tmp,0,sizeof(tmp));
1900  tmp.rtyp=COMMAND;
1901  tmp.data=D;
1902
1903  if (type_id == PACKAGE_CMD)
1904  {
1905    // do not dump Top
1906    if (strcmp(IDID(h), "Top") == 0) return FALSE;
1907    package p=(package)IDDATA(h);
1908    // dump Singular-packages as load("...");
1909    if (p->language==LANG_SINGULAR)
1910    {
1911      D->op=LOAD_CMD;
1912      D->argc=1;
1913      D->arg1.rtyp=STRING_CMD;
1914      D->arg1.data=p->libname;
1915      ssiWrite(l,&tmp);
1916      omFree(D);
1917      return FALSE;
1918    }
1919  }
1920
1921  // handle qrings separately
1922  //if (type_id == QRING_CMD)
1923  //  return DumpSsiQringQring(l, h);
1924
1925  // put type and name
1926  //Print("generic dump:%s,%s\n",IDID(h),Tok2Cmdname(IDTYP(h)));
1927  D->op='=';
1928  D->argc=2;
1929  D->arg1.rtyp=DEF_CMD;
1930  D->arg1.name=IDID(h);
1931  D->arg2.rtyp=IDTYP(h);
1932  D->arg2.data=IDDATA(h);
1933  ssiWrite(l,&tmp);
1934  omFree(D);
1935  return FALSE;
1936}
1937static BOOLEAN ssiDumpIter(si_link l, idhdl h)
1938{
1939  if (h == NULL) return FALSE;
1940
1941  if (ssiDumpIter(l, IDNEXT(h))) return TRUE;
1942
1943  // need to set the ring before writing it, otherwise we get in
1944  // trouble with minpoly
1945  if (IDTYP(h) == RING_CMD || IDTYP(h) == QRING_CMD)
1946    rSetHdl(h);
1947
1948  if (DumpSsiIdhdl(l, h)) return TRUE;
1949
1950  // do not dump ssi internal rings: ssiRing*
1951  // but dump objects of all other rings
1952  if ((IDTYP(h) == RING_CMD || IDTYP(h) == QRING_CMD)
1953  && (strncmp(IDID(h),"ssiRing",7)!=0))
1954    return ssiDumpIter(l, IDRING(h)->idroot);
1955  else
1956    return FALSE;
1957}
1958BOOLEAN ssiDump(si_link l)
1959{
1960  idhdl h = IDROOT, rh = currRingHdl;
1961  BOOLEAN status = ssiDumpIter(l, h);
1962
1963  //if (! status ) status = DumpAsciiMaps(fd, h, NULL);
1964
1965  if (currRingHdl != rh) rSetHdl(rh);
1966  //fprintf(fd, "option(set, intvec(%d, %d));\n", si_opt_1, si_opt_2);
1967
1968  return status;
1969}
1970BOOLEAN ssiGetDump(si_link l)
1971{
1972  ssiInfo *d=(ssiInfo*)l->data;
1973  loop
1974  {
1975    if (!SI_LINK_OPEN_P(l)) break;
1976    if (s_iseof(d->f_read)) break;
1977    leftv h=ssiRead1(l); /*contains an exit.... */
1978    if (feErrors != NULL && *feErrors != '\0')
1979    {
1980      // handle errors:
1981      PrintS(feErrors); /* currently quite simple */
1982      return TRUE;
1983      *feErrors = '\0';
1984    }
1985    h->CleanUp();
1986    omFreeBin(h, sleftv_bin);
1987  }
1988  return FALSE;
1989}
1990// ----------------------------------------------------------------
1991// format
1992// 1 int %d
1993// 2 string <len> %s
1994// 3 number
1995// 4 bigint 4 %d or 3 <mpz_t>
1996// 5 ring
1997// 6 poly
1998// 7 ideal
1999// 8 matrix
2000// 9 vector
2001// 10 module
2002// 11 command
2003// 12 def <len> %s
2004// 13 proc <len> %s
2005// 14 list %d <elem1> ....
2006// 15 setring .......
2007// 16 nothing
2008// 17 intvec <len> ...
2009// 18 intmat
2010// 19 bigintmat <r> <c> ...
2011//
2012// 20 blackbox <name> ...
2013//
2014// 98: verify version: <ssi-version> <MAX_TOK> <OPT1> <OPT2>
2015// 99: quit Singular
Note: See TracBrowser for help on using the repository browser.