source: git/Singular/ssiLink.cc @ 6d672b4

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