source: git/Singular/ssiLink.cc @ 854405

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