Changeset 896561 in git


Ignore:
Timestamp:
Jul 4, 2011, 4:48:27 PM (12 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
Children:
0312c51932007766f7474511cc741c2233b62ef9
Parents:
5f36297071b90604201cce58a71b12514e338582
Message:
experiment with threads: p_Mult_nn

git-svn-id: file:///usr/local/Singular/svn/trunk@14320 2c84dea3-7e68-4137-9b89-c4e89433aadc
Files:
1 added
17 edited

Legend:

Unmodified
Added
Removed

  • Singular/interpolation.cc

    r5f3629 r896561  
    17541754          if (mpz_sgn(temp->polycoef[a])!=0)
    17551755          {
    1756              number n=(number)omAllocBin(rnumber_bin);
     1756             number n=ALLOC_RNUMBER();
    17571757#ifdef LDEBUG
    17581758             n->debug=123456;

  • Singular/misc_ip.cc

    r5f3629 r896561  
    3232number mpz2number(mpz_t m)
    3333{
    34   number z = (number)omAllocBin(rnumber_bin);
     34  number z = ALLOC_RNUMBER();
    3535  mpz_init_set(z->z, m);
    3636  mpz_init_set_ui(z->n, 1);

  • Singular/mpsr_Get.cc

    r5f3629 r896561  
    172172static inline void InitApIntLeftv(mpsr_leftv mlv, mpz_ptr apint)
    173173{
    174   number n = (number) omAllocBin(rnumber_bin);
     174  number n = (number) ALLOC_RNUMBER();
    175175#if defined(LDEBUG)
    176176    n->debug=123456;

  • Singular/mpsr_GetPoly.cc

    r5f3629 r896561  
    230230  {
    231231    mpz_ptr gnum;
    232     y =  (number) omAlloc0Bin(rnumber_bin);
     232    y =  ALLOC0_RNUMBER();
    233233#if defined(LDEBUG)
    234234    y->debug = 123456;
    … …  
    251251      if (req) return mpsr_SetError(mpsr_ReqAnnotSkip);
    252252    }
    253     *x =  (number) omAlloc0Bin(rnumber_bin);
     253    *x =  (number) ALLOC0_RNUMBER();
    254254    y = (number) *x;
    255255#if defined(LDEBUG)
    … …  
    276276    {
    277277      // otherwise, make an apint out of it
    278       *x =  (number) omAlloc0Bin(rnumber_bin);
     278      *x =  (number) ALLOC0_RNUMBER();
    279279      y = (number) *x;
    280280#if defined(LDEBUG)
    … …  
    353353  {
    354354    // single number
    355     b = (lnumber) omAlloc0Bin(rnumber_bin);
     355    b = (lnumber) ALLOC0_RNUMBER();
    356356    *a = (number) b;
    357357    failr(GetAlgPoly(link, &(b->z)));

  • Singular/tesths.cc

    r5f3629 r896561  
    4747
    4848#if ! defined(LIBSINGULAR)
    49 int mmInit( void )
     49int mmInit2( void )
    5050{
    5151#if defined(OMALLOC_USES_MALLOC) || defined(X_OMALLOC)
    … …  
    5959  return 1;
    6060}
     61int mmInit( void )
     62{
     63#ifndef SI_THREADS
     64  return mmInit2();
     65#else
     66  return 1;
     67#endif
     68}
    6169
    6270/*0 implementation*/
    … …  
    6573    char** argv)   /* parameter array */
    6674{
    67   mmInit();
     75  //mmInit();
    6876  // Don't worry: ifdef OM_NDEBUG, then all these calls are undef'ed
    6977  omInitRet_2_Info(argv[0]);

  • kernel/clapconv.cc

    r5f3629 r896561  
    8282  else
    8383  {
    84     number z=(number)omAllocBin(rnumber_bin);
     84    number z=ALLOC_RNUMBER();
    8585#if defined(LDEBUG)
    8686    z->debug=123456;
    … …  
    144144    else
    145145    {
    146       number z=(number)omAllocBin(rnumber_bin);
     146      number z=ALLOC_RNUMBER();
    147147#if defined(LDEBUG)
    148148      z->debug=123456;
    … …  
    300300          p_AddExp(z,i,exp[i-par_start],r->algring);
    301301      }
    302       pGetCoeff(term)=(number)omAlloc0Bin(rnumber_bin);
     302      pGetCoeff(term)=(number)ALLOC0_LNUMBER();
    303303      ((lnumber)pGetCoeff(term))->z=z;
    304304      p_Setm( term,r );
    … …  
    358358  else
    359359  {
    360     number z=(number)omAllocBin(rnumber_bin);
     360    number z=ALLOC_RNUMBER();
    361361#if defined(LDEBUG)
    362362    z->debug=123456;
    … …  
    471471      p_SetExp( term, i ,exp[i], r);
    472472    //if (rRing_has_Comp(currRing)) p_SetComp(term, 0, currRing); // done by pInit
    473     pGetCoeff(term)=(number)omAlloc0Bin(rnumber_bin);
     473    pGetCoeff(term)=(number)ALLOC0_LNUMBER();
    474474    ((lnumber)pGetCoeff(term))->z=convFactoryPSingP( f, r->algring );
    475475    p_Setm( term,r );

  • kernel/ideals.cc

    r5f3629 r896561  
    41574157  if (rField_has_simple_inverse()) return; /* Z/p, GF(p,n), R, long R/C */
    41584158  int i;
    4159   poly p;
    41604159  for(i=IDELEMS(I)-1;i>=0;i--)
    41614160  {
    4162     p=I->m[i] ;
    4163     while(p!=NULL)
    4164     {
    4165       nNormalize(pGetCoeff(p));
    4166       pIter(p);
    4167     }
     4161    pNormalize(I->m[i]);
    41684162  }
    41694163}

  • kernel/longalg.cc

    r5f3629 r896561  
    3939omBin   snaIdeal_bin = omGetSpecBin(sizeof(snaIdeal));
    4040number  (*naMap)(number from);
     41//omBin lnumber_bin = omGetSpecBin(sizeof(slnumber));
     42//omBin rnumber_bin = omGetSpecBin(sizeof(snumber));
    4143
    4244void redefineFunctionPointers()
    … …  
    200202      poly z=p_Init(r->algring);
    201203      pSetCoeff0(z,c);
    202       lnumber l = (lnumber)omAllocBin(rnumber_bin);
     204      lnumber l = ALLOC_LNUMBER();
    203205      l->z = z;
    204206      l->s = 2;
    … …  
    213215number  naPar(int i)
    214216{
    215   lnumber l = (lnumber)omAllocBin(rnumber_bin);
     217  lnumber l = ALLOC_LNUMBER();
    216218  l->s = 2;
    217219  l->z = p_ISet(1,nacRing);
    … …  
    270272    p_Delete(&(l->z),r->algring);
    271273    p_Delete(&(l->n),r->algring);
    272     omFreeBin((ADDRESS)l,  rnumber_bin);
     274    FREE_LNUMBER(l);
    273275  }
    274276  *p = NULL;
    … …  
    284286  lnumber erg;
    285287  lnumber src = (lnumber)p;
    286   erg = (lnumber)omAlloc0Bin(rnumber_bin);
     288  erg = ALLOC_LNUMBER();
    287289  erg->z = p_Copy(src->z, nacRing);
    288290  erg->n = p_Copy(src->n, nacRing);
    … …  
    295297  lnumber erg;
    296298  lnumber src = (lnumber)p;
    297   erg = (lnumber)omAlloc0Bin(rnumber_bin);
     299  erg = (lnumber)ALLOC_LNUMBER();
    298300  erg->z = p_Copy(src->z,r->algring);
    299301  erg->n = p_Copy(src->n,r->algring);
    … …  
    315317  lnumber b = (lnumber)lb;
    316318  #ifdef LDEBUG
    317   omCheckAddrSize(a,sizeof(snumber));
    318   omCheckAddrSize(b,sizeof(snumber));
     319  omCheckAddrSize(a,sizeof(slnumber));
     320  omCheckAddrSize(b,sizeof(slnumber));
    319321  #endif
    320322  if (b->n!=NULL) x = pp_Mult_qq(a->z, b->n,nacRing);
    … …  
    327329    return (number)NULL;
    328330  }
    329   lu = (lnumber)omAllocBin(rnumber_bin);
     331  lu = ALLOC_LNUMBER();
    330332  lu->z=res;
    331333  if (a->n!=NULL)
    … …  
    383385
    384386  #ifdef LDEBUG
    385   omCheckAddrSize(a,sizeof(snumber));
    386   omCheckAddrSize(b,sizeof(snumber));
     387  omCheckAddrSize(a,sizeof(slnumber));
     388  omCheckAddrSize(b,sizeof(slnumber));
    387389  #endif
    388390
    … …  
    397399    return (number)NULL;
    398400  }
    399   lu = (lnumber)omAllocBin(rnumber_bin);
     401  lu = ALLOC_LNUMBER();
    400402  lu->z=res;
    401403  if (a->n!=NULL)
    … …  
    439441
    440442  #ifdef LDEBUG
    441   omCheckAddrSize(a,sizeof(snumber));
    442   omCheckAddrSize(b,sizeof(snumber));
     443  omCheckAddrSize(a,sizeof(slnumber));
     444  omCheckAddrSize(b,sizeof(slnumber));
    443445  #endif
    444446  naTest(la);
    445447  naTest(lb);
    446448
    447   lo = (lnumber)omAllocBin(rnumber_bin);
     449  lo = ALLOC_LNUMBER();
    448450  lo->z = pp_Mult_qq(a->z, b->z,nacRing);
    449451
    … …  
    493495  if(lo->z==NULL)
    494496  {
    495     omFreeBin((ADDRESS)lo, rnumber_bin);
     497    FREE_LNUMBER(lo);
    496498    lo=NULL;
    497499  }
    … …  
    525527  assume(a->z!=NULL && b->z!=NULL);
    526528  assume(a->n==NULL && b->n==NULL);
    527   res = (lnumber)omAllocBin(rnumber_bin);
     529  res = ALLOC_LNUMBER();
    528530  res->z = napCopy(a->z);
    529531  res->n = napCopy(b->z);
    … …  
    556558  }
    557559  #ifdef LDEBUG
    558   omCheckAddrSize(a,sizeof(snumber));
    559   omCheckAddrSize(b,sizeof(snumber));
     560  omCheckAddrSize(a,sizeof(slnumber));
     561  omCheckAddrSize(b,sizeof(slnumber));
    560562  #endif
    561   lo = (lnumber)omAllocBin(rnumber_bin);
     563  lo = ALLOC_LNUMBER();
    562564  if (b->n!=NULL)
    563565    lo->z = pp_Mult_qq(a->z, b->n,nacRing);
    … …  
    633635  }
    634636  #ifdef LDEBUG
    635   omCheckAddrSize(b,sizeof(snumber));
     637  omCheckAddrSize(b,sizeof(slnumber));
    636638  #endif
    637   lo = (lnumber)omAlloc0Bin(rnumber_bin);
     639  lo = ALLOC0_LNUMBER();
    638640  lo->s = b->s;
    639641  if (b->n!=NULL)
    … …  
    784786    return s;
    785787  }
    786   *p = (number)omAlloc0Bin(rnumber_bin);
     788  *p = (number)ALLOC0_LNUMBER();
    787789  a = (lnumber)*p;
    788790  if ((naMinimalPoly!=NULL)
    … …  
    799801  if(a->z==NULL)
    800802  {
    801     omFreeBin((ADDRESS)*p, rnumber_bin);
     803    FREE_LNUMBER(a);
    802804    *p=NULL;
    803805  }
    … …  
    879881  if (a==NULL) return FALSE;
    880882#ifdef LDEBUG
    881   omCheckAddrSize(a,sizeof(snumber));
     883  omCheckAddrSize(a,sizeof(slnumber));
    882884  if (a->z==NULL)
    883885  {
    … …  
    937939  if (a==NULL) return FALSE;
    938940#ifdef LDEBUG
    939   omCheckAddrSize(a,sizeof(snumber));
     941  omCheckAddrSize(a,sizeof(slnumber));
    940942  if (a->z==NULL)
    941943  {
    … …  
    976978
    977979  lnumber x, y;
    978   lnumber result = (lnumber)omAlloc0Bin(rnumber_bin);
     980  lnumber result = ALLOC0_LNUMBER();
    979981
    980982  x = (lnumber)a;
    … …  
    13761378  lnumber a = (lnumber)la;
    13771379  lnumber b = (lnumber)lb;
    1378   result = (lnumber)omAlloc0Bin(rnumber_bin);
     1380  result = ALLOC0_LNUMBER();
    13791381  naTest(la);
    13801382  naTest(lb);
    … …  
    14621464{
    14631465  if (npIsZero(c)) return NULL;
    1464   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     1466  lnumber l=ALLOC_LNUMBER();
    14651467  l->s=2;
    14661468  l->z=(napoly)p_Init(nacRing);
    … …  
    14781480{
    14791481  if (nlIsZero(c)) return NULL;
    1480   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     1482  lnumber l=ALLOC_LNUMBER();
    14811483  l->s=0;
    14821484  l->z=(napoly)p_Init(nacRing);
    … …  
    14921494{
    14931495  if (npIsZero(c)) return NULL;
    1494   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     1496  lnumber l=ALLOC_LNUMBER();
    14951497  l->s=2;
    14961498  l->z=(napoly)p_Init(nacRing);
    … …  
    15101512  number n=npInit(i,ntMapRing);
    15111513  if (npIsZero(n)) return NULL;
    1512   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     1514  lnumber l=ALLOC_LNUMBER();
    15131515  l->s=2;
    15141516  l->z=(napoly)p_Init(nacRing);
    … …  
    15271529  if (npIsZero(n)) return NULL;
    15281530  npTest(n);
    1529   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     1531  lnumber l=ALLOC_LNUMBER();
    15301532  l->s=2;
    15311533  l->z=(napoly)p_Init(nacRing);
    … …  
    15411543{
    15421544  if (c==NULL) return NULL;
    1543   lnumber erg= (lnumber)omAlloc0Bin(rnumber_bin);
     1545  lnumber erg= ALLOC_LNUMBER();
    15441546  lnumber src =(lnumber)c;
    15451547  erg->s=src->s;
    … …  
    16521654    return TRUE;
    16531655  #ifdef LDEBUG
    1654   omCheckAddrSize(a, sizeof(snumber));
     1656  omCheckAddrSize(a, sizeof(slnumber));
    16551657  #endif
    16561658  napoly p = x->z;

  • kernel/longalg.h

    r5f3629 r896561  
    100100extern number (*naMap)(number from);
    101101
     102//extern omBin lnumber_bin;
     103//#define ALLOC_LNUMBER() (lnumber)omAllocBin(lnumber_bin)
     104//#define ALLOC0_LNUMBER() (lnumber)omAlloc0Bin(lnumber_bin)
     105//#define FREE_LNUMBER(x) omFreeBin((ADDRESS)x, lnumber_bin)
     106extern omBin rnumber_bin;
     107#define ALLOC_LNUMBER() (lnumber)omAllocBin(rnumber_bin)
     108#define ALLOC0_LNUMBER() (lnumber)omAlloc0Bin(rnumber_bin)
     109#define FREE_LNUMBER(x) omFreeBin((ADDRESS)x, rnumber_bin)
    102110#endif
    103111

  • kernel/longrat.cc

    r5f3629 r896561  
    3333    {
    3434      mpz_clear(x->z);
    35       omFreeBin((ADDRESS)x, rnumber_bin);
     35      FREE_RNUMBER(x);
    3636      return INT_TO_SR(ui);
    3737    }
    … …  
    120120}
    121121
     122//#ifndef SI_THREADS
    122123omBin rnumber_bin = omGetSpecBin(sizeof(snumber));
     124//#endif
    123125
    124126number nlOne=INT_TO_SR(1);
    … …  
    154156number nlMapGMP(number from)
    155157{
    156   number z=(number)omAllocBin(rnumber_bin);
     158  number z=ALLOC_RNUMBER();
    157159#if defined(LDEBUG)
    158160  z->debug=123456;
    … …  
    170172number nlMapMachineInt(number from)
    171173{
    172   number z=(number)omAllocBin(rnumber_bin);
     174  number z=ALLOC_RNUMBER();
    173175#if defined(LDEBUG)
    174176  z->debug=123456;
    … …  
    248250    return FALSE;
    249251  }
    250   omCheckAddrSize(a->z[0]._mp_d,a->z[0]._mp_alloc*BYTES_PER_MP_LIMB);
     252  // only, if gmp memeory routine are from omalloc:
     253  //omCheckAddrSize(a->z[0]._mp_d,a->z[0]._mp_alloc*BYTES_PER_MP_LIMB);
    251254  if (a->z[0]._mp_alloc==0)
    252255    Print("!!longrat:z->alloc=0 in %s:%d\n",f,l);
    … …  
    254257  if (a->s<2)
    255258  {
    256     omCheckIf(omCheckAddrSize(a->n[0]._mp_d,a->n[0]._mp_alloc*BYTES_PER_MP_LIMB), return FALSE);
     259    // only, if gmp memeory routine are from omalloc:
     260    //omCheckIf(omCheckAddrSize(a->n[0]._mp_d,a->n[0]._mp_alloc*BYTES_PER_MP_LIMB), return FALSE);
    257261    if (a->z[0]._mp_alloc==0)
    258262      Print("!!longrat:n->alloc=0 in %s:%d\n",f,l);
    … …  
    349353
    350354  e=(*f)[0]._mp_exp-size;
    351   res = (number)omAllocBin(rnumber_bin);
     355  res = ALLOC_RNUMBER();
    352356#if defined(LDEBUG)
    353357  res->debug=123456;
    … …  
    551555      return INT_TO_SR(0);
    552556    }
    553     n=(number)omAllocBin(rnumber_bin);
     557    n=ALLOC_RNUMBER();
    554558#if defined(LDEBUG)
    555559    n->debug=123456;
    … …  
    569573    return n;
    570574  }
    571   n=(number)omAllocBin(rnumber_bin);
     575  n=ALLOC_RNUMBER();
    572576#if defined(LDEBUG)
    573577  n->debug=123456;
    … …  
    642646    b=bb;
    643647  }
    644   u=(number)omAllocBin(rnumber_bin);
     648  u=ALLOC_RNUMBER();
    645649#if defined(LDEBUG)
    646650  u->debug=123456;
    … …  
    656660    bb->debug=654324;
    657661#endif
    658     omFreeBin((ADDRESS)bb, rnumber_bin);
     662    FREE_RNUMBER(bb);
    659663  }
    660664  u=nlShort3(u);
    … …  
    706710    b=bb;
    707711  }
    708   u=(number)omAllocBin(rnumber_bin);
     712  u=ALLOC_RNUMBER();
    709713#if defined(LDEBUG)
    710714  u->debug=123456;
    … …  
    722726    bb->debug=654324;
    723727#endif
    724     omFreeBin((ADDRESS)bb, rnumber_bin);
     728    FREE_RNUMBER(bb);
    725729  }
    726730  u=nlShort3(u);
    … …  
    786790    b=bb;
    787791  }
    788   u=(number)omAllocBin(rnumber_bin);
     792  u=ALLOC_RNUMBER();
    789793#if defined(LDEBUG)
    790794  u->debug=123456;
    … …  
    799803    bb->debug=654324;
    800804#endif
    801     omFreeBin((ADDRESS)bb, rnumber_bin);
     805    FREE_RNUMBER(bb);
    802806  }
    803807  if (mpz_isNeg(u->z))
    … …  
    824828    return INT_TO_SR(0);
    825829  }
    826   u=(number)omAllocBin(rnumber_bin);
     830  u=ALLOC_RNUMBER();
    827831  u->s=0;
    828832#if defined(LDEBUG)
    … …  
    836840    if ((i==-POW_2_28) && (j== -1L))
    837841    { 
    838       omFreeBin((ADDRESS)u, rnumber_bin);
     842      FREE_RNUMBER(u);
    839843      return nlRInit(POW_2_28);
    840844    }
    … …  
    842846    if (r==0)
    843847    {
    844       omFreeBin((ADDRESS)u, rnumber_bin);
     848      FREE_RNUMBER(u);
    845849      return INT_TO_SR(i/j);
    846850    }
    … …  
    867871      {
    868872        mpz_clear(u->z);
    869         omFreeBin((ADDRESS)u, rnumber_bin);
     873        FREE_RNUMBER(u);
    870874        return INT_TO_SR(1);
    871875      }
    … …  
    936940    else if (x->s==0)
    937941      nlNormalize(x);
    938     *u=(number)omAllocBin(rnumber_bin);
     942    *u=ALLOC_RNUMBER();
    939943#if defined(LDEBUG)
    940944    (*u)->debug=123456;
    … …  
    960964    {
    961965      mpz_clear(aa->z);
    962       omFreeBin((ADDRESS)aa, rnumber_bin);
     966      FREE_RNUMBER(aa);
    963967    }
    964968  }
    … …  
    10751079  else
    10761080  {
    1077     result=(number)omAllocBin(rnumber_bin);
     1081    result=ALLOC_RNUMBER();
    10781082    mpz_init(result->z);
    10791083    mpz_gcd(result->z,a->z,b->z);
    … …  
    11691173    return nlCopy(a);
    11701174  }
    1171   result=(number)omAllocBin(rnumber_bin);
     1175  result=ALLOC_RNUMBER();
    11721176#if defined(LDEBUG)
    11731177  result->debug=123456;
    … …  
    12611265      if (n->s!=3)
    12621266      {
    1263         number u=(number)omAllocBin(rnumber_bin);
     1267        number u=ALLOC_RNUMBER();
    12641268        u->s=3;
    12651269#if defined(LDEBUG)
    … …  
    12881292    if (!(SR_HDL(n) & SR_INT))
    12891293    {
    1290       number u=(number)omAllocBin(rnumber_bin);
     1294      number u=ALLOC_RNUMBER();
    12911295#if defined(LDEBUG)
    12921296      u->debug=123456;
    … …  
    13601364  assume(!((SR_HDL(a) & SR_INT)||(a==NULL)));
    13611365  nlTest(a);
    1362   number b=(number)omAllocBin(rnumber_bin);
     1366  number b=ALLOC_RNUMBER();
    13631367#if defined(LDEBUG)
    13641368  b->debug=123456;
    … …  
    13921396#endif
    13931397    }
    1394     omFreeBin((ADDRESS) *a, rnumber_bin);
     1398    FREE_RNUMBER(*a);
    13951399  }
    13961400}
    … …  
    14111415number _nlAdd_aNoImm_OR_bNoImm(number a, number b)
    14121416{
    1413   number u=(number)omAllocBin(rnumber_bin);
     1417  number u=ALLOC_RNUMBER();
    14141418#if defined(LDEBUG)
    14151419  u->debug=123456;
    … …  
    14371441        {
    14381442          mpz_clear(u->z);
    1439           omFreeBin((ADDRESS)u, rnumber_bin);
     1443          FREE_RNUMBER(u);
    14401444          return INT_TO_SR(0);
    14411445        }
    … …  
    14431447        {
    14441448          mpz_clear(u->z);
    1445           omFreeBin((ADDRESS)u, rnumber_bin);
     1449          FREE_RNUMBER(u);
    14461450          return INT_TO_SR(1);
    14471451        }
    … …  
    14591463        {
    14601464          mpz_clear(u->z);
    1461           omFreeBin((ADDRESS)u, rnumber_bin);
     1465          FREE_RNUMBER(u);
    14621466          return INT_TO_SR(0);
    14631467        }
    … …  
    14911495            {
    14921496              mpz_clear(u->z);
    1493               omFreeBin((ADDRESS)u, rnumber_bin);
     1497              FREE_RNUMBER(u);
    14941498              return INT_TO_SR(0);
    14951499            }
    … …  
    15001504               mpz_clear(u->z);
    15011505               mpz_clear(u->n);
    1502                omFreeBin((ADDRESS)u, rnumber_bin);
     1506               FREE_RNUMBER(u);
    15031507               return INT_TO_SR(1);
    15041508            }
    … …  
    15131517            {
    15141518              mpz_clear(u->z);
    1515               omFreeBin((ADDRESS)u, rnumber_bin);
     1519              FREE_RNUMBER(u);
    15161520              return INT_TO_SR(0);
    15171521            }
    … …  
    15191523            {
    15201524              mpz_clear(u->z);
    1521               omFreeBin((ADDRESS)u, rnumber_bin);
     1525              FREE_RNUMBER(u);
    15221526              return INT_TO_SR(1);
    15231527            }
    … …  
    15411545            {
    15421546              mpz_clear(u->z);
    1543               omFreeBin((ADDRESS)u, rnumber_bin);
     1547              FREE_RNUMBER(u);
    15441548              return INT_TO_SR(0);
    15451549            }
    … …  
    15471551            {
    15481552              mpz_clear(u->z);
    1549               omFreeBin((ADDRESS)u, rnumber_bin);
     1553              FREE_RNUMBER(u);
    15501554              return INT_TO_SR(1);
    15511555            }
    … …  
    15601564            {
    15611565              mpz_clear(u->z);
    1562               omFreeBin((ADDRESS)u, rnumber_bin);
     1566              FREE_RNUMBER(u);
    15631567              return INT_TO_SR(0);
    15641568            }
    … …  
    15781582number _nlSub_aNoImm_OR_bNoImm(number a, number b)
    15791583{
    1580   number u=(number)omAllocBin(rnumber_bin);
     1584  number u=ALLOC_RNUMBER();
    15811585#if defined(LDEBUG)
    15821586  u->debug=123456;
    … …  
    15981602        {
    15991603          mpz_clear(u->z);
    1600           omFreeBin((ADDRESS)u, rnumber_bin);
     1604          FREE_RNUMBER(u);
    16011605          return INT_TO_SR(0);
    16021606        }
    … …  
    16041608        {
    16051609          mpz_clear(u->z);
    1606           omFreeBin((ADDRESS)u, rnumber_bin);
     1610          FREE_RNUMBER(u);
    16071611          return INT_TO_SR(1);
    16081612        }
    … …  
    16261630        {
    16271631          mpz_clear(u->z);
    1628           omFreeBin((ADDRESS)u, rnumber_bin);
     1632          FREE_RNUMBER(u);
    16291633          return INT_TO_SR(0);
    16301634        }
    … …  
    16501654        {
    16511655          mpz_clear(u->z);
    1652           omFreeBin((ADDRESS)u, rnumber_bin);
     1656          FREE_RNUMBER(u);
    16531657          return INT_TO_SR(0);
    16541658        }
    … …  
    16561660        {
    16571661          mpz_clear(u->z);
    1658           omFreeBin((ADDRESS)u, rnumber_bin);
     1662          FREE_RNUMBER(u);
    16591663          return INT_TO_SR(1);
    16601664        }
    … …  
    16761680        {
    16771681          mpz_clear(u->z);
    1678           omFreeBin((ADDRESS)u, rnumber_bin);
     1682          FREE_RNUMBER(u);
    16791683          return INT_TO_SR(0);
    16801684        }
    … …  
    17091713            {
    17101714              mpz_clear(u->z);
    1711               omFreeBin((ADDRESS)u, rnumber_bin);
     1715              FREE_RNUMBER(u);
    17121716              return INT_TO_SR(0);
    17131717            }
    … …  
    17181722              mpz_clear(u->z);
    17191723              mpz_clear(u->n);
    1720               omFreeBin((ADDRESS)u, rnumber_bin);
     1724              FREE_RNUMBER(u);
    17211725              return INT_TO_SR(1);
    17221726            }
    … …  
    17341738            {
    17351739              mpz_clear(u->z);
    1736               omFreeBin((ADDRESS)u, rnumber_bin);
     1740              FREE_RNUMBER(u);
    17371741              return INT_TO_SR(0);
    17381742            }
    … …  
    17401744            {
    17411745              mpz_clear(u->z);
    1742               omFreeBin((ADDRESS)u, rnumber_bin);
     1746              FREE_RNUMBER(u);
    17431747              return INT_TO_SR(1);
    17441748            }
    … …  
    17651769            {
    17661770              mpz_clear(u->z);
    1767               omFreeBin((ADDRESS)u, rnumber_bin);
     1771              FREE_RNUMBER(u);
    17681772              return INT_TO_SR(0);
    17691773            }
    … …  
    17711775            {
    17721776              mpz_clear(u->z);
    1773               omFreeBin((ADDRESS)u, rnumber_bin);
     1777              FREE_RNUMBER(u);
    17741778              return INT_TO_SR(1);
    17751779            }
    … …  
    17841788            {
    17851789              mpz_clear(u->z);
    1786               omFreeBin((ADDRESS)u, rnumber_bin);
     1790              FREE_RNUMBER(u);
    17871791              return INT_TO_SR(0);
    17881792            }
    … …  
    18031807number _nlMult_aImm_bImm_rNoImm(number a, number b)
    18041808{
    1805   number u=(number)omAllocBin(rnumber_bin);
     1809  number u=ALLOC_RNUMBER();
    18061810#if defined(LDEBUG)
    18071811  u->debug=123456;
    … …  
    18181822{
    18191823  assume(! (SR_HDL(a) & SR_HDL(b) & SR_INT));
    1820   number u=(number)omAllocBin(rnumber_bin);
     1824  number u=ALLOC_RNUMBER();
    18211825#if defined(LDEBUG)
    18221826  u->debug=123456;
    … …  
    18561860      {
    18571861        mpz_clear(u->z);
    1858         omFreeBin((ADDRESS)u, rnumber_bin);
     1862        FREE_RNUMBER(u);
    18591863        return INT_TO_SR(1);
    18601864      }
    … …  
    18811885        {
    18821886          mpz_clear(u->z);
    1883           omFreeBin((ADDRESS)u, rnumber_bin);
     1887          FREE_RNUMBER(u);
    18841888          return INT_TO_SR(1);
    18851889        }
    … …  
    18941898        {
    18951899          mpz_clear(u->z);
    1896           omFreeBin((ADDRESS)u, rnumber_bin);
     1900          FREE_RNUMBER(u);
    18971901          return INT_TO_SR(1);
    18981902        }
    … …  
    19071911          mpz_clear(u->z);
    19081912          mpz_clear(u->n);
    1909           omFreeBin((ADDRESS)u, rnumber_bin);
     1913          FREE_RNUMBER(u);
    19101914          return INT_TO_SR(1);
    19111915        }
    … …  
    19221926number nlRInit (long i)
    19231927{
    1924   number z=(number)omAllocBin(rnumber_bin);
     1928  number z=ALLOC_RNUMBER();
    19251929#if defined(LDEBUG)
    19261930  z->debug=123456;
    … …  
    19361940number nlInit2 (int i, int j)
    19371941{
    1938   number z=(number)omAllocBin(rnumber_bin);
     1942  number z=ALLOC_RNUMBER();
    19391943#if defined(LDEBUG)
    19401944  z->debug=123456;
    … …  
    19491953number nlInit2gmp (mpz_t i, mpz_t j)
    19501954{
    1951   number z=(number)omAllocBin(rnumber_bin);
     1955  number z=ALLOC_RNUMBER();
    19521956#if defined(LDEBUG)
    19531957  z->debug=123456;
    … …  
    21322136  if (SR_HDL(a) & SR_INT)
    21332137  {
    2134     number u=(number)omAllocBin(rnumber_bin);
     2138    number u=ALLOC_RNUMBER();
    21352139    #if defined(LDEBUG)
    21362140    u->debug=123456;
    … …  
    22832287LINLINE void nlInpMult(number &a, number b, const ring r)
    22842288{
    2285   if (((SR_HDL(b)|SR_HDL(a))&SR_INT))
    2286   {
    2287     number n=nlMult(a,b);
     2289  number aa=a;
     2290  if (((SR_HDL(b)|SR_HDL(aa))&SR_INT))
     2291  {
     2292    number n=nlMult(aa,b);
    22882293    nlDelete(&a,r);
    22892294    a=n;
    … …  
    22912296  else
    22922297  {
    2293     mpz_mul(a->z,a->z,b->z);
    2294     if (a->s==3)
     2298    mpz_mul(aa->z,a->z,b->z);
     2299    if (aa->s==3)
    22952300    {
    22962301      if(b->s!=3)
    … …  
    23902395    {
    23912396       // return N/B
    2392        z=(number)omAllocBin(rnumber_bin);
     2397       z=ALLOC_RNUMBER();
    23932398       #ifdef LDEBUG
    23942399       z->debug=123456;
    … …  
    24572462    {
    24582463      mpz_clear(&r->z);
    2459       omFreeBin((ADDRESS)r, rnumber_bin);
     2464      FREE_RNUMBER(r);
    24602465      r=INT_TO_SR(ui);
    24612466    }

  • kernel/longrat.h

    r5f3629 r896561  
    9393#endif
    9494
     95//#define SI_THREADS 1
     96#ifdef SI_THREADS
     97#define ALLOC_RNUMBER() (number)malloc(sizeof(snumber));
     98static inline number ALLOC0_RNUMBER() { number n=(number)malloc(sizeof(snumber));memset(n,0,sizeof(snumber));return n; }
     99#define FREE_RNUMBER(x) free(x)
     100#else
    95101extern omBin rnumber_bin;
     102#define ALLOC_RNUMBER() (number)omAllocBin(rnumber_bin)
     103#define ALLOC0_RNUMBER() (number)omAlloc0Bin(rnumber_bin)
     104#define FREE_RNUMBER(x) omFreeBin((ADDRESS)x, rnumber_bin)
     105#endif
    96106
    97107// in-place operations

  • kernel/longrat0.cc

    r5f3629 r896561  
    5555    return s;
    5656  }
    57   *a=(number)omAllocBin(rnumber_bin);
     57  *a=(number)ALLOC_RNUMBER();
    5858  {
    5959    (*a)->s = 3;
    … …  
    8686    {
    8787      mpz_clear(z);
    88       omFreeBin(*a,rnumber_bin);
     88      FREE_RNUMBER(*a);
    8989      *a=INT_TO_SR(0);
    9090    }
    … …  
    9696    }
    9797    else
    98       nlNormalize(*a);
     98    {
     99      number aa=*a;
     100      nlNormalize(aa);
     101      *a=aa;
     102    }
    99103  }
    100104  return s;

  • kernel/longtrans.cc

    r5f3629 r896561  
    5858#endif
    5959
    60 static number ntdGcd( number a, number b, const ring r) { return nacInit(1,r); }
     60//static number ntdGcd( number a, number b, const ring r) { return nacInit(1,r); }
    6161/*2
    6262*  sets the appropriate operators
    … …  
    943943    p = pInit();
    944944    pNext(p)=NULL;
    945     nNew(&pGetCoeff(p));
     945    //nNew(&pGetCoeff(p));
    946946    int i;
    947     for(i=pVariables;i;i--)
    948        pSetExp(p,i, 0);
    949     if (rRing_has_Comp(currRing)) pSetComp(p, 0);
     947    //for(i=pVariables;i;i--) pSetExp(p,i, 0); // done by pInit
     948    //if (rRing_has_Comp(currRing)) pSetComp(p, 0); // done by pInit
    950949    napoly pa=NULL;
    951950    lnumber pan;
    … …  
    953952    {
    954953      assume(oldRing->algring!=NULL);
    955       pGetCoeff(p)=(number)omAlloc0Bin(rnumber_bin);
     954      pGetCoeff(p)=(number)ALLOC0_LNUMBER();
    956955      pan=(lnumber)pGetCoeff(p);
    957956      pan->s=2;
    … …  
    10281027  if (x->n!=NULL)
    10291028  {
    1030     lnumber rr=(lnumber)omAlloc0Bin(rnumber_bin);
     1029    lnumber rr=ALLOC0_LNUMBER();
    10311030    rr->z=p_Copy(x->n,r->algring);
    10321031    rr->s = 2;
    … …  
    10391038{
    10401039  lnumber x=(lnumber)n;
    1041   lnumber rr=(lnumber)omAlloc0Bin(rnumber_bin);
     1040  lnumber rr=ALLOC0_LNUMBER();
    10421041  rr->z=p_Copy(x->z,r->algring);
    10431042  rr->s = 2;
    … …  
    10591058      poly z=p_Init(r->algring);
    10601059      pSetCoeff0(z,c);
    1061       lnumber l = (lnumber)omAllocBin(rnumber_bin);
     1060      lnumber l = (lnumber)ALLOC_LNUMBER();
    10621061      l->z = z;
    10631062      l->s = 2;
    … …  
    11031102number  ntPar(int i)
    11041103{
    1105   lnumber l = (lnumber)omAllocBin(rnumber_bin);
     1104  lnumber l = ALLOC_LNUMBER();
    11061105  l->s = 2;
    11071106  l->z = p_ISet(1,nacRing);
    … …  
    11601159    p_Delete(&(l->z),r->algring);
    11611160    p_Delete(&(l->n),r->algring);
    1162     omFreeBin((ADDRESS)l,  rnumber_bin);
     1161    FREE_LNUMBER(l);
    11631162  }
    11641163  *p = NULL;
    … …  
    11741173  lnumber erg;
    11751174  lnumber src = (lnumber)p;
    1176   erg = (lnumber)omAlloc0Bin(rnumber_bin);
     1175  erg = ALLOC_LNUMBER();
    11771176  erg->z = p_Copy(src->z, nacRing);
    11781177  erg->n = p_Copy(src->n, nacRing);
    … …  
    11851184  lnumber erg;
    11861185  lnumber src = (lnumber)p;
    1187   erg = (lnumber)omAlloc0Bin(rnumber_bin);
     1186  erg = ALLOC_LNUMBER();
    11881187  erg->z = p_Copy(src->z,r->algring);
    11891188  erg->n = p_Copy(src->n,r->algring);
    … …  
    12051204  lnumber b = (lnumber)lb;
    12061205  #ifdef LDEBUG
    1207   omCheckAddrSize(a,sizeof(snumber));
    1208   omCheckAddrSize(b,sizeof(snumber));
     1206  omCheckAddrSize(a,sizeof(slnumber));
     1207  omCheckAddrSize(b,sizeof(slnumber));
    12091208  #endif
    12101209  if (b->n!=NULL) x = pp_Mult_qq(a->z, b->n,nacRing);
    … …  
    12171216    return (number)NULL;
    12181217  }
    1219   lu = (lnumber)omAllocBin(rnumber_bin);
     1218  lu = ALLOC_LNUMBER();
    12201219  lu->z=res;
    12211220  if (a->n!=NULL)
    … …  
    12731272
    12741273  #ifdef LDEBUG
    1275   omCheckAddrSize(a,sizeof(snumber));
    1276   omCheckAddrSize(b,sizeof(snumber));
     1274  omCheckAddrSize(a,sizeof(slnumber));
     1275  omCheckAddrSize(b,sizeof(slnumber));
    12771276  #endif
    12781277
    … …  
    12871286    return (number)NULL;
    12881287  }
    1289   lu = (lnumber)omAllocBin(rnumber_bin);
     1288  lu = ALLOC_LNUMBER();
    12901289  lu->z=res;
    12911290  if (a->n!=NULL)
    … …  
    13281327
    13291328  #ifdef LDEBUG
    1330   omCheckAddrSize(a,sizeof(snumber));
    1331   omCheckAddrSize(b,sizeof(snumber));
     1329  omCheckAddrSize(a,sizeof(slnumber));
     1330  omCheckAddrSize(b,sizeof(slnumber));
    13321331  #endif
    13331332  ntTest(la);
    13341333  ntTest(lb);
    13351334
    1336   lo = (lnumber)omAllocBin(rnumber_bin);
     1335  lo = ALLOC_LNUMBER();
    13371336  lo->z = pp_Mult_qq(a->z, b->z,nacRing);
    13381337
    … …  
    13611360  if(lo->z==NULL)
    13621361  {
    1363     omFreeBin((ADDRESS)lo, rnumber_bin);
     1362    FREE_LNUMBER(lo);
    13641363    lo=NULL;
    13651364  }
    … …  
    13941393  }
    13951394#ifdef LDEBUG
    1396   omCheckAddrSize(a,sizeof(snumber));
    1397   omCheckAddrSize(b,sizeof(snumber));
     1395  omCheckAddrSize(a,sizeof(slnumber));
     1396  omCheckAddrSize(b,sizeof(slnumber));
    13981397#endif
    13991398  assume(a->z!=NULL && b->z!=NULL);
    14001399  assume(a->n==NULL && b->n==NULL);
    1401   res = (lnumber)omAllocBin(rnumber_bin);
     1400  res = ALLOC_LNUMBER();
    14021401  res->z = napCopy(a->z);
    14031402  res->n = napCopy(b->z);
    … …  
    14301429  }
    14311430  #ifdef LDEBUG
    1432   omCheckAddrSize(a,sizeof(snumber));
    1433   omCheckAddrSize(b,sizeof(snumber));
     1431  omCheckAddrSize(a,sizeof(slnumber));
     1432  omCheckAddrSize(b,sizeof(slnumber));
    14341433  #endif
    1435   lo = (lnumber)omAllocBin(rnumber_bin);
     1434  lo = ALLOC_LNUMBER();
    14361435  if (b->n!=NULL)
    14371436    lo->z = pp_Mult_qq(a->z, b->n,nacRing);
    … …  
    14881487  }
    14891488  #ifdef LDEBUG
    1490   omCheckAddrSize(b,sizeof(snumber));
     1489  omCheckAddrSize(b,sizeof(slnumber));
    14911490  #endif
    1492   lo = (lnumber)omAlloc0Bin(rnumber_bin);
     1491  lo = ALLOC0_LNUMBER();
    14931492  lo->s = b->s;
    14941493  if (b->n!=NULL)
    … …  
    16231622    return s;
    16241623  }
    1625   *p = (number)omAlloc0Bin(rnumber_bin);
     1624  *p = (number)ALLOC0_LNUMBER();
    16261625  a = (lnumber)*p;
    16271626  a->z = x;
    16281627  if(a->z==NULL)
    16291628  {
    1630     omFreeBin((ADDRESS)*p, rnumber_bin);
     1629    FREE_LNUMBER(a);
    16311630    *p=NULL;
    16321631  }
    … …  
    17041703{
    17051704  lnumber a = (lnumber)za;
    1706   number t;
    17071705  if (a==NULL) return FALSE;
    17081706#ifdef LDEBUG
    1709   omCheckAddrSize(a,sizeof(snumber));
     1707  omCheckAddrSize(a,sizeof(slnumber));
    17101708  if (a->z==NULL)
    17111709  {
    … …  
    17231721  }
    17241722#if 0
     1723  number t;
    17251724  x = a->z;
    17261725  y = a->n;
    … …  
    17631762  if (a==NULL) return FALSE;
    17641763#ifdef LDEBUG
    1765   omCheckAddrSize(a,sizeof(snumber));
     1764  omCheckAddrSize(a,sizeof(slnumber));
    17661765  if (a->z==NULL)
    17671766  {
    … …  
    18021801
    18031802  lnumber x, y;
    1804   lnumber result = (lnumber)omAlloc0Bin(rnumber_bin);
     1803  lnumber result = ALLOC0_LNUMBER();
    18051804
    18061805  x = (lnumber)a;
    … …  
    21632162  lnumber a = (lnumber)la;
    21642163  lnumber b = (lnumber)lb;
    2165   result = (lnumber)omAlloc0Bin(rnumber_bin);
     2164  result = ALLOC0_LNUMBER();
    21662165  ntTest(la);
    21672166  ntTest(lb);
    … …  
    22052204{
    22062205  if (npIsZero(c)) return NULL;
    2207   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     2206  lnumber l=ALLOC_LNUMBER();
    22082207  l->s=2;
    22092208  l->z=(napoly)p_Init(nacRing);
    … …  
    22212220{
    22222221  if (nlIsZero(c)) return NULL;
    2223   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     2222  lnumber l=ALLOC_LNUMBER();
    22242223  l->s=0;
    22252224  l->z=(napoly)p_Init(nacRing);
    … …  
    22352234{
    22362235  if (npIsZero(c)) return NULL;
    2237   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     2236  lnumber l=ALLOC_LNUMBER();
    22382237  l->s=2;
    22392238  l->z=(napoly)p_Init(nacRing);
    … …  
    22532252  number n=npInit(i,ntMapRing);
    22542253  if (npIsZero(n)) return NULL;
    2255   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     2254  lnumber l=ALLOC_LNUMBER();
    22562255  l->s=2;
    22572256  l->z=(napoly)p_Init(nacRing);
    … …  
    22702269  if (npIsZero(n)) return NULL;
    22712270  npTest(n);
    2272   lnumber l=(lnumber)omAllocBin(rnumber_bin);
     2271  lnumber l=ALLOC_LNUMBER();
    22732272  l->s=2;
    22742273  l->z=(napoly)p_Init(nacRing);
    … …  
    22842283{
    22852284  if (c==NULL) return NULL;
    2286   lnumber erg= (lnumber)omAlloc0Bin(rnumber_bin);
     2285  lnumber erg= ALLOC0_LNUMBER();
    22872286  lnumber src =(lnumber)c;
    22882287  erg->s=src->s;
    … …  
    23982397    {
    23992398      assume(oldRing->algring!=NULL);
    2400       pGetCoeff(p)=(number)omAlloc0Bin(rnumber_bin);
     2399      pGetCoeff(p)=(number)ALLOC0_LNUMBER();
    24012400      pan=(lnumber)pGetCoeff(p);
    24022401      pan->s=2;
    … …  
    24682467  if (x->n!=NULL)
    24692468  {
    2470     lnumber rr=(lnumber)omAlloc0Bin(rnumber_bin);
     2469    lnumber rr=ALLOC0_LNUMBER();
    24712470    rr->z=p_Copy(x->n,r->algring);
    24722471    rr->s = 2;
    … …  
    24792478{
    24802479  lnumber x=(lnumber)n;
    2481   lnumber rr=(lnumber)omAlloc0Bin(rnumber_bin);
     2480  lnumber rr=ALLOC0_LNUMBER();
    24822481  rr->z=p_Copy(x->z,r->algring);
    24832482  rr->s = 2;
    … …  
    24922491    return TRUE;
    24932492  #ifdef LDEBUG
    2494   omCheckAddrSize(a, sizeof(snumber));
     2493  omCheckAddrSize(a, sizeof(slnumber));
    24952494  #endif
    24962495  napoly p = x->z;

  • kernel/modulop.cc

    r5f3629 r896561  
    449449  else
    450450    e=0;
    451   res = (number)omAllocBin(rnumber_bin);
     451  res = ALLOC_RNUMBER();
    452452#if defined(LDEBUG)
    453453  res->debug=123456;
    … …  
    486486  iz=mpz_fdiv_ui(dest,npPrimeM);
    487487  mpz_clear(dest);
    488   omFreeBin((ADDRESS)res, rnumber_bin);
    489488  if(res->s==0)
    490489    iz=(long)npDiv((number)iz,(number)in);
     490  FREE_RNUMBER(res);
    491491  return (number)iz;
    492492}

  • kernel/p_Procs_Set.h

    r5f3629 r896561  
    152152
    153153  SetProcs(field, length, ord);
     154  extern poly p_Mult_nn_pthread(poly p, const number n, const ring r);
    154155  #ifdef NV_OPS
    155156  if ((field==FieldZp) && (r->ch>NV_MAX_PRIME))
    … …  
    160161    SetProcs_nv(FieldZp, length,ord); // p_Delete, p_ShallowCopyDelete...
    161162  }
     163  #endif
     164  #ifdef SI_THREADS
     165  else if (field==FieldQ) p_Procs->p_Mult_nn=p_Mult_nn_pthread;
    162166  #endif
    163167  CheckProc(p_Copy);

  • kernel/polys1.cc

    r5f3629 r896561  
    498498        poly h = pNext(ph);
    499499        p_Mult_nn(ph,k,currRing);
    500         if (!rField_has_simple_inverse()) pNormalize(ph);
     500        pNormalize(ph);
    501501      }
    502502      nDelete(&k);

  • kernel/ring.cc

    r5f3629 r896561  
    30553055    assume (r != NULL );
    30563056    assume (exp_limit > 1);
    3057     BOOLEAN omitted_degree = FALSE;
    30583057    int bits;
    30593058
    … …  
    48394838{
    48404839  int r_blocks = rBlocks(currRing);
    4841   int i;
    48424840
    48434841  assume(b1 == ringorder_c || b1 == ringorder_C ||
    … …  
    48844882}
    48854883
    4886 
    48874884ring rAssure_InducedSchreyerOrdering(const ring r, BOOLEAN complete = TRUE, int sgn = 1)
    48884885{ // TODO: ???? Add leading Syz-comp ordering here...????
    … …  
    53245321  return size;
    53255322}
    5326 static int rReallocM1(ring r, ring src, int size, int pos)
    5327 {
    5328   r->order=(int*)omReallocSize(r->order, size*sizeof(int), (size-1)*sizeof(int));
    5329   r->block0=(int*)omReallocSize(r->block0, size*sizeof(int), (size-1)*sizeof(int));
    5330   r->block1=(int*)omReallocSize(r->block1, size*sizeof(int), (size-1)*sizeof(int));
    5331   r->wvhdl=(int_ptr*)omReallocSize(r->wvhdl,size*sizeof(int_ptr), (size-1)*sizeof(int_ptr));
    5332   for(int k=pos+1; k<size; k++) r->wvhdl[k]=r->wvhdl[k+1];
    5333   size--;
    5334   return size;
    5335 }
     5323//static int rReallocM1(ring r, ring src, int size, int pos)
     5324//{
     5325//  r->order=(int*)omReallocSize(r->order, size*sizeof(int), (size-1)*sizeof(int));
     5326//  r->block0=(int*)omReallocSize(r->block0, size*sizeof(int), (size-1)*sizeof(int));
     5327//  r->block1=(int*)omReallocSize(r->block1, size*sizeof(int), (size-1)*sizeof(int));
     5328//  r->wvhdl=(int_ptr*)omReallocSize(r->wvhdl,size*sizeof(int_ptr), (size-1)*sizeof(int_ptr));
     5329//  for(int k=pos+1; k<size; k++) r->wvhdl[k]=r->wvhdl[k+1];
     5330//  size--;
     5331//  return size;
     5332//}
    53365333static void rOppWeight(int *w, int l)
    53375334{
    … …  
    55205517        r->block0[j]=rOppVar(r, src->block1[i]);
    55215518        r->block1[j]=rOppVar(r, src->block0[i]);
    5522         r->wvhdl[j]=r->wvhdl[j+1]; r->wvhdl[j+1]=r->wvhdl[j+1]=NULL;
     5519        r->wvhdl[j]=r->wvhdl[j+1]; r->wvhdl[j+1]=NULL;
    55235520        rOppWeight(r->wvhdl[j], r->block1[j]-r->block0[j]);
    55245521        j++;
    … …  
    57295726  matrix D0 = src->GetNC()->D;
    57305727
    5731 
    5732   poly p = NULL;
    5733   number n = NULL;
    5734 
    57355728  // map C and D into dest
    57365729  for (int i = 1; i < N; i++)
Note: See TracChangeset for help on using the changeset viewer.