Changeset dd012ca in git for libpolys/polys/monomials


Ignore:
Timestamp:
Jan 17, 2011, 12:16:00 PM (13 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')
Children:
7eb7b589b58bd7d0f3394de9432ef2b9a9d03d57
Parents:
8e0bf8b4ab959ccb20087537bf84bdc77bdee721
git-author:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2011-01-17 12:16:00+01:00
git-committer:
Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 11:59:25+01:00
Message:
syntax fixes for libpoly, part 6: char_ptr removed
File:
1 edited

Legend:

Unmodified
Added
Removed
  • libpolys/polys/monomials/ring.cc

    r8e0bf8 rdd012ca  
    8686
    8787/*0 implementation*/
    88 //BOOLEAN rField_is_R(ring r=currRing)
     88//BOOLEAN rField_is_R(ring r)
    8989//{
    9090//  if (r->ch== -1)
     
    110110    rTest(r);
    111111    /*------------ set global ring vars --------------------------------*/
    112     currQuotient=r->qideal;
     112    //currQuotient=r->qideal;
    113113
    114114    /*------------ global variables related to coefficients ------------*/
     
    133133  /*r->P     = 0; Alloc0 */
    134134  /*names*/
    135   r->names = (char **) omAlloc0(N * sizeof(char_ptr));
     135  r->names = (char **) omAlloc0(N * sizeof(char *));
    136136  int i;
    137137  for(i=0;i<N;i++)
     
    253253  omCheckAddrSize(r->block1,nblocks*sizeof(int));
    254254  omCheckAddrSize(r->wvhdl,nblocks*sizeof(int *));
    255   omCheckAddrSize(r->names,r->N*sizeof(char_ptr));
     255  omCheckAddrSize(r->names,r->N*sizeof(char *));
    256256
    257257  nblocks--;
     
    260260  if (rField_is_GF(r))
    261261  {
    262     Print("//   # ground field : %d\n",rInternalChar(r));
     262    Print("//   # ground field : %d\n",r->ch);
    263263    Print("//   primitive element : %s\n", r->parameter[0]);
    264264    StringSetS("//   minpoly        : ");
     
    396396  {
    397397    PrintS("\n//   noncommutative relations:");
    398     if (r==currRing)
    399     {
    400       poly pl=NULL;
    401       int nl;
    402       int i,j;
    403       for (i = 1; i<r->N; i++)
    404       {
    405         for (j = i+1; j<=r->N; j++)
    406         {
    407           nl = nIsOne(p_GetCoeff(MATELEM(r->GetNC()->C,i,j),r->GetNC()->basering));
    408           if ( (MATELEM(r->GetNC()->D,i,j)!=NULL) || (!nl) )
    409           {
    410             Print("\n//    %s%s=",r->names[j-1],r->names[i-1]);
    411             pl = MATELEM(r->GetNC()->MT[UPMATELEM(i,j,r->N)],1,1);
    412             p_Write0(pl, r, r);
    413           }
    414         }
    415       }
    416     }
    417     else PrintS(" ...");
     398    poly pl=NULL;
     399    int nl;
     400    int i,j;
     401    for (i = 1; i<r->N; i++)
     402    {
     403      for (j = i+1; j<=r->N; j++)
     404      {
     405        nl = nIsOne(p_GetCoeff(MATELEM(r->GetNC()->C,i,j),r->GetNC()->basering));
     406        if ( (MATELEM(r->GetNC()->D,i,j)!=NULL) || (!nl) )
     407        {
     408          Print("\n//    %s%s=",r->names[j-1],r->names[i-1]);
     409          pl = MATELEM(r->GetNC()->MT[UPMATELEM(i,j,r->N)],1,1);
     410          p_Write0(pl, r, r);
     411        }
     412      }
     413    }
    418414#if 0  /*Singularg should not differ from Singular except in error case*/
    419415    Print("\n//   noncommutative type:%d", (int)ncRingType(r));
     
    489485      if (r->names[i] != NULL) omFree((ADDRESS)r->names[i]);
    490486    }
    491     omFreeSize((ADDRESS)r->names,r->N*sizeof(char_ptr));
     487    omFreeSize((ADDRESS)r->names,r->N*sizeof(char *));
    492488  }
    493489
     
    503499      j++;
    504500    }
    505     omFreeSize((ADDRESS)r->parameter,rPar(r)*sizeof(char_ptr));
     501    omFreeSize((ADDRESS)r->parameter,rPar(r)*sizeof(char *));
    506502  }
    507503#ifdef HAVE_RINGS
     
    796792int rSumInternal(ring r1, ring r2, ring &sum, BOOLEAN vartest, BOOLEAN dp_dp)
    797793{
    798   ring save=currRing;
    799794  ip_sring tmpR;
    800795  memset(&tmpR,0,sizeof(tmpR));
     
    810805        if (!vartest || (strcmp(r1->parameter[0],r2->parameter[0])==0)) /* 1 par */
    811806        {
    812           tmpR.parameter=(char **)omAllocBin(char_ptr_bin);
     807          tmpR.parameter=(char **)omAllocBin(char *);
    813808          tmpR.parameter[0]=omStrDup(r1->parameter[0]);
    814809          tmpR.P=1;
     
    832827              && n_Equal(r1->minpoly,r2->minpoly, r1))
    833828          {
    834             tmpR.parameter=(char **)omAllocBin(char_ptr_bin);
     829            tmpR.parameter=(char **)omAllocBin(char *);
    835830            tmpR.parameter[0]=omStrDup(r1->parameter[0]);
    836831            tmpR.minpoly=n_Copy(r1->minpoly, r1);
    837832            tmpR.P=1;
    838833            // HANNES: TODO: delete nSetChar
    839             rChangeCurrRing(save);
    840834          }
    841835          else
    842836          {
    843837            // HANNES: TODO: delete nSetChar
    844             rChangeCurrRing(save);
    845838            WerrorS("different minpolys");
    846839            return -1;
     
    852845              && (rPar(r2)==1))
    853846          {
    854             tmpR.parameter=(char **)omAllocBin(char_ptr_bin);
     847            tmpR.parameter=(char **)omAllocBin(char *);
    855848            tmpR.parameter[0]=omStrDup(r1->parameter[0]);
    856849            tmpR.P=1;
     
    871864              && (rPar(r1)==1))
    872865          {
    873             tmpR.parameter=(char **)omAllocBin(char_ptr_bin);
     866            tmpR.parameter=(char **)omAllocBin(char *);
    874867            tmpR.parameter[0]=omStrDup(r1->parameter[0]);
    875868            tmpR.P=1;
     
    885878        {
    886879          int len=rPar(r1)+rPar(r2);
    887           tmpR.parameter=(char **)omAlloc0(len*sizeof(char_ptr));
     880          tmpR.parameter=(char **)omAlloc0(len*sizeof(char *));
    888881          int i;
    889882          for (i=0;i<rPar(r1);i++)
     
    913906          {
    914907            tmpR.parameter=(char**)omReallocSize(tmpR.parameter,
    915                                                  len*sizeof(char_ptr),
    916                                                  i*sizeof(char_ptr));
     908                                                 len*sizeof(char *),
     909                                                 i*sizeof(char *));
    917910          }
    918911          tmpR.P=i;
     
    960953        tmpR.ch=rInternalChar(r1);
    961954        tmpR.P=rPar(r1);
    962         tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char_ptr));
     955        tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
    963956        int i;
    964957        for (i=0;i<rPar(r1);i++)
     
    988981        tmpR.ch=rInternalChar(r2);
    989982        tmpR.P=rPar(r2);
    990         tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char_ptr));
     983        tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
    991984        int i;
    992985        for (i=0;i<rPar(r2);i++)
     
    1004997        if (r2->parameter!=NULL)
    1005998        {
    1006           tmpR.parameter=(char **)omAllocBin(char_ptr_bin);
     999          tmpR.parameter=(char **)omAllocBin(char *);
    10071000          tmpR.P=1;
    10081001          tmpR.parameter[0]=omStrDup(r2->parameter[0]);
     
    10211014        tmpR.ch=rInternalChar(r1);
    10221015        tmpR.P=rPar(r1);
    1023         tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char_ptr));
     1016        tmpR.parameter=(char **)omAlloc(rPar(r1)*sizeof(char *));
    10241017        int i;
    10251018        for(i=0;i<rPar(r1);i++)
     
    10481041        tmpR.ch=rInternalChar(r2);
    10491042        tmpR.P=rPar(r2);
    1050         tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char_ptr));
     1043        tmpR.parameter=(char **)omAlloc(rPar(r2)*sizeof(char *));
    10511044        int i;
    10521045        for(i=0;i<rPar(r2);i++)
     
    10691062  int i,j,k;
    10701063  int l=r1->N+r2->N;
    1071   char **names=(char **)omAlloc0(l*sizeof(char_ptr));
     1064  char **names=(char **)omAlloc0(l*sizeof(char *));
    10721065  k=0;
    10731066
     
    13191312    {
    13201313      for(i=0;i<k;i++) omFree((ADDRESS)tmpR.names[i]);
    1321       omFreeSize((ADDRESS)names,tmpR.N*sizeof(char_ptr));
     1314      omFreeSize((ADDRESS)names,tmpR.N*sizeof(char *));
    13221315      Werror("difficulties with variables: %d,%d -> %d",rVar(r1),rVar(r2),k);
    13231316      return -1;
     
    14961489  if (r2->qideal!=NULL)
    14971490  {
    1498     if (currRing!=sum)
    1499       rChangeCurrRing(sum);
     1491    //if (currRing!=sum)
     1492    //  rChangeCurrRing(sum);
    15001493    int *perm2 = (int *)omAlloc0((rVar(r2)+1)*sizeof(int));
    15011494    int *par_perm2 = NULL;
     
    15071500    Q2 = idInit(IDELEMS(r2->qideal),1);
    15081501    for (int for_i=0;for_i<IDELEMS(r2->qideal);for_i++)
    1509       Q2->m[for_i] = pPermPoly(r2->qideal->m[for_i],perm2,r2,nMap2,par_perm2,rPar(r2));
     1502      Q2->m[for_i] = p_PermPoly(r2->qideal->m[for_i],perm2,r2,nMap2,par_perm2,rPar(r2),sum);
    15101503    omFree((ADDRESS)perm2);
    15111504  }
     
    15131506  {
    15141507    Q = idSimpleAdd(Q1,Q2);
    1515     rChangeCurrRing(old_ring2);
     1508    //rChangeCurrRing(old_ring2);
    15161509  }
    15171510  sum->qideal = Q;
     
    16441637  if (r->parameter!=NULL)
    16451638  {
    1646     if (r->minpoly!=NULL) res->minpoly=nCopy(r->minpoly);
     1639    if (r->minpoly!=NULL) res->minpoly=n_Copy(r->minpoly,r->cf);
    16471640    int l=rPar(r);
    1648     res->parameter=(char **)omAlloc(l*sizeof(char_ptr));
     1641    res->parameter=(char **)omAlloc(l*sizeof(char *));
    16491642    int i;
    16501643    for(i=0;i<rPar(r);i++)
     
    16851678  //memset: }
    16861679
    1687   res->names   = (char **)omAlloc0(rVar(r) * sizeof(char_ptr));
     1680  res->names   = (char **)omAlloc0(rVar(r) * sizeof(char *));
    16881681  for (i=0; i<rVar(res); i++)
    16891682  {
     
    17471740  if (r1 == NULL || r2 == NULL) return FALSE;
    17481741
    1749   if ((rInternalChar(r1) != rInternalChar(r2))
     1742  if ((r1->cf->type != r2->cf->type)
    17501743  || (r1->float_len != r2->float_len)
    17511744  || (r1->float_len2 != r2->float_len2)
     
    18161809        m2 = id2->m;
    18171810        for (i=0; i<n; i++)
    1818           if (! pEqualPolys(m1[i],m2[i])) return FALSE;
     1811          if (! p_EqualPolys(m1[i],m2[i],r1)) return FALSE;
    18191812      }
    18201813    }
     
    18361829  if (r1 == NULL || r2 == NULL) return FALSE;
    18371830
    1838   if ((rInternalChar(r1) != rInternalChar(r2))
     1831  if ((r1->cf->type != r2->cf->type)
    18391832  || (r1->float_len != r2->float_len)
    18401833  || (r1->float_len2 != r2->float_len2)
     
    51595152n_coeffType rFieldType(ring r)
    51605153{
     5154  return (r->cf->type);
    51615155  if (rField_is_Zp(r))     return n_Zp;
    51625156  if (rField_is_Q(r))      return n_Q;
Note: See TracChangeset for help on using the changeset viewer.