source: git/Singular/ipassign.cc @ b8fdd5

spielwiese
Last change on this file since b8fdd5 was b8fdd5, checked in by Hans Schönemann <hannes@…>, 26 years ago
* hannes: allowed assignment of ring to list elements git-svn-id: file:///usr/local/Singular/svn/trunk@1839 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 32.4 KB
Line 
1/****************************************
2*  Computer Algebra System SINGULAR     *
3****************************************/
4/* $Id: ipassign.cc,v 1.29 1998-05-18 17:23:00 Singular Exp $ */
5
6/*
7* ABSTRACT: interpreter:
8*           assignment of expressions and lists to objects or lists
9*/
10
11#include <stdlib.h>
12#include <string.h>
13#include <ctype.h>
14
15#include "mod2.h"
16#include "tok.h"
17#include "ipid.h"
18#include "intvec.h"
19#include "mmemory.h"
20#include "febase.h"
21#include "polys.h"
22#include "ideals.h"
23#include "matpol.h"
24#include "kstd1.h"
25#include "timer.h"
26#include "ring.h"
27#include "subexpr.h"
28#include "numbers.h"
29#include "lists.h"
30#include "longalg.h"
31#include "stairc.h"
32#include "maps.h"
33#include "syz.h"
34#include "weight.h"
35#include "ipconv.h"
36#include "attrib.h"
37#include "silink.h"
38#include "ipshell.h"
39
40/*=================== proc =================*/
41#ifdef SRING
42static BOOLEAN jjALTVARS(leftv res, leftv a)
43{
44  pAltVars=(int)a->Data();
45  pSRING=(pAltVars<=pVariables);
46  pSDRING=pSDRING||pSRING;
47  currRing->partN=pAltVars;
48  return FALSE;
49}
50#endif
51static BOOLEAN jjECHO(leftv res, leftv a)
52{
53  si_echo=(int)a->Data();
54  return FALSE;
55}
56static BOOLEAN jjPAGELENGTH(leftv res, leftv a)
57{
58  pagelength=(int)a->Data();
59  return FALSE;
60}
61static BOOLEAN jjPRINTLEVEL(leftv res, leftv a)
62{
63  printlevel=(int)a->Data();
64  return FALSE;
65}
66static BOOLEAN jjCOLMAX(leftv res, leftv a)
67{
68  colmax=(int)a->Data();
69  return FALSE;
70}
71static BOOLEAN jjTIMER(leftv res, leftv a)
72{
73  timerv=(int)a->Data();
74  initTimer();
75  return FALSE;
76}
77#ifdef HAVE_RTIMER
78static BOOLEAN jjRTIMER(leftv res, leftv a)
79{
80  rtimerv=(int)a->Data();
81  initRTimer();
82  return FALSE;
83}
84#endif
85static BOOLEAN jjMAXDEG(leftv res, leftv a)
86{
87  Kstd1_deg=(int)a->Data();
88  if (Kstd1_deg!=0)
89    test |=Sy_bit(24);
90  else
91    test &=(~Sy_bit(24));
92  return FALSE;
93}
94static BOOLEAN jjMAXMULT(leftv res, leftv a)
95{
96  Kstd1_mu=(int)a->Data();
97  if (Kstd1_mu!=0)
98    test |=Sy_bit(23);
99  else
100    test &=(~Sy_bit(23));
101  return FALSE;
102}
103static BOOLEAN jjTRACE(leftv res, leftv a)
104{
105  traceit=(int)a->Data();
106  return FALSE;
107}
108static BOOLEAN jjSHORTOUT(leftv res, leftv a)
109{
110#ifdef HAVE_TCL
111  if (!tclmode)
112#endif
113    pShortOut=(int)a->Data();
114  return FALSE;
115}
116static BOOLEAN jjMINPOLY(leftv res, leftv a)
117{
118  if ((currRing->parameter!=NULL)
119  && ((currRing->P>1)
120    || (currRing->ch>1/* GF(q) <=> parameter!=NULL && ch=q>1*/)))
121  {
122    WerrorS("no minpoly allowed");
123    return TRUE;
124  }
125  if (currRing->minpoly!=NULL)
126  {
127    WerrorS("minpoly already set");
128    return TRUE;
129  }
130  number p=(number)a->CopyD(NUMBER_CMD);
131  if (!nIsZero(p))
132  {
133    nNormalize(p);
134    currRing->minpoly=p;
135    naMinimalPoly=((lnumber)currRing->minpoly)->z;
136  }
137  else
138  {
139    naMinimalPoly=NULL;
140  }
141  return FALSE;
142}
143static BOOLEAN jjNOETHER(leftv res, leftv a)
144{
145  poly p=(poly)a->CopyD(POLY_CMD);
146  pDelete(&ppNoether);
147  ppNoether=p;
148  return FALSE;
149}
150/*=================== proc =================*/
151static void jiAssignAttr(leftv l,leftv r)
152{
153  // get the attribute of th right side
154  // and set it to l
155  leftv rv=r->LData();
156  if (rv!=NULL)
157  {
158    if (rv->e==NULL)
159    {
160      if (rv->attribute!=NULL)
161      {
162        attr la;
163        if (r->rtyp!=IDHDL)
164        {
165          la=rv->attribute;
166          rv->attribute=NULL;
167        }
168        else
169        {
170          la=rv->attribute->Copy();
171        }
172        l->attribute=la;
173      }
174      l->flag=rv->flag;
175    }
176  }
177  if (l->rtyp==IDHDL)
178  {
179    idhdl h=(idhdl)l->data;
180    IDATTR(h)=l->attribute;
181    IDFLAG(h)=l->flag;
182  }
183}
184static BOOLEAN jiA_INT(leftv res, leftv a, Subexpr e)
185{
186  if (e==NULL)
187  {
188    res->data=(void *)a->Data();
189    jiAssignAttr(res,a);
190  }
191  else
192  {
193    int i=e->start-1;
194    if (i<0)
195    {
196      Werror("index[%d] must be positive",i+1);
197      return TRUE;
198    }
199    intvec *iv=(intvec *)res->data;
200    if (e->next==NULL)
201    {
202      if (i>=iv->length())
203      {
204        intvec *iv1=new intvec(i+1);
205        (*iv1)[i]=(int)a->Data();
206        intvec *ivn=ivAdd(iv,iv1);
207        delete iv;
208        delete iv1;
209        res->data=(void *)ivn;
210      }
211      else
212        (*iv)[i]=(int)a->Data();
213    }
214    else
215    {
216      int c=e->next->start;
217      if ((i>=iv->rows())||(c<1)||(c>iv->cols()))
218      {
219        Werror("wrong range [%d,%d] in intmat (%d,%d)",i+1,c,iv->rows(),iv->cols());
220        return TRUE;
221      }
222      else
223        IMATELEM(*iv,i+1,c) = (int)a->Data();
224    }
225  }
226  return FALSE;
227}
228static BOOLEAN jiA_NUMBER(leftv res, leftv a, Subexpr e)
229{
230  number p=(number)a->CopyD(NUMBER_CMD);
231  if (res->data!=NULL) nDelete((number *)&res->data);
232  nNormalize(p);
233  res->data=(void *)p;
234  jiAssignAttr(res,a);
235  return FALSE;
236}
237static BOOLEAN jiA_LIST_RES(leftv res, leftv a,Subexpr e)
238{
239  syStrategy r=(syStrategy)a->CopyD();
240  if (res->data!=NULL) ((lists)res->data)->Clean();
241  res->data=(void *)syConvRes(r);
242  //jiAssignAttr(res,a);
243  return FALSE;
244}
245static BOOLEAN jiA_LIST(leftv res, leftv a,Subexpr e)
246{
247  lists l=(lists)a->CopyD();
248  if (res->data!=NULL) ((lists)res->data)->Clean();
249  res->data=(void *)l;
250  jiAssignAttr(res,a);
251  return FALSE;
252}
253static BOOLEAN jiA_POLY(leftv res, leftv a,Subexpr e)
254{
255  poly p=(poly)a->CopyD(POLY_CMD);
256  pNormalize(p);
257  if (e==NULL)
258  {
259    if (res->data!=NULL) pDelete((poly*)&res->data);
260    res->data=(void*)p;
261    jiAssignAttr(res,a);
262  }
263  else
264  {
265    int i,j;
266    matrix m=(matrix)res->data;
267    i=e->start;
268    if (e->next==NULL)
269    {
270      j=i; i=1;
271      // for ideal/module: check indices
272      if (res->rtyp!=MAP_CMD)
273      {
274        if (j>MATCOLS(m))
275        {
276          pEnlargeSet(&(m->m),MATCOLS(m),j-MATCOLS(m));
277          MATCOLS(m)=j;
278        }
279        else if (j<=0)
280        {
281          Werror("index[%d] must be positive",j/*e->start*/);
282          return TRUE;
283        }
284      }
285    }
286    else
287    {
288      // for matrices: indices are correct (see ipExprArith3(..,'['..) )
289      j=e->next->start;
290    }
291    pDelete(&MATELEM(m,i,j));
292    MATELEM(m,i,j)=p;
293    /* for module: update rank */
294    if ((p!=NULL) && (pGetComp(p)!=0))
295    {
296      m->rank=max(m->rank,pMaxComp(p));
297    }
298  }
299  return FALSE;
300}
301static BOOLEAN jiA_1x1MATRIX(leftv res, leftv a,Subexpr e)
302{
303  if ((res->rtyp!=MATRIX_CMD) /*|| (e!=NULL) - TRUE because of type poly */)
304     return TRUE;
305  matrix am=(matrix)a->CopyD(MATRIX_CMD);
306  if ((MATROWS(am)!=1) || (MATCOLS(am)!=1))
307  {
308    idDelete((ideal *)&am);
309    return TRUE;
310  }
311  matrix m=(matrix)res->data;
312  // indices are correct (see ipExprArith3(..,'['..) )
313  int i=e->start;
314  int j=e->next->start;
315  pDelete(&MATELEM(m,i,j));
316  pNormalize(MATELEM(am,1,1));
317  MATELEM(m,i,j)=MATELEM(am,1,1);
318  return FALSE;
319}
320static BOOLEAN jiA_STRING(leftv res, leftv a, Subexpr e)
321{
322  if (e==NULL)
323  {
324    FreeL((ADDRESS)res->data);
325    res->data=(void *)a->CopyD(STRING_CMD);
326    jiAssignAttr(res,a);
327  }
328  else
329  {
330    char *s=(char *)res->data;
331    if ((e->start>0)&&(e->start<=(int)strlen(s)))
332      s[e->start-1]=(char)(*((char *)a->Data()));
333    else
334    {
335      Werror("string index %d out of range 1..%d",e->start,strlen(s));
336      return TRUE;
337    }
338  }
339  return FALSE;
340}
341static BOOLEAN jiA_PROC(leftv res, leftv a, Subexpr e)
342{
343  extern procinfo *iiInitSingularProcinfo(procinfo *pi, char *libname,
344                                          char *procname, int line,
345                                          long pos, BOOLEAN pstatic=FALSE);
346  extern void piCleanUp(procinfov pi);
347
348  if(res->data!=NULL) piCleanUp((procinfo *)res->data);
349  if(a->rtyp==STRING_CMD)
350  {
351    res->data = (void *)Alloc0(sizeof(procinfo));
352    ((procinfo *)(res->data))->language=LANG_NONE;
353    iiInitSingularProcinfo((procinfo *)res->data,"",res->name,0,0);
354    ((procinfo *)res->data)->data.s.body=(char *)a->CopyD(STRING_CMD);
355  }
356  else
357    res->data=(void *)a->CopyD(PROC_CMD);
358  jiAssignAttr(res,a);
359  return FALSE;
360}
361static BOOLEAN jiA_INTVEC(leftv res, leftv a, Subexpr e)
362{
363  if (res->data!=NULL) delete ((intvec *)res->data);
364  res->data=(void *)a->CopyD(INTVEC_CMD);
365  jiAssignAttr(res,a);
366  return FALSE;
367}
368static BOOLEAN jiA_IDEAL(leftv res, leftv a, Subexpr e)
369{
370  if (res->data!=NULL) idDelete((ideal*)&res->data);
371  res->data=(void *)a->CopyD(MATRIX_CMD);
372  jiAssignAttr(res,a);
373  return FALSE;
374}
375static BOOLEAN jiA_RESOLUTION(leftv res, leftv a, Subexpr e)
376{
377  if (res->data!=NULL) syKillComputation((syStrategy)res->data);
378  res->data=(void *)a->CopyD(RESOLUTION_CMD);
379  jiAssignAttr(res,a);
380  return FALSE;
381}
382static BOOLEAN jiA_MODUL_P(leftv res, leftv a, Subexpr e)
383{
384  if (res->data!=NULL) idDelete((ideal*)&res->data);
385  ideal I=idInit(1,1);
386  I->m[0]=(poly)a->CopyD();
387  if (I->m[0]!=NULL) pSetComp(I->m[0],1);
388  res->data=(void *)I;
389  return FALSE;
390}
391static BOOLEAN jiA_IDEAL_M(leftv res, leftv a, Subexpr e)
392{
393  if (res->data!=NULL) idDelete((ideal*)&res->data);
394  matrix m=(matrix)a->CopyD(MATRIX_CMD);
395  IDELEMS((ideal)m)=MATROWS(m)*MATCOLS(m);
396  ((ideal)m)->rank=1;
397  MATROWS(m)=1;
398  res->data=(void *)m;
399  return FALSE;
400}
401static BOOLEAN jiA_LINK(leftv res, leftv a, Subexpr e)
402{
403  si_link l=(si_link)res->data;
404
405  if (l!=NULL) slCleanUp(l);
406
407  if (a->Typ() == STRING_CMD)
408  {
409    if (l == NULL)
410    {
411      l = (si_link) Alloc0(sizeof(sip_link));
412      res->data = (void *) l;
413    }
414    return slInit(l, (char *) a->Data());
415  }
416  else if (a->Typ() == LINK_CMD)
417  {
418    if (l != NULL) Free(l, sizeof(sip_link));
419    res->data = slCopy((si_link)a->Data());
420    return FALSE;
421  }
422  return TRUE;
423}
424static BOOLEAN jiA_MAP(leftv res, leftv a, Subexpr e)
425{
426  if (res->data!=NULL)
427  {
428    FreeL((ADDRESS)((map)res->data)->preimage);
429    ((map)res->data)->preimage=NULL;
430    idDelete((ideal*)&res->data);
431  }
432  res->data=(void *)a->CopyD();
433  jiAssignAttr(res,a);
434  return FALSE;
435}
436static BOOLEAN jiA_MAP_ID(leftv res, leftv a, Subexpr e)
437{
438  map f=(map)res->data;
439  char *rn=f->preimage; // save the old/already assigned preimage ring name
440  f->preimage=NULL;
441  idDelete((ideal *)&f);
442  res->data=(void *)a->CopyD(IDEAL_CMD);
443  f=(map)res->data;
444  f->preimage = rn;
445  return FALSE;
446}
447static BOOLEAN jiA_QRING(leftv res, leftv a,Subexpr e)
448{
449  // the follwing can only happen, if:
450  //   - the left side is of type qring AND not an id
451  if ((e!=NULL)||(res->rtyp!=IDHDL))
452  {
453    WerrorS("qring_id expected");
454    return TRUE;
455  }
456  ring qr;
457  int i,j;
458  int *pi;
459
460  assumeStdFlag(a);
461  qr=(ring)res->Data();
462  ring qrr=rCopy(currRing);
463  // hier fehlt noch: evtl. vorhandenen alten Ring streichen
464  // vorerst nur:
465  memcpy4(qr,qrr,sizeof(ip_sring));
466  Free(qrr,sizeof(ip_sring));
467  if (qr->qideal!=NULL) idDelete(&qr->qideal);
468  qr->qideal = (ideal)a->CopyD(IDEAL_CMD);
469  currRing=qr;
470  currRingHdl=(idhdl)res->data;
471  currQuotient=qr->qideal;
472  return FALSE;
473}
474static BOOLEAN jiA_RING(leftv res, leftv a, Subexpr e)
475{
476  BOOLEAN have_id=TRUE;
477  if ((e!=NULL)||(res->rtyp!=IDHDL))
478  {
479    //WerrorS("id expected");
480    //return TRUE;
481    have_id=FALSE;
482  }
483  ring r=(ring)a->Data();
484  if (have_id)
485  {
486    idhdl rl=(idhdl)res->data;
487    if (&IDRING(rl)!=NULL) rKill(rl);
488    IDRING(rl)=r;
489    if ((IDLEV((idhdl)a->data)!=myynest) && (r==currRing))
490      currRingHdl=(idhdl)res->data;
491  }
492  else
493  {
494    if (e==NULL) res->data=(char *)r;
495    else
496    {
497      WerrorS("id expected");
498      return TRUE;
499    }
500  }
501  r->ref++;
502  return FALSE;
503}
504/*=================== table =================*/
505struct sValAssign dAssign[]=
506{
507// proc         res             arg
508 {jiA_IDEAL,    IDEAL_CMD,      IDEAL_CMD }
509,{jiA_IDEAL_M,  IDEAL_CMD,      MATRIX_CMD }
510,{jiA_RESOLUTION,RESOLUTION_CMD,RESOLUTION_CMD }
511,{jiA_INT,      INT_CMD,        INT_CMD }
512,{jiA_IDEAL,    MATRIX_CMD,     MATRIX_CMD }
513,{jiA_MAP_ID,   MAP_CMD,        IDEAL_CMD }
514,{jiA_MAP,      MAP_CMD,        MAP_CMD }
515,{jiA_IDEAL,    MODUL_CMD,      MODUL_CMD }
516,{jiA_MODUL_P,  MODUL_CMD,      POLY_CMD }
517,{jiA_POLY,     POLY_CMD,       POLY_CMD }
518,{jiA_1x1MATRIX,POLY_CMD,       MATRIX_CMD }
519,{jiA_QRING,    QRING_CMD,      IDEAL_CMD }
520,{jiA_RING,     RING_CMD,       RING_CMD }
521,{jiA_RING,     QRING_CMD,      QRING_CMD }
522,{jiA_STRING,   STRING_CMD,     STRING_CMD }
523,{jiA_PROC,     PROC_CMD,       STRING_CMD }
524,{jiA_PROC,     PROC_CMD,       PROC_CMD }
525,{jiA_POLY,     VECTOR_CMD,     VECTOR_CMD }
526,{jiA_INTVEC,   INTVEC_CMD,     INTVEC_CMD }
527,{jiA_INTVEC,   INTMAT_CMD,     INTMAT_CMD }
528,{jiA_NUMBER,   NUMBER_CMD,     NUMBER_CMD }
529,{jiA_LIST_RES, LIST_CMD,       RESOLUTION_CMD }
530,{jiA_LIST,     LIST_CMD,       LIST_CMD }
531,{jiA_LINK,     LINK_CMD,       STRING_CMD }
532,{jiA_LINK,     LINK_CMD,       LINK_CMD }
533,{NULL,         0,              0 }
534};
535struct sValAssign_sys dAssign_sys[]=
536{
537// sysvars:
538 {jjECHO,       VECHO,          INT_CMD }
539#ifdef SRING
540,{jjALTVARS,    VALTVARS,       INT_CMD }
541#endif
542,{jjPAGELENGTH, VPAGELENGTH,    INT_CMD }
543,{jjPRINTLEVEL, VPRINTLEVEL,    INT_CMD }
544,{jjCOLMAX,     VCOLMAX,        INT_CMD }
545,{jjTIMER,      VTIMER,         INT_CMD }
546#ifdef HAVE_RTIMER
547,{jjRTIMER,     VRTIMER,        INT_CMD }
548#endif
549,{jjMAXDEG,     VMAXDEG,        INT_CMD }
550,{jjMAXMULT,    VMAXMULT,       INT_CMD }
551,{jjTRACE,      TRACE,          INT_CMD }
552,{jjSHORTOUT,   VSHORTOUT,      INT_CMD }
553,{jjMINPOLY,    VMINPOLY,       NUMBER_CMD }
554,{jjNOETHER,    VNOETHER,       POLY_CMD }
555,{NULL,         0,              0 }
556};
557/*=================== operations ============================*/
558/*2
559* assign a = b
560*/
561static BOOLEAN jiAssign_1(leftv l, leftv r)
562{
563  int rt=r->Typ();
564  if (rt==0)
565  {
566    if (!errorreported) Werror("`%s` is undefined",r->Name());
567    return TRUE;
568  }
569
570  int lt=l->Typ();
571  if((lt==0)/*&&(l->name!=NULL)*/)
572  {
573    if (!errorreported) Werror("left side `%s` is undefined",l->Name());
574    return TRUE;
575  }
576  if((rt==DEF_CMD)||(rt==NONE))
577  {
578    if (!errorreported) WerrorS("right side is not a datum");
579    return TRUE;
580  }
581
582  int i=0;
583  BOOLEAN nok=FALSE;
584
585  if (lt==DEF_CMD)
586  {
587    if (l->rtyp==IDHDL)
588    {
589      IDTYP((idhdl)l->data)=rt;
590    }
591    else if (l->name!=NULL)
592    {
593      sleftv ll;
594      iiDeclCommand(&ll,l,myynest,rt,&idroot);
595      memcpy(l,&ll,sizeof(sleftv));
596    }
597    else
598    {
599      l->rtyp=rt;
600    }
601    lt=rt;
602  }
603  else
604  {
605    if ((l->data==r->data)&&(l->e==NULL)&&(r->e==NULL))
606      return FALSE;
607  }
608  leftv ld=l;
609  if ((l->rtyp==IDHDL)&&(lt!=QRING_CMD)&&(lt!=RING_CMD))
610    ld=(leftv)l->data;
611  while (((dAssign[i].res!=lt)
612      || (dAssign[i].arg!=rt))
613    && (dAssign[i].res!=0)) i++;
614  if (dAssign[i].res!=0)
615  {
616    BOOLEAN b;
617    b=dAssign[i].p(ld,r,l->e);
618    if(l!=ld) /* i.e. l is IDHDL, l->data is ld */
619    {
620      l->flag=ld->flag;
621      l->attribute=ld->attribute;
622    }
623    return b;
624  }
625  // implicite type conversion ----------------------------------------------
626  if (dAssign[i].res==0)
627  {
628    int ri;
629    leftv rn = (leftv)Alloc0(sizeof(sleftv));
630    BOOLEAN failed=FALSE;
631    i=0;
632    while ((dAssign[i].res!=lt)
633      && (dAssign[i].res!=0)) i++;
634    while (dAssign[i].res==lt)
635    {
636      if ((ri=iiTestConvert(rt,dAssign[i].arg))!=0)
637      {
638        failed= iiConvert(rt,dAssign[i].arg,ri,r,rn);
639        if(!failed)
640        {
641          failed= dAssign[i].p(ld,rn,l->e);
642        }
643        // everything done, clean up temp. variables
644        rn->CleanUp();
645        Free((ADDRESS)rn,sizeof(sleftv));
646        if (failed)
647        {
648          // leave loop, goto error handling
649          break;
650        }
651        else
652        {
653          if(l!=ld) /* i.e. l is IDHDL, l->data is ld */
654          {
655            l->flag=ld->flag;
656            l->attribute=ld->attribute;
657          }
658          // everything ok, return
659          return FALSE;
660        }
661     }
662     i++;
663    }
664    // error handling ---------------------------------------------------
665    if (!errorreported)
666    {
667      Werror("`%s` = `%s` is not supported"
668             ,Tok2Cmdname(lt),Tok2Cmdname(rt));
669      if (BVERBOSE(V_SHOW_USE))
670      {
671        i=0;
672        while ((dAssign[i].res!=lt)
673          && (dAssign[i].res!=0)) i++;
674        while (dAssign[i].res==lt)
675        {
676          Werror("expected `%s` = `%s`"
677              ,Tok2Cmdname(lt),Tok2Cmdname(dAssign[i].arg));
678          i++;
679        }
680      }
681    }
682  }
683  return TRUE;
684}
685/*2
686* assign sys_var = val
687*/
688static BOOLEAN iiAssign_sys(leftv l, leftv r)
689{
690  int rt=r->Typ();
691
692  if (rt==0)
693  {
694    if (!errorreported) Werror("`%s` is undefined",r->Name());
695    return TRUE;
696  }
697  int i=0;
698  int lt=l->rtyp;
699  while (((dAssign_sys[i].res!=lt)
700      || (dAssign_sys[i].arg!=rt))
701    && (dAssign_sys[i].res!=0)) i++;
702  if (dAssign_sys[i].res!=0)
703  {
704    if (!dAssign_sys[i].p(l,r))
705    {
706      // everything ok, clean up
707      return FALSE;
708    }
709  }
710  // implicite type conversion ----------------------------------------------
711  if (dAssign_sys[i].res==0)
712  {
713    int ri;
714    leftv rn = (leftv)Alloc0(sizeof(sleftv));
715    BOOLEAN failed=FALSE;
716    i=0;
717    while ((dAssign_sys[i].res!=lt)
718      && (dAssign_sys[i].res!=0)) i++;
719    while (dAssign_sys[i].res==lt)
720    {
721      if ((ri=iiTestConvert(rt,dAssign_sys[i].arg))!=0)
722      {
723        failed= ((iiConvert(rt,dAssign_sys[i].arg,ri,r,rn))
724            || (dAssign_sys[i].p(l,rn)));
725        // everything done, clean up temp. variables
726        rn->CleanUp();
727        Free((ADDRESS)rn,sizeof(sleftv));
728        if (failed)
729        {
730          // leave loop, goto error handling
731          break;
732        }
733        else
734        {
735          // everything ok, return
736          return FALSE;
737        }
738     }
739     i++;
740    }
741    // error handling ---------------------------------------------------
742    if(!errorreported)
743    {
744      Werror("`%s` = `%s` is not supported"
745             ,Tok2Cmdname(lt),Tok2Cmdname(rt));
746      if (BVERBOSE(V_SHOW_USE))
747      {
748        i=0;
749        while ((dAssign_sys[i].res!=lt)
750          && (dAssign_sys[i].res!=0)) i++;
751        while (dAssign_sys[i].res==lt)
752        {
753          Werror("expected `%s` = `%s`"
754              ,Tok2Cmdname(lt),Tok2Cmdname(dAssign_sys[i].arg));
755          i++;
756        }
757      }
758    }
759  }
760  return TRUE;
761}
762static BOOLEAN jiA_INTVEC_L(leftv l,leftv r)
763{
764  /* right side is intvec, left side is list (of int)*/
765  BOOLEAN nok;
766  int i=0;
767  leftv l1=l;
768  leftv h;
769  sleftv t;
770  intvec *iv=(intvec *)r->Data();
771  memset(&t,0,sizeof(sleftv));
772  t.rtyp=INT_CMD;
773  while ((i<iv->length())&&(l!=NULL))
774  {
775    t.data=(char *)(*iv)[i];
776    h=l->next;
777    l->next=NULL;
778    nok=jiAssign_1(l,&t);
779    if (nok) return TRUE;
780    i++;
781    l=h;
782  }
783  l1->CleanUp();
784  r->CleanUp();
785  return FALSE;
786}
787static BOOLEAN jiA_VECTOR_L(leftv l,leftv r)
788{
789  /* right side is vector, left side is list (of poly)*/
790  BOOLEAN nok;
791  leftv l1=l;
792  ideal I=idVec2Ideal((poly)r->Data());
793  leftv h;
794  sleftv t;
795  int i=0;
796  while (l!=NULL)
797  {
798    memset(&t,0,sizeof(sleftv));
799    t.rtyp=POLY_CMD;
800    if (i>=IDELEMS(I))
801    {
802      t.data=NULL;
803    }
804    else
805    {
806      t.data=(char *)I->m[i];
807      I->m[i]=NULL;
808    }
809    h=l->next;
810    l->next=NULL;
811    nok=jiAssign_1(l,&t);
812    t.CleanUp();
813    if (nok)
814    {
815      idDelete(&I);
816      return TRUE;
817    }
818    i++;
819    l=h;
820  }
821  idDelete(&I);
822  l1->CleanUp();
823  r->CleanUp();
824  return FALSE;
825}
826static BOOLEAN jjA_L_LIST(leftv l, leftv r)
827/* left side: list
828*  right side: expression list
829*/
830{
831  int sl = r->listLength();
832  lists L=(lists)Alloc(sizeof(slists));
833  leftv h=NULL,o_r=r;
834  int i;
835  int rt;
836
837  L->Init(sl);
838  for (i=0;i<sl;i++)
839  {
840    if (h!=NULL) { /* e.g. not in the first step:
841                   * h is the pointer to the old sleftv,
842                   * r is the pointer to the next sleftv
843                   * (in this moment) */
844                   h->next=r;
845                 }
846    h=r;
847    r=r->next;
848    h->next=NULL;
849    rt=h->Typ();
850    if ((rt==0)||(rt==NONE))
851    {
852      L->Clean();
853      Werror("`%s` is undefined",h->Name());
854      return TRUE;
855    }
856    //if ((rt==RING_CMD)||(rt==QRING_CMD))
857    //{
858    //  L->m[i].rtyp=rt;
859    //  L->m[i].data=h->Data();
860    //  ((ring)L->m[i].data)->ref++;
861    //}
862    //else
863      L->m[i].Copy(h);
864  }
865  IDLIST((idhdl)l->data)->Clean();
866  IDLIST((idhdl)l->data)=L;
867  ipMoveId((idhdl)l->data);
868  o_r->CleanUp();
869  return FALSE;
870}
871static BOOLEAN jiA_L_LIST(leftv l, leftv r)
872/* left side: list
873*  right side: expression list
874*/
875{
876  int sl = r->listLength();
877  lists L=(lists)Alloc(sizeof(slists));
878  leftv h=NULL;
879  int i;
880  int rt;
881
882  L->Init(sl);
883  for (i=0;i<sl;i++)
884  {
885    if (h!=NULL) { /* e.g. not in the first step:
886                   * h is the pointer to the old sleftv,
887                   * r is the pointer to the next sleftv
888                   * (in this moment) */
889                   h->next=r;
890                 }
891    h=r;
892    r=r->next;
893    h->next=NULL;
894    rt=h->Typ();
895    if (rt==0)
896    {
897      L->Clean();
898      Werror("`%s` is undefined",h->Name());
899      return TRUE;
900    }
901    if ((rt==RING_CMD)||(rt==QRING_CMD))
902    {
903      L->m[i].rtyp=rt;
904      L->m[i].data=h->Data();
905      ((ring)L->m[i].data)->ref++;
906    }
907    else
908      L->m[i].Copy(h);
909  }
910  IDLIST((idhdl)l->data)->Clean();
911  IDLIST((idhdl)l->data)=L;
912  return FALSE;
913}
914static BOOLEAN jjA_L_INTVEC(leftv l,leftv r,intvec *iv)
915{
916  /* left side is intvec/intmat, right side is list (of int,intvec,intmat)*/
917  leftv hh=r;
918  int i = 0;
919  while (hh!=NULL)
920  {
921    if (i>=iv->length()) break;
922    if (hh->Typ() == INT_CMD)
923    {
924      (*iv)[i++] = (int)(hh->Data());
925    }
926    else if ((hh->Typ() == INTVEC_CMD)
927            ||(hh->Typ() == INTMAT_CMD))
928    {
929      intvec *ivv = (intvec *)(hh->Data());
930      int ll = 0,l = min(ivv->length(),iv->length());
931      for (; l>0; l--)
932      {
933        (*iv)[i++] = (*ivv)[ll++];
934      }
935    }
936    else
937    {
938      delete iv;
939      return TRUE;
940    }
941    hh = hh->next;
942  }
943  if (IDINTVEC((idhdl)l->data)!=NULL) delete IDINTVEC((idhdl)l->data);
944  IDINTVEC((idhdl)l->data)=iv;
945  return FALSE;
946}
947static BOOLEAN jjA_L_STRING(leftv l,leftv r)
948{
949  /* left side is string, right side is list of string*/
950  leftv hh=r;
951  int sl = 1;
952  char *s;
953  char *t;
954  int tl;
955  /* find the length */
956  while (hh!=NULL)
957  {
958    if (hh->Typ()!= STRING_CMD)
959    {
960      return TRUE;
961    }
962    sl += strlen((char *)hh->Data());
963    hh = hh->next;
964  }
965  s = (char * )AllocL(sl);
966  sl=0;
967  hh = r;
968  while (hh!=NULL)
969  {
970    t=(char *)hh->Data();
971    tl=strlen(t);
972    memcpy(s+sl,t,tl);
973    sl+=tl;
974    hh = hh->next;
975  }
976  s[sl]='\0';
977  FreeL((ADDRESS)IDDATA((idhdl)(l->data)));
978  IDDATA((idhdl)(l->data))=s;
979  return FALSE;
980}
981static BOOLEAN jjA_LIST_L(leftv l,leftv r)
982{
983  /*left side are something, right side are lists*/
984  /*e.g. a,b,c=l */
985  //int ll=l->listLength();
986  if (l->listLength()==1) return jiAssign_1(l,r);
987  BOOLEAN nok;
988  sleftv t;
989  leftv h;
990  lists L=(lists)r->Data();
991  int rl=L->nr;
992  int i=0;
993
994  memset(&t,0,sizeof(sleftv));
995  while ((i<=rl)&&(l!=NULL))
996  {
997    memset(&t,0,sizeof(sleftv));
998    t.Copy(&L->m[i]);
999    h=l->next;
1000    l->next=NULL;
1001    nok=jiAssign_1(l,&t);
1002    if (nok) return TRUE;
1003    i++;
1004    l=h;
1005  }
1006  r->CleanUp();
1007  return FALSE;
1008}
1009static BOOLEAN jiA_MATRIX_L(leftv l,leftv r)
1010{
1011  /* right side is matrix, left side is list (of poly)*/
1012  BOOLEAN nok=FALSE;
1013  int i;
1014  matrix m=(matrix)r->CopyD(MATRIX_CMD);
1015  leftv h;
1016  leftv ol=l;
1017  leftv o_r=r;
1018  sleftv t;
1019  memset(&t,0,sizeof(sleftv));
1020  t.rtyp=POLY_CMD;
1021  loop
1022  {
1023    i=0;
1024    while ((i<MATROWS(m)*MATCOLS(m))&&(l!=NULL))
1025    {
1026      t.data=(char *)m->m[i];
1027      m->m[i]=NULL;
1028      h=l->next;
1029      l->next=NULL;
1030      nok=jiAssign_1(l,&t);
1031      l->next=h;
1032      if (nok)
1033      {
1034        idDelete((ideal *)&m);
1035        goto ende;
1036      }
1037      i++;
1038      l=h;
1039    }
1040    idDelete((ideal *)&m);
1041    h=r;
1042    r=r->next;
1043    if (l==NULL)
1044    {
1045      if (r!=NULL)
1046      {
1047        Warn("list length mismatch in assign (l>r)");
1048        nok=TRUE;
1049      }
1050      break;
1051    }
1052    else if (r==NULL)
1053    {
1054      Warn("list length mismatch in assign (l<r)");
1055      nok=TRUE;
1056      break;
1057    }
1058    if ((r->Typ()==IDEAL_CMD)||(r->Typ()==MATRIX_CMD))
1059      m=(matrix)r->CopyD(MATRIX_CMD);
1060    else if (r->Typ()==POLY_CMD)
1061    {
1062      m=mpNew(1,1);
1063      MATELEM(m,1,1)=(poly)r->CopyD(POLY_CMD);
1064    }
1065    else
1066    {
1067      nok=TRUE;
1068      break;
1069    }
1070  }
1071ende:
1072  o_r->CleanUp();
1073  ol->CleanUp();
1074  return nok;
1075}
1076static BOOLEAN jiA_STRING_L(leftv l,leftv r)
1077{
1078  /*left side are strings, right side is a string*/
1079  /*e.g. s[2..3]="12" */
1080  /*the case s=t[1..4] is handled in iiAssign,
1081  * the case s[2..3]=t[3..4] is handled in iiAssgn_rec*/
1082  int ll=l->listLength();
1083  int rl=r->listLength();
1084  BOOLEAN nok=FALSE;
1085  sleftv t;
1086  leftv h,l1=l;
1087  int i=0;
1088  char *ss;
1089  char *s=(char *)r->Data();
1090  int sl=strlen(s);
1091
1092  memset(&t,0,sizeof(sleftv));
1093  t.rtyp=STRING_CMD;
1094  while ((i<sl)&&(l!=NULL))
1095  {
1096    ss=(char *)AllocL(2);
1097    ss[1]='\0';
1098    ss[0]=s[i];
1099    t.data=ss;
1100    h=l->next;
1101    l->next=NULL;
1102    nok=jiAssign_1(l,&t);
1103    if (nok)
1104    {
1105      break;
1106    }
1107    i++;
1108    l=h;
1109  }
1110  r->CleanUp();
1111  l1->CleanUp();
1112  return nok;
1113}
1114static BOOLEAN jiAssign_list(leftv l, leftv r)
1115{
1116  int i=l->e->start-1;
1117  if (i<0)
1118  {
1119    Werror("index[%d] must be positive",i+1);
1120    return TRUE;
1121  }
1122  if(l->attribute!=NULL)
1123  {
1124    atKillAll((idhdl)l);
1125    l->attribute=NULL;
1126  }
1127  l->flag=0;
1128  lists li;
1129  if (l->rtyp==IDHDL)
1130  {
1131    li=IDLIST((idhdl)l->data);
1132  }
1133  else
1134  {
1135    li=(lists)l->data;
1136  }
1137  if (i>li->nr)
1138  {
1139    li->m=(leftv)ReAlloc(li->m,(li->nr+1)*sizeof(sleftv),(i+1)*sizeof(sleftv));
1140    memset(&(li->m[li->nr+1]),0,(i-li->nr)*sizeof(sleftv));
1141    int j=li->nr+1;
1142    for(;j<=i;j++)
1143      li->m[j].rtyp=DEF_CMD;
1144    li->nr=i;
1145  }
1146  leftv ld=&(li->m[i]);
1147  ld->e=l->e->next;
1148  BOOLEAN b;
1149  if ((ld->rtyp!=LIST_CMD)
1150  &&(ld->e==NULL)
1151  &&(ld->Typ()!=r->Typ()))
1152  {
1153    sleftv tmp;
1154    memset(&tmp,0,sizeof(sleftv));
1155    tmp.rtyp=DEF_CMD;
1156    b=iiAssign(&tmp,r);
1157    ld->CleanUp();
1158    memcpy(ld,&tmp,sizeof(sleftv));
1159  }
1160  else
1161  {
1162    b=iiAssign(ld,r);
1163    l->e->next=ld->e;
1164  }
1165  return b;
1166}
1167static BOOLEAN jiAssign_rec(leftv l, leftv r)
1168{
1169  leftv l1=l;
1170  leftv r1=r;
1171  leftv lrest;
1172  leftv rrest;
1173  BOOLEAN b;
1174  do
1175  {
1176    lrest=l->next;
1177    rrest=r->next;
1178    l->next=NULL;
1179    r->next=NULL;
1180    b=iiAssign(l,r);
1181    l->next=lrest;
1182    r->next=rrest;
1183    l=lrest;
1184    r=rrest;
1185  } while  ((!b)&&(l!=NULL));
1186  l1->CleanUp();
1187  r1->CleanUp();
1188  return b;
1189}
1190BOOLEAN iiAssign(leftv l, leftv r)
1191{
1192  int ll=l->listLength();
1193  int rl;
1194  int lt=l->Typ();
1195  int rt=NONE;
1196  BOOLEAN b;
1197
1198  if(l->attribute!=NULL)
1199  {
1200    if (l->rtyp==IDHDL)
1201    {
1202      atKillAll((idhdl)l->data);
1203      l->attribute=NULL;
1204    }
1205    else
1206      atKillAll((idhdl)l);
1207  }
1208  if(l->rtyp==IDHDL)
1209  {
1210    IDFLAG((idhdl)l->data)=0;
1211  }
1212  l->flag=0;
1213  if (ll==1)
1214  {
1215    /* l[..] = ... */
1216    if((l->e!=NULL)
1217    && (((l->rtyp==IDHDL) && (IDTYP((idhdl)l->data)==LIST_CMD))
1218      || (l->rtyp==LIST_CMD)))
1219    {
1220       if(r->next!=NULL)
1221         b=jiA_L_LIST(l,r);
1222       else
1223         b=jiAssign_list(l,r);
1224       if((l->rtyp==IDHDL) && (l->data!=NULL))
1225       {
1226         ipMoveId((idhdl)l->data);
1227         l->attribute=IDATTR((idhdl)l->data);
1228         l->flag=IDFLAG((idhdl)l->data);
1229       }
1230       r->CleanUp();
1231       Subexpr h;
1232       while (l->e!=NULL)
1233       {
1234         h=l->e->next;
1235         Free((ADDRESS)l->e,sizeof(*(l->e)));
1236         l->e=h;
1237       }
1238       return b;
1239    }
1240    rl=r->listLength();
1241    if (rl==1)
1242    {
1243      /* system variables = ... */
1244      if(((l->rtyp>=VECHO)&&(l->rtyp<=VPRINTLEVEL))
1245      ||((l->rtyp>=VALTVARS)&&(l->rtyp<=VMINPOLY)))
1246      {
1247        b=iiAssign_sys(l,r);
1248        r->CleanUp();
1249        //l->CleanUp();
1250        return b;
1251      }
1252      rt=r->Typ();
1253      /* a = ... */
1254      if ((lt!=MATRIX_CMD)
1255      &&(lt!=INTMAT_CMD)
1256      &&((lt==rt)||(lt!=LIST_CMD)))
1257      {
1258        b=jiAssign_1(l,r);
1259        if (l->rtyp==IDHDL)
1260        {
1261          if ((lt==DEF_CMD)||(lt=LIST_CMD)) ipMoveId((idhdl)l->data);
1262          l->attribute=IDATTR((idhdl)l->data);
1263          l->flag=IDFLAG((idhdl)l->data);
1264          l->CleanUp();
1265        }
1266        r->CleanUp();
1267        return b;
1268      }
1269      if (((lt!=LIST_CMD)
1270        &&((rt==MATRIX_CMD)
1271          ||(rt==INTMAT_CMD)
1272          ||(rt==INTVEC_CMD)
1273          ||(rt==MODUL_CMD)))
1274      ||((lt==LIST_CMD)
1275        &&(rt==RESOLUTION_CMD))
1276      )
1277      {
1278        b=jiAssign_1(l,r);
1279        if((l->rtyp==IDHDL)&&(l->data!=NULL))
1280        {
1281          if ((lt==DEF_CMD) || (lt==LIST_CMD))
1282            ipMoveId((idhdl)l->data);
1283          l->attribute=IDATTR((idhdl)l->data);
1284          l->flag=IDFLAG((idhdl)l->data);
1285        }
1286        r->CleanUp();
1287        Subexpr h;
1288        while (l->e!=NULL)
1289        {
1290          h=l->e->next;
1291          Free((ADDRESS)l->e,sizeof(*(l->e)));
1292          l->e=h;
1293        }
1294        return b;
1295      }
1296    }
1297    if (rt==NONE) rt=r->Typ();
1298  }
1299  else if (ll==(rl=r->listLength()))
1300  {
1301    b=jiAssign_rec(l,r);
1302    return b;
1303  }
1304  else
1305  {
1306    if (rt==NONE) rt=r->Typ();
1307    if (rt==INTVEC_CMD)
1308      return jiA_INTVEC_L(l,r);
1309    else if (rt==VECTOR_CMD)
1310      return jiA_VECTOR_L(l,r);
1311    else if ((rt==IDEAL_CMD)||(rt==MATRIX_CMD))
1312      return jiA_MATRIX_L(l,r);
1313    else if ((rt==STRING_CMD)&&(rl==1))
1314      return jiA_STRING_L(l,r);
1315    Werror("length of lists in assignment does not match (l:%d,r:%d)",
1316      ll,rl);
1317    return TRUE;
1318  }
1319
1320  leftv hh=r;
1321  BOOLEAN nok=FALSE;
1322  BOOLEAN map_assign=FALSE;
1323  switch (lt)
1324  {
1325    case INTVEC_CMD:
1326      nok=jjA_L_INTVEC(l,r,new intvec(exprlist_length(r)));
1327      break;
1328    case INTMAT_CMD:
1329    {
1330      nok=jjA_L_INTVEC(l,r,new intvec(IDINTVEC((idhdl)l->data)));
1331      break;
1332    }
1333    case MAP_CMD:
1334    {
1335      // first element in the list sl (r) must be a ring
1336      if (((rt == RING_CMD)||(rt == QRING_CMD))&&(r->e==NULL))
1337      {
1338        FreeL((ADDRESS)IDMAP((idhdl)l->data)->preimage);
1339        IDMAP((idhdl)l->data)->preimage = mstrdup (r->Name());
1340        /* advance the expressionlist to get the next element after the ring */
1341        hh = r->next;
1342        //r=hh;
1343      }
1344      else
1345      {
1346        WerrorS("expected ring-name");
1347        nok=TRUE;
1348        break;
1349      }
1350      if (hh==NULL) /* map-assign: map f=r; */
1351      {
1352        WerrorS("expected image ideal");
1353        nok=TRUE;
1354        break;
1355      }
1356      if ((hh->next==NULL)&&(hh->Typ()==IDEAL_CMD))
1357        return jiAssign_1(l,hh); /* map-assign: map f=r,i; */
1358      //no break, handle the rest like an ideal:
1359      map_assign=TRUE;
1360    }
1361    case MATRIX_CMD:
1362    case IDEAL_CMD:
1363    case MODUL_CMD:
1364    {
1365      sleftv t;
1366      matrix olm = (matrix)l->Data();
1367      int rk=olm->rank;
1368      char *pr=((map)olm)->preimage;
1369      BOOLEAN module_assign=(/*l->Typ()*/ lt==MODUL_CMD);
1370      matrix lm ;
1371      matrix rm;
1372      int  num;
1373      int j,k;
1374      int i=0;
1375      int mtyp=MATRIX_CMD; /*Type of left side object*/
1376      int etyp=POLY_CMD;   /*Type of elements of left side object*/
1377
1378      if (lt /*l->Typ()*/==MATRIX_CMD)
1379      {
1380        num=olm->cols()*olm->rows();
1381        lm=mpNew(olm->rows(),olm->cols());
1382      }
1383      else /* IDEAL_CMD or MODUL_CMD */
1384      {
1385        num=exprlist_length(hh);
1386        lm=(matrix)idInit(num,1);
1387        rk=1;
1388        if (module_assign)
1389        {
1390          mtyp=MODUL_CMD;
1391          etyp=VECTOR_CMD;
1392        }
1393      }
1394
1395      int ht;
1396      loop
1397      {
1398        if (hh==NULL)
1399          break;
1400        else
1401        {
1402          ht=hh->Typ();
1403          if ((j=iiTestConvert(ht,etyp))!=0)
1404          {
1405            nok=iiConvert(ht,etyp,j,hh,&t);
1406            hh->next=t.next;
1407            if (nok) break;
1408            lm->m[i]=(poly)t.CopyD(etyp);
1409            if (module_assign) rk=max(rk,pMaxComp(lm->m[i]));
1410            i++;
1411          }
1412          else
1413          if ((j=iiTestConvert(ht,mtyp))!=0)
1414          {
1415            nok=iiConvert(ht,mtyp,j,hh,&t);
1416            hh->next=t.next;
1417            if (nok) break;
1418            rm = (matrix)t.CopyD(mtyp);
1419            if (module_assign)
1420            {
1421              j = min(num,rm->cols());
1422              rk=max(rk,rm->rank);
1423            }
1424            else
1425              j = min(num-i,rm->rows() * rm->cols());
1426            for(k=0;k<j;k++,i++)
1427            {
1428              lm->m[i]=rm->m[k];
1429              rm->m[k]=NULL;
1430            }
1431            idDelete((ideal *)&rm);
1432          }
1433          else
1434          {
1435            nok=TRUE;
1436            break;
1437          }
1438          t.next=NULL;t.CleanUp();
1439          if (i==num) break;
1440          hh=hh->next;
1441        }
1442      }
1443      if (nok)
1444        idDelete((ideal *)&lm);
1445      else
1446      {
1447        idDelete((ideal *)&olm);
1448        if (module_assign)   lm->rank=rk;
1449        else if (map_assign) ((map)lm)->preimage=pr;
1450#ifdef DRING
1451        else if (pDRING)
1452        {
1453          int i=IDELEMS(lm)-1;
1454          while (i>=0)
1455          {
1456            pdSetDFlag(lm->m[i],1);
1457            i--;
1458          }
1459        }
1460#endif
1461        l=l->LData();
1462        if (l->rtyp==IDHDL)
1463          IDMATRIX((idhdl)l->data)=lm;
1464        else
1465          l->data=(char *)lm;
1466      }
1467      break;
1468    }
1469    case STRING_CMD:
1470      nok=jjA_L_STRING(l,r);
1471      break;
1472    case LIST_CMD:
1473      nok=jjA_L_LIST(l,r);
1474      break;
1475    case NONE:
1476    case 0:
1477      Werror("cannot assign to %s",l->Name());
1478      nok=TRUE;
1479      break;
1480    default:
1481      WerrorS("assign not impl.");
1482      nok=TRUE;
1483      break;
1484  } /* end switch: typ */
1485  if (nok && (!errorreported)) WerrorS("incompatible type in list assignment");
1486  r->CleanUp();
1487  return nok;
1488}
Note: See TracBrowser for help on using the repository browser.