Changeset 069fcf in git


Ignore:
Timestamp:
Jun 12, 1998, 8:11:03 PM (26 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b52fc4b2495505785981d640dcf7eb3e456778ef')
Children:
24c10ac334406bcf0c360f28d2f85e95a6a58946
Parents:
7a5e7c01d23dede5210a4d5e11e3938902d7fc5f
Message:
* hannes: minor cleanups


git-svn-id: file:///usr/local/Singular/svn/trunk@2136 2c84dea3-7e68-4137-9b89-c4e89433aadc
Location:
Singular
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • Singular/weight.cc

    r7a5e7c r069fcf  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: weight.cc,v 1.7 1998-04-27 09:11:56 pohl Exp $ */
     4/* $Id: weight.cc,v 1.8 1998-06-12 18:11:02 Singular Exp $ */
    55
    66/*
     
    3232extern "C" double wFunctionalBuch(int *degw, int *lpol, int npol,
    3333       double *rel, double wx);
    34 extern "C" double * wDouble(void *adr);       
     34extern "C" double * wDouble(void *adr);
    3535extern "C" void wAdd(int *A, int mons, int kn, int xx);
    3636extern "C" void wNorm(int *degw, int *lpol, int npol, double *rel);
     
    3939extern "C" void wSecondSearch(int *A, int *x, int *lpol,
    4040        int npol, int mons, double *rel, double *fk);
    41 extern "C" void wGcd(int *x, int n);       
    42 extern double nsqr;       
     41extern "C" void wGcd(int *x, int n);
     42extern double wNsqr;
    4343#else
    4444short * ecartWeights=NULL;
    45 extern int pVariables;
    46 
     45
     46#ifdef ALIGN_8
     47#define wDouble(A) ((double *)A)
     48#else
    4749static double * wDouble(void *adr)
    4850{
     
    5052  return (double *)((i+7)&(~7));
    5153}
    52 
    53 double nsqr;
     54#endif
     55
     56double wNsqr;
    5457double (*wFunctional)(int *degw, int *lpol, int npol,
    5558       double *rel, double wx);
     
    7275  {
    7376    ecl = ecu = e1 = *ex++;
    74     for (j = lpol[i] - 1; j; j--)
     77    for (j = lpol[i] - 1; j!=0; j--)
    7578    {
    7679      ec = *ex++;
     
    9699    gecart *= ((double)5.0 - ghom);
    97100  }
    98   return (gfmax * gecart) / pow(wx, nsqr);
     101  return (gfmax * gecart) / pow(wx, wNsqr);
    99102}
    100103
     
    115118  {
    116119    ecu = ecl = *ex++;
    117     for (j = lpol[i] - 1; j; j--)
     120    for (j = lpol[i] - 1; j!=0; j--)
    118121    {
    119122      ec = *ex++;
     
    130133  if (ghom > (double)0.5)
    131134    gfmax *= ((double)1.0 - (ghom * ghom)) / (double)0.75;
    132   return gfmax / pow(wx, nsqr);
     135  return gfmax / pow(wx, wNsqr);
    133136}
    134137
     
    142145  if (xx == 1)
    143146  {
    144     for (i = mons; i; i--)
     147    for (i = mons; i!=0; i--)
    145148      *ex++ -= *B++;
    146149  }
    147150  else
    148151  {
    149     for (i = mons; i; i--)
     152    for (i = mons; i!=0; i--)
    150153      *ex++ -= (*B++) * xx;
    151154  }
     
    161164  if (xx == 1)
    162165  {
    163     for (i = mons; i; i--)
     166    for (i = mons; i!=0; i--)
    164167      *ex++ += *B++;
    165168  }
    166169  else
    167170  {
    168     for (i = mons; i; i--)
     171    for (i = mons; i!=0; i--)
    169172      *ex++ += (*B++) * xx;
    170173  }
     
    251254
    252255  y = (double)x[n];
    253   for (i = n - 1; i; i--)
     256  for (i = n - 1; i!=0; i--)
    254257    y *= (double)x[i];
    255258  return y;
     
    267270  degw = A + (n * mons);
    268271  fo2 = fo1 = (double)1.0e10;
    269   for (i1 = n; i1; i1--)
     272  for (i1 = n; i1!=0; i1--)
    270273  {
    271274    if (x[i1] > 1)
     
    280283        k0 = i1;
    281284      }
    282       for (i2 = i1; i2; i2--)
     285      for (i2 = i1; i2!=0; i2--)
    283286      {
    284287        if (x[i2] > 1)
     
    416419  kopt = k = g = 0;
    417420  min = 1000000;
    418   for (i = n; i; i--)
     421  for (i = n; i!=0; i--)
    419422  {
    420423    c = xopt[i];
     
    438441  {
    439442    s2 = s1 = (double)0.0;
    440     for(i = n; i; i--)
     443    for(i = n; i!=0; i--)
    441444    {
    442445      c = xopt[i];
     
    464467    }
    465468  }
    466   for(i = n; i; i--)
     469  for(i = n; i!=0; i--)
    467470  {
    468471    x[i] = 1;
     
    494497  {
    495498    ecu = *ex++;
    496     for (j = lpol[i] - 1; j; j--)
     499    for (j = lpol[i] - 1; j!=0; j--)
    497500    {
    498501      ec = *ex++;
     
    537540
    538541
    539 static void wInit(polyset s, int sl, int mons, int *A) 
     542static void wInit(polyset s, int sl, int mons, int *A)
    540543{
    541544  int  n, a, i, j, *B, *C;
     
    574577          C += mons;
    575578        }
    576         q = pNext(q);
     579        pIter(q);
    577580      }
    578581    }
     
    592595  wDimensions(s, sl, lpol, &npol, &mons);
    593596  xopt = x + (n + 1);
    594   for (i = n; i; i--)
     597  for (i = n; i!=0; i--)
    595598    xopt[i] = 1;
    596599  if (mons==0)
     
    606609  degw = A + (n * mons);
    607610  memset(degw, 0, mons * sizeof(int));
    608   for (i = n; i; i--)
     611  for (i = n; i!=0; i--)
    609612    wAdd(A, mons, i, 1);
    610613  wNorm(degw, lpol, npol, rel);
     
    618621  memcpy(x + 1, xopt + 1, n * sizeof(int));
    619622  memset(degw, 0, mons * sizeof(int));
    620   for (i = n; i; i--)
     623  for (i = n; i!=0; i--)
    621624  {
    622625    x[i] *= 16;
     
    627630  if (fx >= eps)
    628631  {
    629     for (i = n; i; i--)
     632    for (i = n; i!=0; i--)
    630633      xopt[i] = 1;
    631634  }
     
    638641//      wSimple(x, n);
    639642//      memset(degw, 0, mons * sizeof(int));
    640 //      for (i = n; i; i--)
     643//      for (i = n; i!=0; i--)
    641644//        wAdd(A, mons, i, x[i]);
    642645//      eps = wPrWeight(x, n);
     
    662665  *eweight = 0;
    663666  n = pVariables;
    664   nsqr = (double)2.0 / (double)n;
     667  wNsqr = (double)2.0 / (double)n;
    665668  if (pOrdSgn == -1)
    666669    wFunctional = wFunctionalMora;
     
    669672  x = (int * )Alloc(2 * (n + 1) * sizeof(int));
    670673  wCall(s, sl, x);
    671   for (i = n; i; i--)
     674  for (i = n; i!=0; i--)
    672675    eweight[i] = x[i + n + 1];
    673676  Free((ADDRESS)x, 2 * (n + 1) * sizeof(int));
     
    687690  sl = IDELEMS(F) - 1;
    688691  n = pVariables;
    689   nsqr = (double)2.0 / (double)n;
     692  wNsqr = (double)2.0 / (double)n;
    690693  wFunctional = wFunctionalBuch;
    691694  x = (int * )Alloc(2 * (n + 1) * sizeof(int));
    692695  wCall(s, sl, x);
    693   for (i = n; i; i--)
     696  for (i = n; i!=0; i--)
    694697    (*iv)[i-1] = x[i + n + 1];
    695698  Free((ADDRESS)x, 2 * (n + 1) * sizeof(int));
     
    711714  int i;
    712715
    713   for (i=1;i<=len;i++)
     716  for (i=pVariables;i>len;i--)
     717    s[i]= 1;
     718  for (;i>0;i--)
    714719    s[i]= (*iv)[i-1];
    715   for (;i<=pVariables;i++)
    716     s[i]= 1;
    717720  return s;
    718721}
  • Singular/weight0.c

    r7a5e7c r069fcf  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: weight0.c,v 1.3 1998-01-27 15:40:48 pohl Exp $ */
     4/* $Id: weight0.c,v 1.4 1998-06-12 18:11:03 Singular Exp $ */
    55
    66/*
     
    1818double wFunctionalBuch(int *degw, int *lpol, int npol,
    1919       double *rel, double wx);
    20 double * wDouble(void *adr);       
    2120void wAdd(int *A, int mons, int kn, int xx);
    2221void wNorm(int *degw, int *lpol, int npol, double *rel);
     
    2524void wSecondSearch(int *A, int *x, int *lpol,
    2625        int npol, int mons, double *rel, double *fk);
    27 void wGcd(int *x, int n);       
     26void wGcd(int *x, int n);
    2827/*0 implementation*/
    2928
     
    3130extern int pVariables;
    3231
     32#ifdef ALIGN_8
     33#define wDouble(A) ((double *)A)
     34#else
    3335double * wDouble(void *adr)
    3436{
     
    3638  return (double *)((i+7)&(~7));
    3739}
    38 
    39 double nsqr;
     40#endif
     41
     42double wNsqr;
    4043double (*wFunctional)(int *degw, int *lpol, int npol,
    4144       double *rel, double wx);
     
    5861  {
    5962    ecl = ecu = e1 = *ex++;
    60     for (j = lpol[i] - 1; j; j--)
     63    for (j = lpol[i] - 1; j!=0; j--)
    6164    {
    6265      ec = *ex++;
     
    8285    gecart *= ((double)5.0 - ghom);
    8386  }
    84   return (gfmax * gecart) / pow(wx, nsqr);
     87  return (gfmax * gecart) / pow(wx, wNsqr);
    8588}
    8689
     
    101104  {
    102105    ecu = ecl = *ex++;
    103     for (j = lpol[i] - 1; j; j--)
     106    for (j = lpol[i] - 1; j!=0 ; j--)
    104107    {
    105108      ec = *ex++;
     
    116119  if (ghom > (double)0.5)
    117120    gfmax *= ((double)1.0 - (ghom * ghom)) / (double)0.75;
    118   return gfmax / pow(wx, nsqr);
     121  return gfmax / pow(wx, wNsqr);
    119122}
    120123
     
    126129  B = A + ((kn - 1) * mons);
    127130  ex = A + (pVariables * mons);
     131  i = mons;
    128132  if (xx == 1)
    129133  {
    130     for (i = mons; i; i--)
     134    for (/* i=mons */; i!=0 ; i--)
    131135      *ex++ -= *B++;
    132136  }
    133137  else
    134138  {
    135     for (i = mons; i; i--)
     139    for (/* i=mons */; i!=0 ; i--)
    136140      *ex++ -= (*B++) * xx;
    137141  }
     
    145149  B = A + ((kn - 1) * mons);
    146150  ex = A + (pVariables * mons);
     151  i = mons;
    147152  if (xx == 1)
    148153  {
    149     for (i = mons; i; i--)
     154    for (/* i=mons */; i!=0 ; i--)
    150155      *ex++ += *B++;
    151156  }
    152157  else
    153158  {
    154     for (i = mons; i; i--)
     159    for (/* i=mons */; i!=0 ; i--)
    155160      *ex++ += (*B++) * xx;
    156161  }
     
    237242
    238243  y = (double)x[n];
    239   for (i = n - 1; i; i--)
     244  for (i = n - 1; i!=0 ; i--)
    240245    y *= (double)x[i];
    241246  return y;
     
    253258  degw = A + (n * mons);
    254259  fo2 = fo1 = (double)1.0e10;
    255   for (i1 = n; i1; i1--)
     260  for (i1 = n; i1!=0 ; i1--)
    256261  {
    257262    if (x[i1] > 1)
     
    266271        k0 = i1;
    267272      }
    268       for (i2 = i1; i2; i2--)
     273      for (i2 = i1; i2!=0 ; i2--)
    269274      {
    270275        if (x[i2] > 1)
     
    388393      return;
    389394  }
    390   for (i = n; i; i--)
     395  for (i = n; i!=0 ; i--)
    391396    x[i] /= b;
    392397}
     
    402407  kopt = k = g = 0;
    403408  min = 1000000;
    404   for (i = n; i; i--)
     409  for (i = n; i!=0 ; i--)
    405410  {
    406411    c = xopt[i];
     
    424429  {
    425430    s2 = s1 = (double)0.0;
    426     for(i = n; i; i--)
     431    for(i = n; i!=0 ; i--)
    427432    {
    428433      c = xopt[i];
     
    450455    }
    451456  }
    452   for(i = n; i; i--)
     457  for(i = n; i!=0 ; i--)
    453458  {
    454459    x[i] = 1;
     
    480485  {
    481486    ecu = *ex++;
    482     for (j = lpol[i] - 1; j; j--)
     487    for (j = lpol[i] - 1; j!=0 ; j--)
    483488    {
    484489      ec = *ex++;
Note: See TracChangeset for help on using the changeset viewer.