source: git/Singular/ipassign.cc @ 358fe8e

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