Changeset 06c40b in git


Ignore:
Timestamp:
Nov 13, 2009, 4:46:02 PM (15 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '3720ae8bfcff4a4649ee98a15552089151d2d59b')
Children:
a86cda223ebf7fa50fae36321ec45f3377bf4fb9
Parents:
3a7ee744399f1b56c9bf08b3b16e524f0f2f0de4
Message:
*hannes: minor cleanup in fast_map

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

Legend:

Unmodified
Added
Removed
  • kernel/fast_maps.cc

    r3a7ee74 r06c40b  
    511511      if ((p->f1!=NULL)&&(p->f2!=NULL))
    512512      {
    513 #if 0
    514         printf("found prod:");
    515         p_wrp(p->src,src_r);printf("=");
    516         p_wrp(p->f1->src,src_r);printf(" * ");
    517         p_wrp(p->f2->src,src_r);printf("\n");
    518 #endif
    519513        poly f1=p->f1->dest;
    520514        poly f2=p->f2->dest;
     
    537531      else
    538532      {
    539         //printf("compute "); p_wrp(p->src,src_r);printf("\n");
    540533        assume((p->f1==NULL) && (p->f2==NULL));
    541         //if(p->f1!=NULL) { printf(" but f1="); p_wrp(p->f1->src,src_r);printf("\n"); }
    542         //if(p->f2!=NULL) { printf(" but f2="); p_wrp(p->f2->src,src_r);printf("\n"); }
    543534        // no factorization provided, use the classical method:
    544535        p->dest=maPoly_EvalMon(p->src,src_r,dest_id->m,dest_r);
    545         //p_wrp(p->dest, dest_r); printf(" is p->dest\n");
    546536      }
    547537    } /* p->dest==NULL */
     
    585575  int i;
    586576  int dg = 0;
    587   poly ggt = NULL;
     577  poly ggt = p_Init(r);
    588578  q1 = p_Init(r);
    589579  q2 = p_Init(r);
    590   ggt=p_Init(r);
    591580
    592581  for (i=1;i<=r->N;i++)
     
    649638    if (q != NULL)
    650639    {
     640      int tmp_deg;
    651641      assume((q1!=NULL)&&(q2!=NULL));
    652       if (pDeg(q,r) > ggt_deg)
     642      if ((tmp_deg=pDeg(q,r)) > ggt_deg)
    653643      {
    654644        choice=iter;
     
    659649          p_LmFree(fq_p, r);
    660650        }
    661         ggt_deg = pDeg(q, r);
     651        ggt_deg = tmp_deg ; /*pDeg(q, r);*/
    662652        ggT = q;
    663653        fp_p = q1;
     
    729719        }
    730720      }
     721      else assume(ggT==NULL);
    731722    }
    732723    iter=iter->next;
Note: See TracChangeset for help on using the changeset viewer.