Changeset 985043 in git


Ignore:
Timestamp:
May 9, 2005, 3:58:54 PM (18 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a657104b677b4c461d018cbf3204d72d34ad66a9')
Children:
9e0d84e07fb549d6dcaea4d7952cb9e3eba0b0f9
Parents:
d14343a61859db08e4e9a78b16211d1bbc44405e
Message:
*hannes: code cleanup


git-svn-id: file:///usr/local/Singular/svn/trunk@8111 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/extra.cc

    rd14343a r985043  
    1 /*****************************************
     1e*****************************************
    22*  Computer Algebra System SINGULAR      *
    33*****************************************/
    4 /* $Id: extra.cc,v 1.223 2005-05-04 14:09:45 Singular Exp $ */
     4/* $Id: extra.cc,v 1.224 2005-05-09 13:58:54 Singular Exp $ */
    55/*
    66* ABSTRACT: general interface to internals of Singular ("system" command)
     
    807807      if (h->Typ()==MATRIX_CMD)
    808808      {
    809         currRing->nc->type=nc_undef; /* to analyze later ! */
    810         //      currRing->nc->IsSkewConstant=NULL;
    811       }
    812       else
    813       {
    814         nN=pGetCoeff(pN); // pN is not NULL anyway
    815         if (nIsZero(nN))
    816         {
    817           Werror("zero coefficients are not allowed");
    818           return TRUE;
    819         }
    820         if (nIsOne(nN)) currRing->nc->type=nc_lie;
    821         else currRing->nc->type=nc_skew;
    822         currRing->nc->IsSkewConstant=1;
    823         /* create matrix C */
    824         C=mpNew(currRing->N,currRing->N);
    825         for(i=1;i<currRing->N;i++)
    826         {
    827           for(j=i+1;j<=currRing->N;j++)
    828           {
    829             MATELEM(C,i,j) = nc_p_CopyPut(pN,currRing);
    830             //  MATELEM(C,i,j)=pCopy(pN);
    831           }
    832         }
     809        currRing->nc->type=nc_undef; /* to analyze later ! */
     810        //        currRing->nc->IsSkewConstant=NULL;
     811      }
     812      else
     813      {
     814        nN=pGetCoeff(pN); // pN is not NULL anyway
     815        if (nIsZero(nN))
     816        {
     817          Werror("zero coefficients are not allowed");
     818          return TRUE;
     819        }
     820        if (nIsOne(nN)) currRing->nc->type=nc_lie;
     821        else currRing->nc->type=nc_skew;
     822        currRing->nc->IsSkewConstant=1;
     823        /* create matrix C */
     824        C=mpNew(currRing->N,currRing->N);
     825        for(i=1;i<currRing->N;i++)
     826        {
     827          for(j=i+1;j<=currRing->N;j++)
     828          {
     829            MATELEM(C,i,j) = nc_p_CopyPut(pN,currRing);
     830            //  MATELEM(C,i,j)=pCopy(pN);
     831          }
     832        }
    833833      }
    834834      pN=NULL;
     
    862862      if (pN==NULL)
    863863      {
    864         if (currRing->nc->type==nc_lie)
    865         {
    866           currRing->nc->type=nc_skew; /* even commutative! */
    867         }
    868       }
    869       else 
    870       { 
    871         if (currRing->nc->type==nc_skew) currRing->nc->type=nc_general;
     864        if (currRing->nc->type==nc_lie)
     865        {
     866          currRing->nc->type=nc_skew; /* even commutative! */
     867        }
     868      }
     869      else
     870      {
     871        if (currRing->nc->type==nc_skew) currRing->nc->type=nc_general;
    872872      } /* end pN==NULL */
    873873      if (h==NULL)
     
    951951        {
    952952          assume(pN!=NULL);
    953           if ((tmpIsSkewConstant==1) && (nIsOne(pGetCoeff(pN)))) currRing->nc->type=nc_lie;
    954           else currRing->nc->type=nc_general;
    955         }
    956         if (IsNonComm==0)
    957         {
    958           currRing->nc->type=nc_skew; /* could be also commutative */
    959           currRing->nc->IsSkewConstant=tmpIsSkewConstant;
    960         }
     953          if ((tmpIsSkewConstant==1) && (nIsOne(pGetCoeff(pN)))) currRing->nc->type=nc_lie;
     954          else currRing->nc->type=nc_general;
     955        }
     956        if (IsNonComm==0)
     957        {
     958          currRing->nc->type=nc_skew; /* could be also commutative */
     959          currRing->nc->IsSkewConstant=tmpIsSkewConstant;
     960        }
    961961      }
    962962      currRing->nc->COM=COM;
     
    10071007        h   = h->next;
    10081008      }
    1009       if ((h!=NULL)) 
    1010       {
    1011         idhdl w;
    1012         if ((w=Rop->idroot->get(h->Name(),myynest))!=NULL)
    1013         {
    1014           poly p = (poly)IDDATA(w);
    1015           res->data = pOppose(Rop,p);
    1016           res->rtyp = POLY_CMD;
    1017           return FALSE;
    1018         }
     1009      if ((h!=NULL))
     1010      {
     1011        idhdl w;
     1012        if ((w=Rop->idroot->get(h->Name(),myynest))!=NULL)
     1013        {
     1014          poly p = (poly)IDDATA(w);
     1015          res->data = pOppose(Rop,p);
     1016          res->rtyp = POLY_CMD;
     1017          return FALSE;
     1018        }
    10191019       }
    10201020      else
     
    10771077    {
    10781078      if (h == NULL || h->Typ() != INTVEC_CMD ||
    1079           h->next == NULL || h->next->Typ() != INTVEC_CMD)
     1079          h->next == NULL || h->next->Typ() != INTVEC_CMD)
    10801080      {
    10811081        WerrorS("system(\"walkAddIntVec\", intvec, intvec) expected");
     
    10911091    }
    10921092    else
    1093 #endif 
     1093#endif
    10941094#ifdef MwaklNextWeight
    10951095    if (strcmp(sys_cmd, "MwalkNextWeight") == 0)
     
    11181118      res->rtyp = INTVEC_CMD;
    11191119      res->data =  result;
    1120      
     1120
    11211121      return FALSE;
    11221122    }
     
    11371137      }
    11381138      int arg1 = (int) h->Data();
    1139    
     1139
    11401140      intvec* result = (intvec*) Mivdp(arg1);
    11411141
    11421142      res->rtyp = INTVEC_CMD;
    11431143      res->data =  result;
    1144      
    1145       return FALSE;
    1146     }   
     1144
     1145      return FALSE;
     1146    }
    11471147
    11481148    else if(strcmp(sys_cmd, "Mivlp") == 0)
     
    11601160      }
    11611161      int arg1 = (int) h->Data();
    1162    
     1162
    11631163      intvec* result = (intvec*) Mivlp(arg1);
    11641164
    11651165      res->rtyp = INTVEC_CMD;
    11661166      res->data =  result;
    1167      
     1167
    11681168      return FALSE;
    11691169    }
     
    11711171#ifdef MpDiv
    11721172      if(strcmp(sys_cmd, "MpDiv") == 0)
    1173       {       
     1173      {
    11741174        if(h==NULL || h->Typ() != POLY_CMD ||
    11751175           h->next == NULL || h->next->Typ() != POLY_CMD)
     
    11801180        poly arg1 = (poly) h->Data();
    11811181        poly arg2 = (poly) h->next->Data();
    1182        
     1182
    11831183        poly result = MpDiv(arg1, arg2);
    11841184
     
    11911191#ifdef MpMult
    11921192      if(strcmp(sys_cmd, "MpMult") == 0)
    1193       {       
     1193      {
    11941194        if(h==NULL || h->Typ() != POLY_CMD ||
    11951195           h->next == NULL || h->next->Typ() != POLY_CMD)
     
    12001200        poly arg1 = (poly) h->Data();
    12011201        poly arg2 = (poly) h->next->Data();
    1202        
     1202
    12031203        poly result = MpMult(arg1, arg2);
    12041204        res->rtyp = POLY_CMD;
     
    12171217      }
    12181218      /*
    1219       if (((intvec*) h->Data())->length() != currRing->N || 
    1220           ((intvec*) h->next->Data())->length() != currRing->N) 
    1221       {
    1222         Werror("system(\"MivSame\" ...) intvecs not of length %d\n", 
    1223                currRing->N); 
    1224         return TRUE; 
    1225       } 
     1219      if (((intvec*) h->Data())->length() != currRing->N ||
     1220          ((intvec*) h->next->Data())->length() != currRing->N)
     1221      {
     1222        Werror("system(\"MivSame\" ...) intvecs not of length %d\n",
     1223               currRing->N);
     1224        return TRUE;
     1225      }
    12261226      */
    12271227      intvec* arg1 = (intvec*) h->Data();
    12281228      intvec* arg2 = (intvec*) h->next->Data();
    1229       /*   
     1229      /*
    12301230      poly result = (poly) MivSame(arg1, arg2);
    12311231
     
    12411241      if(h == NULL || h->Typ() != INTVEC_CMD ||
    12421242         h->next == NULL || h->next->Typ() != INTVEC_CMD ||
    1243         h->next->next == NULL || h->next->next->Typ() != INTVEC_CMD  )
     1243        h->next->next == NULL || h->next->next->Typ() != INTVEC_CMD  )
    12441244      {
    12451245        Werror("system(\"M3ivSame\", intvec, intvec, intvec) expected");
     
    12471247      }
    12481248      /*
    1249       if (((intvec*) h->Data())->length() != currRing->N || 
     1249      if (((intvec*) h->Data())->length() != currRing->N ||
    12501250          ((intvec*) h->next->Data())->length() != currRing->N ||
    1251           ((intvec*) h->next->next->Data())->length() != currRing->N ) 
    1252       {
    1253         Werror("system(\"M3ivSame\" ...) intvecs not of length %d\n", 
    1254                currRing->N); 
    1255         return TRUE; 
    1256       } 
     1251          ((intvec*) h->next->next->Data())->length() != currRing->N )
     1252      {
     1253        Werror("system(\"M3ivSame\" ...) intvecs not of length %d\n",
     1254               currRing->N);
     1255        return TRUE;
     1256      }
    12571257      */
    12581258      intvec* arg1 = (intvec*) h->Data();
    12591259      intvec* arg2 = (intvec*) h->next->Data();
    12601260      intvec* arg3 = (intvec*) h->next->next->Data();
    1261       /* 
     1261      /*
    12621262      poly result = (poly) M3ivSame(arg1, arg2, arg3);
    12631263
     
    12941294    /************** Perturbation walk **********/
    12951295     if(strcmp(sys_cmd, "MivMatrixOrder") == 0)
    1296       {         
     1296      {
    12971297        if(h==NULL || h->Typ() != INTVEC_CMD)
    12981298        {
     
    13101310    else
    13111311     if(strcmp(sys_cmd, "MivMatrixOrderdp") == 0)
    1312       {         
     1312      {
    13131313        if(h==NULL || h->Typ() != INT_CMD)
    13141314        {
     
    13171317        }
    13181318        int arg1 = (int) h->Data();
    1319        
     1319
    13201320        intvec* result = (intvec*) MivMatrixOrderdp(arg1);
    13211321
     
    13271327    if(strcmp(sys_cmd, "MPertVectors") == 0)
    13281328      {
    1329        
     1329
    13301330        if(h==NULL || h->Typ() != IDEAL_CMD ||
    13311331           h->next == NULL || h->next->Typ() != INTVEC_CMD ||
     
    13351335          return TRUE;
    13361336        }
    1337        
     1337
    13381338        ideal arg1 = (ideal) h->Data();
    13391339        intvec* arg2 = (intvec*) h->next->Data();
    13401340        int arg3 = (int) h->next->next->Data();
    1341                
     1341
    13421342        intvec* result = (intvec*) MPertVectors(arg1, arg2, arg3);
    13431343
     
    13491349    if(strcmp(sys_cmd, "MPertVectorslp") == 0)
    13501350      {
    1351        
     1351
    13521352        if(h==NULL || h->Typ() != IDEAL_CMD ||
    13531353           h->next == NULL || h->next->Typ() != INTVEC_CMD ||
     
    13571357          return TRUE;
    13581358        }
    1359        
     1359
    13601360        ideal arg1 = (ideal) h->Data();
    13611361        intvec* arg2 = (intvec*) h->next->Data();
    13621362        int arg3 = (int) h->next->next->Data();
    1363                
     1363
    13641364        intvec* result = (intvec*) MPertVectorslp(arg1, arg2, arg3);
    13651365
     
    13791379        }
    13801380        ideal arg1 = (ideal) h->Data();
    1381         intvec* arg2 = (intvec*) h->next->Data(); 
     1381        intvec* arg2 = (intvec*) h->next->Data();
    13821382        intvec* result = Mfpertvector(arg1, arg2);
    1383  
     1383
    13841384        res->rtyp = INTVEC_CMD;
    13851385        res->data =  result;
     
    13881388    else
    13891389     if(strcmp(sys_cmd, "MivUnit") == 0)
    1390       {         
     1390      {
    13911391        int arg1 = (int) h->Data();
    1392        
     1392
    13931393        intvec* result = (intvec*) MivUnit(arg1);
    13941394
     
    13961396        res->data =  result;
    13971397        return FALSE;
    1398       } 
     1398      }
    13991399     else
    14001400       if(strcmp(sys_cmd, "MivWeightOrderlp") == 0)
     
    14071407        intvec* arg1 = (intvec*) h->Data();
    14081408        intvec* result = MivWeightOrderlp(arg1);
    1409  
     1409
    14101410        res->rtyp = INTVEC_CMD;
    14111411        res->data =  result;
     
    14221422        intvec* arg1 = (intvec*) h->Data();
    14231423        //int arg2 = (int) h->next->Data();
    1424  
     1424
    14251425        intvec* result = MivWeightOrderdp(arg1);
    1426  
     1426
    14271427        res->rtyp = INTVEC_CMD;
    14281428        res->data =  result;
    14291429        return FALSE;
    14301430      }
    1431     else             
     1431    else
    14321432     if(strcmp(sys_cmd, "MivMatrixOrderlp") == 0)
    1433       {         
     1433      {
    14341434        if(h==NULL || h->Typ() != INT_CMD)
    14351435        {
     
    14381438        }
    14391439        int arg1 = (int) h->Data();
    1440        
     1440
    14411441        intvec* result = (intvec*) MivMatrixOrderlp(arg1);
    14421442
     
    14711471      res->rtyp = INTVEC_CMD;
    14721472      res->data =  result;
    1473      
     1473
    14741474      return FALSE;
    14751475    }
     
    15001500      res->rtyp = INTVEC_CMD;
    15011501      res->data =  result;
    1502      
     1502
    15031503      return FALSE;
    15041504    }
     
    15221522      res->rtyp = INTVEC_CMD;
    15231523      res->data =  result;
    1524      
     1524
    15251525      return FALSE;
    15261526    }
     
    15381538
    15391539      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1540           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1540          ((intvec*) h->next->next->Data())->length() != currRing->N )
    15411541      {
    15421542        Werror("system(\"Mwalk\" ...) intvecs not of length %d\n",
     
    15531553      res->rtyp = IDEAL_CMD;
    15541554      res->data =  result;
    1555      
     1555
    15561556      return FALSE;
    15571557    }
     
    15631563          h->next == NULL || h->next->Typ() != INT_CMD ||
    15641564          h->next->next == NULL || h->next->next->Typ() != INT_CMD ||
    1565           h->next->next->next == NULL || 
    1566             h->next->next->next->Typ() != INTVEC_CMD ||
    1567           h->next->next->next->next == NULL || 
    1568             h->next->next->next->next->Typ() != INTVEC_CMD)
     1565          h->next->next->next == NULL ||
     1566            h->next->next->next->Typ() != INTVEC_CMD ||
     1567          h->next->next->next->next == NULL ||
     1568            h->next->next->next->next->Typ() != INTVEC_CMD)
    15691569      {
    15701570        Werror("system(\"Mpwalk\", ideal, int, int, intvec, intvec) expected");
     
    15731573
    15741574      if (((intvec*) h->next->next->next->Data())->length() != currRing->N &&
    1575           ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
     1575          ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
    15761576      {
    15771577        Werror("system(\"Mpwalk\" ...) intvecs not of length %d\n",
     
    15901590      res->rtyp = IDEAL_CMD;
    15911591      res->data =  result;
    1592      
     1592
    15931593      return FALSE;
    15941594    }
     
    16001600          h->next == NULL || h->next->Typ() != INT_CMD ||
    16011601          h->next->next == NULL || h->next->next->Typ() != INT_CMD ||
    1602           h->next->next->next == NULL || 
    1603             h->next->next->next->Typ() != INTVEC_CMD ||
    1604           h->next->next->next->next == NULL || 
    1605             h->next->next->next->next->Typ() != INTVEC_CMD||
    1606           h->next->next->next->next->next == NULL || 
    1607             h->next->next->next->next->next->Typ() != INT_CMD)
     1602          h->next->next->next == NULL ||
     1603            h->next->next->next->Typ() != INTVEC_CMD ||
     1604          h->next->next->next->next == NULL ||
     1605            h->next->next->next->next->Typ() != INTVEC_CMD||
     1606          h->next->next->next->next->next == NULL ||
     1607            h->next->next->next->next->next->Typ() != INT_CMD)
    16081608      {
    16091609        Werror("system(\"Mpwalk\", ideal, int, int, intvec, intvec, int) expected");
     
    16121612
    16131613      if (((intvec*) h->next->next->next->Data())->length() != currRing->N &&
    1614           ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
     1614          ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
    16151615      {
    16161616        Werror("system(\"Mpwalk\" ...) intvecs not of length %d\n",
     
    16301630      res->rtyp = IDEAL_CMD;
    16311631      res->data =  result;
    1632      
     1632
    16331633      return FALSE;
    16341634    }
     
    16391639          h->next == NULL || h->next->Typ() != INT_CMD ||
    16401640          h->next->next == NULL || h->next->next->Typ() != INT_CMD ||
    1641           h->next->next->next == NULL || 
    1642             h->next->next->next->Typ() != INTVEC_CMD ||
    1643           h->next->next->next->next == NULL || 
    1644             h->next->next->next->next->Typ() != INTVEC_CMD)
     1641          h->next->next->next == NULL ||
     1642            h->next->next->next->Typ() != INTVEC_CMD ||
     1643          h->next->next->next->next == NULL ||
     1644            h->next->next->next->next->Typ() != INTVEC_CMD)
    16451645      {
    16461646        Werror("system(\"MAltwalk1\", ideal, int, int, intvec, intvec) expected");
     
    16491649
    16501650      if (((intvec*) h->next->next->next->Data())->length() != currRing->N &&
    1651           ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
     1651          ((intvec*) h->next->next->next->next->Data())->length()!=currRing->N)
    16521652      {
    16531653        Werror("system(\"MAltwalk1\" ...) intvecs not of length %d\n",
     
    16661666      res->rtyp = IDEAL_CMD;
    16671667      res->data =  result;
    1668      
     1668
    16691669      return FALSE;
    16701670    }
     
    16831683
    16841684      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1685           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1685          ((intvec*) h->next->next->Data())->length() != currRing->N )
    16861686      {
    16871687        Werror("system(\"Mfwalk\" ...) intvecs not of length %d\n",
     
    16981698      res->rtyp = IDEAL_CMD;
    16991699      res->data =  result;
    1700      
     1700
    17011701      return FALSE;
    17021702    }
     
    17141714
    17151715      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1716           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1716          ((intvec*) h->next->next->Data())->length() != currRing->N )
    17171717      {
    17181718        Werror("system(\"Mfwalk\" ...) intvecs not of length %d\n",
     
    17281728      res->rtyp = IDEAL_CMD;
    17291729      res->data =  result;
    1730      
    1731       return FALSE;
    1732     }
    1733     else
    1734  
     1730
     1731      return FALSE;
     1732    }
     1733    else
     1734
    17351735#ifdef TRAN_Orig
    17361736    if (strcmp(sys_cmd, "TranMImprovwalk") == 0)
     
    17451745
    17461746      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1747           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1747          ((intvec*) h->next->next->Data())->length() != currRing->N )
    17481748      {
    17491749        Werror("system(\"TranMImprovwalk\" ...) intvecs not of length %d\n",
     
    17601760      res->rtyp = IDEAL_CMD;
    17611761      res->data =  result;
    1762      
     1762
    17631763      return FALSE;
    17641764    }
     
    17761776
    17771777      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1778           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1778          ((intvec*) h->next->next->Data())->length() != currRing->N )
    17791779      {
    17801780        Werror("system(\"MAltwalk2\" ...) intvecs not of length %d\n",
     
    17911791      res->rtyp = IDEAL_CMD;
    17921792      res->data =  result;
    1793      
     1793
    17941794      return FALSE;
    17951795    }
     
    18071807
    18081808      if (((intvec*) h->next->Data())->length() != currRing->N &&
    1809           ((intvec*) h->next->next->Data())->length() != currRing->N )
     1809          ((intvec*) h->next->next->Data())->length() != currRing->N )
    18101810      {
    18111811        Werror("system(\"TranMImprovwalk\" ...) intvecs not of length %d\n",
     
    18221822      res->rtyp = IDEAL_CMD;
    18231823      res->data =  result;
    1824      
     1824
    18251825      return FALSE;
    18261826    }
     
    24792479#endif
    24802480/*==================== t-rep-GB ==================================*/
    2481     if (strcmp(sys_cmd, "trepgb")==0)
    2482     {
    2483       ring r = currRing;
    2484       ideal i = (ideal)h->Data();
    2485       h=h->next;
    2486       if(h)
    2487       {
    2488        
    2489         BOOLEAN b=(BOOLEAN) h->Data();
    2490         res->data=t_rep_gb(r,i,b);
    2491       }
    2492       else
    2493         res->data=t_rep_gb(r,i);
    2494       res->rtyp=IDEAL_CMD;
    2495       setFlag(res,FLAG_STD);
    2496       return(FALSE);
    2497     }
    2498     else
    2499       if (strcmp(sys_cmd, "unifastmult")==0)
     2481    if (strcmp(sys_cmd, "unifastmult")==0)
    25002482    {
    25012483      ring r = currRing;
     
    25082490    }
    25092491    else
    2510       if (strcmp(sys_cmd, "multifastmult")==0)
     2492    if (strcmp(sys_cmd, "multifastmult")==0)
    25112493    {
    25122494      ring r = currRing;
     
    25182500      return(FALSE);
    25192501    }
    2520       else
    2521 
    2522       if (strcmp(sys_cmd, "mults")==0)
    2523     {
    2524 
     2502    else
     2503    if (strcmp(sys_cmd, "mults")==0)
     2504    {
    25252505      res->rtyp=INT_CMD ;
    25262506      res->data=(void*) Mults();
    25272507      return(FALSE);
    25282508    }
    2529       else
    2530         if (strcmp(sys_cmd, "fastpower")==0)
     2509    else
     2510    if (strcmp(sys_cmd, "fastpower")==0)
    25312511    {
    25322512      ring r = currRing;
     
    25382518      return(FALSE);
    25392519    }
    2540         else
    2541                 if (strcmp(sys_cmd, "normalpower")==0)
     2520    else
     2521    if (strcmp(sys_cmd, "normalpower")==0)
    25422522    {
    25432523      ring r = currRing;
     
    25502530    }
    25512531    else
    2552      
    2553                 if (strcmp(sys_cmd, "MCpower")==0)
     2532    if (strcmp(sys_cmd, "MCpower")==0)
    25542533    {
    25552534      ring r = currRing;
     
    25622541    }
    25632542    else
    2564      
    2565                 if (strcmp(sys_cmd, "bit_subst")==0)
     2543    if (strcmp(sys_cmd, "bit_subst")==0)
    25662544    {
    25672545      ring r = currRing;
Note: See TracChangeset for help on using the changeset viewer.