source: git/Singular/ssiLink.cc @ d828d63

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