Ignore:
Timestamp:
Apr 21, 2011, 1:34:10 PM (13 years ago)
Author:
Oleksandr Motsak <motsak@…>
Branches:
(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')
Children:
59a62bd935a023662ba918abaefcc861c4adbe81
Parents:
0b7bf7ebce7bb8127ebbb274e87efef8dce670c7
git-author:
Oleksandr Motsak <motsak@mathematik.uni-kl.de>2011-04-21 13:34:10+02:00
git-committer:
Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 12:31:24+01:00
Message:
FIX: moved P/minpoly/parameters from ring to coeff
File:
1 edited

Legend:

Unmodified
Added
Removed

  • libpolys/polys/monomials/ring.cc

    r0b7bf7 r01c1d0  
    7474};
    7575
    76 #define  rInternalChar(r) ((r)->ch)
     76#define  rInternalChar(r) ((r)->cf->ch)
    7777
    7878const char * rSimpleOrdStr(int ord)
    … …  
    9494//BOOLEAN rField_is_R(ring r)
    9595//{
    96 //  if (r->ch== -1)
     96//  if (r->cf->ch== -1)
    9797//  {
    9898//    if (r->float_len==(short)0) return TRUE;
    … …  
    108108void rChangeCurrRing(ring r)
    109109{
    110  // if ((currRing!=NULL) && (currRing->minpoly!=NULL))
     110 // if ((currRing!=NULL) && (currRing->cf->minpoly!=NULL))
    111111 // {
    112  //   omCheckAddr(currRing->minpoly);
     112 //   omCheckAddr(currRing->cf->minpoly);
    113113 // }
    114114  //------------ set global ring vars --------------------------------
    … …  
    140140  assume( cf != NULL);
    141141  ring r=(ring) omAlloc0Bin(sip_sring_bin);
    142   r->ch    = n_GetChar(cf);
     142  r->cf->ch    = n_GetChar(cf);
    143143  r->N     = N;
    144144  r->cf = cf;
    145   /*r->P     = 0; Alloc0 */
     145  /*rPar(r)  = 0; Alloc0 */
    146146  /*names*/
    147147  r->names = (char **) omAlloc0(N * sizeof(char *));
    … …  
    279279  nblocks--;
    280280
    281 
     281  n_CoeffWrite(r->cf);
     282#if 0
    282283  if (rField_is_GF(r))
    283284  {
    284     Print("//   # ground field : %d\n",r->ch);
     285    Print("//   # ground field : %d\n",r->cf->ch);
    285286    Print("//   primitive element : %s\n", r->parameter[0]);
    286287    StringSetS("//   minpoly        : ");
    … …  
    313314    else
    314315      Print ("%d\n",rChar(r)); /* Fp(a) */
    315     if (r->parameter!=NULL)
     316    if (r->cf->parameter!=NULL)
    316317    {
    317318      Print ("//   %d parameter    : ",rPar(r));
    318       char **sp=r->parameter;
     319      char **sp=r->cf->parameter;
    319320      int nop=0;
    320321      while (nop<rPar(r))
    … …  
    328329      {
    329330        // i^2+1:
    330         Print("(%s^2+1)\n",r->parameter[0]);
    331       }
    332       else if (r->minpoly==NULL)
     331        Print("(%s^2+1)\n",r->cf->parameter[0]);
     332      }
     333      else if (r->cf->minpoly==NULL)
    333334      {
    334335        PrintS("0\n");
    … …  
    336337      else
    337338      {
    338         StringSetS(""); n_Write(r->minpoly,r->cf); PrintS(StringAppendS("\n"));
     339        StringSetS(""); n_Write(r->cf->minpoly,r->cf); PrintS(StringAppendS("\n"));
    339340      }
    340341      //if (r->minideal!=NULL)
    … …  
    345346    }
    346347  }
     348#endif 
    347349  Print("//   number of vars : %d",r->N);
    348350
    … …  
    477479#endif
    478480
    479   nKillChar(r->cf);
     481  nKillChar(r->cf); r->cf = NULL;
    480482  rUnComplete(r);
    481483  // delete order stuff
    … …  
    512514
    513515  // delete parameter
    514   if (r->parameter!=NULL)
    515   {
    516     char **s=r->parameter;
     516  if (r->cf->parameter!=NULL)
     517  {
     518    char **s=r->cf->parameter;
    517519    j = 0;
    518520    while (j < rPar(r))
    … …  
    522524      j++;
    523525    }
    524     omFreeSize((ADDRESS)r->parameter,rPar(r)*sizeof(char *));
    525   }
    526 #ifdef HAVE_RINGS
    527   if (r->cf->modBase != NULL)
    528   {
    529     mpz_clear(r->cf->modBase);
    530     omFree((ADDRESS)r->cf->modBase);
    531   }
    532   if (r->cf->modNumber != NULL)
    533   {
    534     mpz_clear(r->cf->modNumber);
    535     omFree((ADDRESS)r->cf->modNumber);
    536   }
    537 #endif
     526    omFreeSize((ADDRESS)r->cf->parameter,rPar(r)*sizeof(char *));
     527  }
    538528  omFreeBin(r, sip_sring_bin);
    539529}
    … …  
    679669  }
    680670#endif
    681   if (r->parameter==NULL)
    682   {
    683     i=r->ch;
     671  if (r->cf->parameter==NULL)
     672  {
     673    i=r->cf->ch;
    684674    if(i==-1)
    685675      s=omStrDup("real");                    /* R */
    … …  
    693683  if (rField_is_long_C(r))
    694684  {
    695     s=(char *)omAlloc(21+strlen(r->parameter[0]));
    696     sprintf(s,"complex,%d,%s",r->float_len,r->parameter[0]);   /* C */
     685    s=(char *)omAlloc(21+strlen(r->cf->parameter[0]));
     686    sprintf(s,"complex,%d,%s",r->float_len,r->cf->parameter[0]);   /* C */
    697687    return s;
    698688  }
    … …  
    700690  for(i=0; i<rPar(r);i++)
    701691  {
    702     l+=(strlen(r->parameter[i])+1);
     692    l+=(strlen(r->cf->parameter[i])+1);
    703693  }
    704694  s=(char *)omAlloc((long)(l+MAX_INT_LEN+1));
    705695  s[0]='\0';
    706   if (r->ch<0)       sprintf(s,"%d",-r->ch); /* Fp(a) */
    707   else if (r->ch==1) sprintf(s,"0");         /* Q(a)  */
     696  if (r->cf->ch<0)       sprintf(s,"%d",-r->cf->ch); /* Fp(a) */
     697  else if (r->cf->ch==1) sprintf(s,"0");         /* Q(a)  */
    708698  else
    709699  {
    710     sprintf(s,"%d,%s",r->ch,r->parameter[0]); /* GF(q)  */
     700    sprintf(s,"%d,%s",r->cf->ch,r->cf->parameter[0]); /* GF(q)  */
    711701    return s;
    712702  }
    … …  
    717707  {
    718708    strcat(s,tt);
    719     strcat(s,r->parameter[i]);
     709    strcat(s,r->cf->parameter[i]);
    720710  }
    721711  return s;
    … …  
    724714char * rParStr(ring r)
    725715{
    726   if ((r==NULL)||(r->parameter==NULL)) return omStrDup("");
     716  if ((r==NULL)||(r->cf->parameter==NULL)) return omStrDup("");
    727717
    728718  int i;
    … …  
    731721  for (i=0; i<rPar(r); i++)
    732722  {
    733     l+=strlen(r->parameter[i])+1;
     723    l+=strlen(r->cf->parameter[i])+1;
    734724  }
    735725  char *s=(char *)omAlloc((long)l);
    … …  
    737727  for (i=0; i<rPar(r)-1; i++)
    738728  {
    739     strcat(s,r->parameter[i]);
     729    strcat(s,r->cf->parameter[i]);
    740730    strcat(s,",");
    741731  }
    742   strcat(s,r->parameter[i]);
     732  strcat(s,r->cf->parameter[i]);
    743733  return s;
    744734}
    … …  
    759749int  rIsExtension(const ring r)
    760750{
    761   return (r->parameter!=NULL); /* R, Q, Fp: FALSE */
     751  return (r->cf->parameter!=NULL); /* R, Q, Fp: FALSE */
    762752}
    763753
    … …  
    793783    return 0;
    794784  if (!rIsExtension(r)) /* Q, Fp */
    795     return r->ch;
     785    return r->cf->ch;
    796786  if (rField_is_Zp_a(r))  /* Fp(a)  */
    797     return -r->ch;
     787    return -r->cf->ch;
    798788  if (rField_is_Q_a(r))   /* Q(a)  */
    799789    return 0;
    800790  /*else*/               /* GF(p,n) */
    801791  {
    802     if ((r->ch & 1)==0) return 2;
     792    if ((r->cf->ch & 1)==0) return 2;
    803793    int i=3;
    804     while ((r->ch % i)!=0) i+=2;
     794    while ((r->cf->ch % i)!=0) i+=2;
    805795    return i;
    806796  }
    … …  
    843833int rSumInternal(ring r1, ring r2, ring &sum, BOOLEAN vartest, BOOLEAN dp_dp)
    844834{
     835
    845836  ip_sring tmpR;
    846837  memset(&tmpR,0,sizeof(tmpR));
    847838  /* check coeff. field =====================================================*/
     839
     840  TODO(Hans, switch to the use of ->cf);
     841  assume( tmpR.cf !=NULL ); // ???!
     842 
    848843  if ((rFieldType(r1)==rFieldType(r2))
    849844  && (rInternalChar(r1)==rInternalChar(r2)))
    850845  {
    851     tmpR.ch=r1->ch;
    852846    if (rField_is_Q(r1)||rField_is_Zp(r1)||rField_is_GF(r1)) /*Q, Z/p, GF(p,n)*/
    853847    {
    854       if (r1->parameter!=NULL)
    855       {
    856         if (!vartest || (strcmp(r1->parameter[0],r2->parameter[0])==0)) /* 1 par */
     848      if (r1->cf->parameter!=NULL)
     849      {
     850        if (!vartest || (strcmp(r1->cf->parameter[0],r2->cf->parameter[0])==0)) /* 1 par */
    857851        {
    858852          //tmpR.parameter=(char **)omAllocBin(char *);
    859           //tmpR.parameter[0]=omStrDup(r1->parameter[0]);
     853          //tmpR.parameter[0]=omStrDup(r1->cf->parameter[0]);
    860854          //tmpR.P=1;
    861855        }
    … …  
    866860        }
    867861      }
     862      tmpR.cf->ch=r1->cf->ch;
    868863    }
    869864    else if (rField_is_Extension(r1)) /* Q(a),Z/p(a) */
    870865    {
    871       if (r1->minpoly!=NULL)
    872       {
    873         if (r2->minpoly!=NULL)
     866      if (r1->cf->minpoly!=NULL)
     867      {
     868        if (r2->cf->minpoly!=NULL)
    874869        {
    875870          // HANNES: TODO: delete nSetChar
    … …  
    878873//        rChangeCurrRing(r1);
    879874           
    880           if ((!vartest || (strcmp(r1->parameter[0],r2->parameter[0])==0)) /* 1 par */
    881               && n_Equal(r1->minpoly,r2->minpoly, r1->cf))
     875          if ((!vartest || (strcmp(r1->cf->parameter[0],r2->cf->parameter[0])==0)) /* 1 par */
     876              && n_Equal(r1->cf->minpoly,r2->cf->minpoly, r1->cf))
    882877          {
    883878            //tmpR.parameter=(char **)omAllocBin(char *);
    884             //tmpR.parameter[0]=omStrDup(r1->parameter[0]);
    885             //tmpR.minpoly=n_Copy(r1->minpoly, r1->cf);
     879            //tmpR.parameter[0]=omStrDup(r1->cf->parameter[0]);
     880            //tmpR.minpoly=n_Copy(r1->cf->minpoly, r1->cf);
    886881            //tmpR.P=1;
    887882            // HANNES: TODO: delete nSetChar
    … …  
    896891        else
    897892        {
    898           if ((!vartest || (strcmp(r1->parameter[0],r2->parameter[0])==0)) /* 1 par */
     893          if ((!vartest || (strcmp(r1->cf->parameter[0],r2->cf->parameter[0])==0)) /* 1 par */
    899894              && (rPar(r2)==1))
    900895          {
    901896            //tmpR.parameter=(char **)omAllocBin(char *);
    902             //tmpR.parameter[0]=omStrDup(r1->parameter[0]);
     897            //tmpR.parameter[0]=omStrDup(r1->cf->parameter[0]);
    903898            //tmpR.P=1;
    904             //tmpR.minpoly=n_Copy(r1->minpoly, r1->cf);
     899            //tmpR.minpoly=n_Copy(r1->cf->minpoly, r1->cf);
    905900          }
    906901          else
    … …  
    911906        }
    912907      }
    913       else /* r1->minpoly==NULL */
    914       {
    915         if (r2->minpoly!=NULL)
    916         {
    917           if ((!vartest || (strcmp(r1->parameter[0],r2->parameter[0])==0)) /* 1 par */
     908      else /* r1->cf->minpoly==NULL */
     909      {
     910        if (r2->cf->minpoly!=NULL)
     911        {
     912          if ((!vartest || (strcmp(r1->cf->parameter[0],r2->cf->parameter[0])==0)) /* 1 par */
    918913              && (rPar(r1)==1))
    919914          {
    920915            //tmpR.parameter=(char **)omAllocBin(char *);
    921             //tmpR.parameter[0]=omStrDup(r1->parameter[0]);
     916            //tmpR.parameter[0]=omStrDup(r1->cf->parameter[0]);
    922917            //tmpR.P=1;
    923             //tmpR.minpoly=n_Copy(r2->minpoly, r2->cf);
     918            //tmpR.minpoly=n_Copy(r2->cf->minpoly, r2->cf);
    924919          }
    925920          else
    … …  
    932927        {
    933928          int len=rPar(r1)+rPar(r2);
    934           tmpR.parameter=(char **)omAlloc0(len*sizeof(char *));
     929          tmpR.cf->parameter=(char **)omAlloc0(len*sizeof(char *));
    935930          int i;
    936931          for (i=0;i<rPar(r1);i++)
    937932          {
    938             tmpR.parameter[i]=omStrDup(r1->parameter[i]);
     933            tmpR.cf->parameter[i]=omStrDup(r1->cf->parameter[i]);
    939934          }
    940935          int j,l;
    … …  
    945940              for(l=0;l<i;l++)
    946941              {
    947                 if(strcmp(tmpR.parameter[l],r2->parameter[j])==0)
     942                if(strcmp(tmpR.cf->parameter[l],r2->cf->parameter[j])==0)
    948943                  break;
    949944              }
    … …  
    953948            if (l==i)
    954949            {
    955               tmpR.parameter[i]=omStrDup(r2->parameter[j]);
     950              tmpR.cf->parameter[i]=omStrDup(r2->cf->parameter[j]);
    956951              i++;
    957952            }
    … …  
    959954          if (i!=len)
    960955          {
    961             tmpR.parameter=(char**)omReallocSize(tmpR.parameter,
     956            tmpR.cf->parameter=(char**)omReallocSize(tmpR.cf->parameter,
    962957                                                 len*sizeof(char *),
    963958                                                 i*sizeof(char *));
    964959          }
    965           tmpR.P=i;
     960          tmpR.cf->P=i;
    966961        }
    967962      }
    … …  
    998993    #endif
    999994  }
    1000   else /* r1->ch!=r2->ch */
    1001   {
    1002     if (r1->ch<-1) /* Z/p(a) */
    1003     {
    1004       if ((r2->ch==0) /* Q */
    1005           || (r2->ch==-r1->ch)) /* Z/p */
    1006       {
    1007         tmpR.ch=r1->ch;
    1008         tmpR.P=rPar(r1);
    1009         tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
     995  else /* r1->cf->ch!=r2->cf->ch */
     996  {
     997    if (r1->cf->ch<-1) /* Z/p(a) */
     998    {
     999      if ((r2->cf->ch==0) /* Q */
     1000          || (r2->cf->ch==-r1->cf->ch)) /* Z/p */
     1001      {
     1002        tmpR.cf->ch=r1->cf->ch;
     1003        tmpR.cf->P=rPar(r1);
     1004        tmpR.cf->parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
    10101005        int i;
    10111006        for (i=0;i<rPar(r1);i++)
    10121007        {
    1013           tmpR.parameter[i]=omStrDup(r1->parameter[i]);
    1014         }
    1015         if (r1->minpoly!=NULL)
    1016         {
    1017           tmpR.minpoly=n_Copy(r1->minpoly, r1->cf);
     1008          tmpR.cf->parameter[i]=omStrDup(r1->cf->parameter[i]);
     1009        }
     1010        if (r1->cf->minpoly!=NULL)
     1011        {
     1012          tmpR.cf->minpoly=n_Copy(r1->cf->minpoly, r1->cf);
    10181013        }
    10191014      }
    … …  
    10241019      }
    10251020    }
    1026     else if (r1->ch==-1) /* R */
     1021    else if (r1->cf->ch==-1) /* R */
    10271022    {
    10281023      WerrorS("R+..");
    10291024      return -1;
    10301025    }
    1031     else if (r1->ch==0) /* Q */
    1032     {
    1033       if ((r2->ch<-1)||(r2->ch==1)) /* Z/p(a),Q(a) */
    1034       {
    1035         tmpR.ch=r2->ch;
    1036         tmpR.P=rPar(r2);
    1037         tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
     1026    else if (r1->cf->ch==0) /* Q */
     1027    {
     1028      if ((r2->cf->ch<-1)||(r2->cf->ch==1)) /* Z/p(a),Q(a) */
     1029      {
     1030        tmpR.cf->ch=r2->cf->ch;
     1031        tmpR.cf->P=rPar(r2);
     1032        tmpR.cf->parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
    10381033        int i;
    10391034        for (i=0;i<rPar(r2);i++)
    10401035        {
    1041           tmpR.parameter[i]=omStrDup(r2->parameter[i]);
    1042         }
    1043         if (r2->minpoly!=NULL)
    1044         {
    1045           tmpR.minpoly=n_Copy(r2->minpoly, r2->cf);
    1046         }
    1047       }
    1048       else if (r2->ch>1) /* Z/p,GF(p,n) */
    1049       {
    1050         tmpR.ch=r2->ch;
    1051         if (r2->parameter!=NULL)
    1052         {
    1053           tmpR.parameter=(char_ptr *)omAllocBin(char_ptr_bin);
    1054           tmpR.P=1;
    1055           tmpR.parameter[0]=omStrDup(r2->parameter[0]);
     1036          tmpR.cf->parameter[i]=omStrDup(r2->cf->parameter[i]);
     1037        }
     1038        if (r2->cf->minpoly!=NULL)
     1039        {
     1040          tmpR.cf->minpoly=n_Copy(r2->cf->minpoly, r2->cf);
     1041        }
     1042      }
     1043      else if (r2->cf->ch>1) /* Z/p,GF(p,n) */
     1044      {
     1045        tmpR.cf->ch=r2->cf->ch;
     1046        if (r2->cf->parameter!=NULL)
     1047        {
     1048          tmpR.cf->parameter=(char_ptr *)omAllocBin(char_ptr_bin);
     1049          tmpR.cf->P=1;
     1050          tmpR.cf->parameter[0]=omStrDup(r2->cf->parameter[0]);
    10561051        }
    10571052      }
    … …  
    10621057      }
    10631058    }
    1064     else if (r1->ch==1) /* Q(a) */
    1065     {
    1066       if (r2->ch==0) /* Q */
    1067       {
    1068         tmpR.ch=r1->ch;
    1069         tmpR.P=rPar(r1);
    1070         tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
     1059    else if (r1->cf->ch==1) /* Q(a) */
     1060    {
     1061      if (r2->cf->ch==0) /* Q */
     1062      {
     1063        tmpR.cf->ch=r1->cf->ch;
     1064        tmpR.cf->P=rPar(r1);
     1065        tmpR.cf->parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
    10711066        int i;
    10721067        for(i=0;i<rPar(r1);i++)
    10731068        {
    1074           tmpR.parameter[i]=omStrDup(r1->parameter[i]);
    1075         }
    1076         if (r1->minpoly!=NULL)
    1077         {
    1078           tmpR.minpoly=n_Copy(r1->minpoly, r1->cf);
     1069          tmpR.cf->parameter[i]=omStrDup(r1->cf->parameter[i]);
     1070        }
     1071        if (r1->cf->minpoly!=NULL)
     1072        {
     1073          tmpR.cf->minpoly=n_Copy(r1->cf->minpoly, r1->cf);
    10791074        }
    10801075      }
    … …  
    10851080      }
    10861081    }
    1087     else /* r1->ch >=2 , Z/p */
    1088     {
    1089       if (r2->ch==0) /* Q */
    1090       {
    1091         tmpR.ch=r1->ch;
    1092       }
    1093       else if (r2->ch==-r1->ch) /* Z/p(a) */
    1094       {
    1095         tmpR.ch=r2->ch;
    1096         tmpR.P=rPar(r2);
    1097         tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
     1082    else /* r1->cf->ch >=2 , Z/p */
     1083    {
     1084      if (r2->cf->ch==0) /* Q */
     1085      {
     1086        tmpR.cf->ch=r1->cf->ch;
     1087      }
     1088      else if (r2->cf->ch==-r1->cf->ch) /* Z/p(a) */
     1089      {
     1090        tmpR.cf->ch=r2->cf->ch;
     1091        tmpR.cf->P=rPar(r2);
     1092        tmpR.cf->parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
    10981093        int i;
    10991094        for(i=0;i<rPar(r2);i++)
    11001095        {
    1101           tmpR.parameter[i]=omStrDup(r2->parameter[i]);
    1102         }
    1103         if (r2->minpoly!=NULL)
    1104         {
    1105           tmpR.minpoly=n_Copy(r2->minpoly, r2->cf);
     1096          tmpR.cf->parameter[i]=omStrDup(r2->cf->parameter[i]);
     1097        }
     1098        if (r2->cf->minpoly!=NULL)
     1099        {
     1100          tmpR.cf->minpoly=n_Copy(r2->cf->minpoly, r2->cf);
    11061101        }
    11071102      }
    … …  
    11271122    if (*(r1->names[i]) == '\0')
    11281123      b = FALSE;
    1129     else if ((r2->parameter!=NULL) && (strlen(r1->names[i])==1))
     1124    else if ((r2->cf->parameter!=NULL) && (strlen(r1->names[i])==1))
    11301125    {
    11311126      if (vartest)
    … …  
    11331128        for(j=0;j<rPar(r2);j++)
    11341129        {
    1135           if (strcmp(r1->names[i],r2->parameter[j])==0)
     1130          if (strcmp(r1->names[i],r2->cf->parameter[j])==0)
    11361131          {
    11371132            b=FALSE;
    … …  
    11591154    if (*(r2->names[i]) == '\0')
    11601155      b = FALSE;
    1161     else if ((r1->parameter!=NULL) && (strlen(r2->names[i])==1))
     1156    else if ((r1->cf->parameter!=NULL) && (strlen(r2->names[i])==1))
    11621157    {
    11631158      if (vartest)
    … …  
    11651160        for(j=0;j<rPar(r1);j++)
    11661161        {
    1167           if (strcmp(r2->names[i],r1->parameter[j])==0)
     1162          if (strcmp(r2->names[i],r1->cf->parameter[j])==0)
    11681163          {
    11691164            b=FALSE;
    … …  
    14211416      if (rPar(R2)!=0) par_perm2=(int *)omAlloc0((rPar(R2)+1)*sizeof(int));
    14221417
    1423       maFindPerm(R1->names,  rVar(R1),  R1->parameter,  rPar(R1),
    1424                  sum->names, rVar(sum), sum->parameter, rPar(sum),
    1425                  perm1, par_perm1, sum->ch);
    1426 
    1427       maFindPerm(R2->names,  rVar(R2),  R2->parameter,  rPar(R2),
    1428                  sum->names, rVar(sum), sum->parameter, rPar(sum),
    1429                  perm2, par_perm2, sum->ch);
     1418      maFindPerm(R1->names,  rVar(R1),  R1->cf->parameter,  rPar(R1),
     1419                 sum->names, rVar(sum), sum->cf->parameter, rPar(sum),
     1420                 perm1, par_perm1, sum->cf->ch);
     1421
     1422      maFindPerm(R2->names,  rVar(R2),  R2->cf->parameter,  rPar(R2),
     1423                 sum->names, rVar(sum), sum->cf->parameter, rPar(sum),
     1424                 perm2, par_perm2, sum->cf->ch);
    14301425
    14311426
    … …  
    15261521    int *par_perm1 = NULL;
    15271522    if (rPar(r1)!=0) par_perm1=(int *)omAlloc0((rPar(r1)+1)*sizeof(int));
    1528     maFindPerm(r1->names,  rVar(r1),  r1->parameter,  rPar(r1),
    1529                sum->names, rVar(sum), sum->parameter, rPar(sum),
    1530                perm1, par_perm1, sum->ch);
     1523    maFindPerm(r1->names,  rVar(r1),  r1->cf->parameter,  rPar(r1),
     1524               sum->names, rVar(sum), sum->cf->parameter, rPar(sum),
     1525               perm1, par_perm1, sum->cf->ch);
    15311526    nMapFunc nMap1 = n_SetMap(r1->cf,sum->cf);
    15321527    Q1 = idInit(IDELEMS(r1->qideal),1);
    … …  
    15491544    int *par_perm2 = NULL;
    15501545    if (rPar(r2)!=0) par_perm2=(int *)omAlloc0((rPar(r2)+1)*sizeof(int));
    1551     maFindPerm(r2->names,  rVar(r2),  r2->parameter,  rPar(r2),
    1552                sum->names, rVar(sum), sum->parameter, rPar(sum),
    1553                perm2, par_perm2, sum->ch);
     1546    maFindPerm(r2->names,  rVar(r2),  r2->cf->parameter,  rPar(r2),
     1547               sum->names, rVar(sum), sum->cf->parameter, rPar(sum),
     1548               perm2, par_perm2, sum->cf->ch);
    15541549    nMapFunc nMap2 = n_SetMap(r2->cf,sum->cf);
    15551550    Q2 = idInit(IDELEMS(r2->qideal),1);
    … …  
    16201615  //struct omBin   PolyBin; /* Bin from where monoms are allocated */
    16211616  //memset: res->PolyBin=NULL; // rComplete
    1622   res->ch=r->ch;     /* characteristic */
     1617  res->cf->ch=r->cf->ch;     /* characteristic */
    16231618#ifdef HAVE_RINGS
    16241619  res->cf->ringtype = r->cf->ringtype;  // cring = 0 => coefficient field, cring = 1 => coeffs from Z/2^m
    … …  
    16401635  res->float_len2=r->float_len2; /* additional char-flags */
    16411636
    1642   res->N=r->N;      /* number of vars */
    1643   res->P=r->P;      /* number of pars */
     1637  res->N=rVar(r);      /* number of vars */
     1638  res->cf->P=rPar(r);      /* number of pars */
    16441639  res->OrdSgn=r->OrdSgn; /* 1 for polynomial rings, -1 otherwise */
    16451640
    … …  
    16971692  res->algring=r->algring;
    16981693  //memset: res->minideal=NULL;
    1699   if (r->parameter!=NULL)
    1700   {
    1701     if (r->minpoly!=NULL) res->minpoly=n_Copy(r->minpoly,r->cf);
     1694  if (r->cf->parameter!=NULL)
     1695  {
     1696    if (r->cf->minpoly!=NULL) res->cf->minpoly=n_Copy(r->cf->minpoly,r->cf);
    17021697    int l=rPar(r);
    1703     res->parameter=(char **)omAlloc(l*sizeof(char *));
     1698    res->cf->parameter=(char **)omAlloc(l*sizeof(char *));
    17041699    int i;
    17051700    for(i=0;i<rPar(r);i++)
    17061701    {
    1707       res->parameter[i]=omStrDup(r->parameter[i]);
     1702      res->cf->parameter[i]=omStrDup(r->cf->parameter[i]);
    17081703    }
    17091704    if (r->minideal!=NULL)
    … …  
    18451840  for (i=0; i<rPar(r1);i++)
    18461841  {
    1847       if (strcmp(r1->parameter[i], r2->parameter[i])!=0)
     1842      if (strcmp(r1->cf->parameter[i], r2->cf->parameter[i])!=0)
    18481843        return FALSE;
    18491844  }
    18501845
    1851   if (r1->minpoly != NULL)
    1852   {
    1853     if (r2->minpoly == NULL) return FALSE;
    1854     if (! n_Equal(r1->minpoly, r2->minpoly, r1->cf)) return FALSE;
    1855   }
    1856   else if (r2->minpoly != NULL) return FALSE;
     1846  if (r1->cf->minpoly != NULL)
     1847  {
     1848    if (r2->cf->minpoly == NULL) return FALSE;
     1849    if (! n_Equal(r1->cf->minpoly, r2->cf->minpoly, r1->cf)) return FALSE;
     1850  }
     1851  else if (r2->cf->minpoly != NULL) return FALSE;
    18571852
    18581853  if (qr)
    … …  
    22762271    }
    22772272  }
    2278   if (r->minpoly!=NULL)
    2279   {
    2280     omCheckAddr(r->minpoly);
     2273  if (r->cf->minpoly!=NULL)
     2274  {
     2275    omCheckAddr(r->cf->minpoly);
    22812276  }
    22822277  //assume(r->cf!=NULL);
    … …  
    31603155  {
    31613156    int i;
    3162     if (r->parameter!=NULL)
     3157    if (r->cf->parameter!=NULL)
    31633158    {
    31643159      for (i=0;i<rPar(r);i++)
    31653160      {
    3166         if(strlen(r->parameter[i])>1)
     3161        if(strlen(r->cf->parameter[i])>1)
    31673162        {
    31683163          r->ShortOut=FALSE;
    … …  
    55055500
    55065501        assume(MATELEM(src->GetNC()->C,i,j) != NULL);
    5507         MATELEM(C,nj,ni) = p_PermPoly(MATELEM(src->GetNC()->C,i,j),perm,src,r, nMap,par_perm,src->P);
     5502        MATELEM(C,nj,ni) = p_PermPoly(MATELEM(src->GetNC()->C,i,j),perm,src,r, nMap,par_perm,rPar(src));
    55085503
    55095504        if(MATELEM(src->GetNC()->D,i,j) != NULL)
    5510           MATELEM(D,nj,ni) = p_PermPoly(MATELEM(src->GetNC()->D,i,j),perm,src,r, nMap,par_perm,src->P);
     5505          MATELEM(D,nj,ni) = p_PermPoly(MATELEM(src->GetNC()->D,i,j),perm,src,r, nMap,par_perm,rPar(src));
    55115506      }
    55125507    }
Note: See TracChangeset for help on using the changeset viewer.