Changeset 85a342 in git


Ignore:
Timestamp:
Mar 20, 2009, 9:34:43 PM (14 years ago)
Author:
Christian Eder
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
Children:
b99421493f70c5a40b8867d26b69e67529326c2b
Parents:
1122b22f64e9f1131b42cfe12973c0ad95b902b2
Message:
criterion2() for previously computed elements deleted


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

Legend:

Unmodified
Added
Removed
  • kernel/f5gb.cc

    r1122b2 r85a342  
    2525#include "f5data.h"
    2626#include "f5lists.h"
     27#include "timer.h"
    2728int reductionsToZero   =  0;
    2829
     
    232233        // testing both new labels by the F5 Criterion
    233234        //critPairs->print();
    234         if(!criterion2(gPrev->getFirst()->getIndex(), u1, newElement, rules, rTag)
    235            && !criterion2(gPrev->getFirst()->getIndex(), u2, temp, rules, rTag) &&
    236            !criterion1(gPrev,u1,newElement,lTag) && !criterion1(gPrev,u2,temp,lTag)) {
     235        if(!criterion2(gPrev->getFirst()->getIndex(), u2, temp, rules, rTag)
     236           && !criterion2(gPrev->getFirst()->getIndex(), u1, newElement, rules, rTag)
     237           && !criterion1(gPrev,u1,newElement,lTag) && !criterion1(gPrev,u2,temp,lTag)) {
    237238            // if they pass the test, add them to CList critPairs, having the LPoly with greater
    238239            // label as first element in the CPair
     
    300301        //pWrite(ppMult_qq(t,l->getTerm()));
    301302        //Print("%d\n\n",l->getIndex());
    302         /*
    303303        while(testNode->getIndex() < idx) { // && NULL != testNode->getLPoly()) {
    304304            //pWrite(testNode->getPoly());
     
    311311            testNode    =   testNode->getNext();
    312312        }
    313         */
     313        /*
    314314        ideal testId    =   idInit(idx-1,1);
    315315        for(i=0;i<idx-1;i++) {
     
    329329            return true;
    330330        }
     331        */
    331332        return false;
    332333    }
     
    886887*/
    887888ideal F5main(ideal id, ring r) {
     889    int timer   =   initTimer();
     890    startTimer();
    888891    int i,j,k,l;
    889892    int gbLength;
     
    986989                gPrev    =   new LList(pOne,1,gbPrev->m[0]);
    987990                gPrev->insert(pOne,1,gbPrev->m[1]);
    988                 poly tempPoly = pInit();
    989                 pLcm(pHead(gbPrev->m[0]),pHead(gbPrev->m[1]),tempPoly);
    990                 tempPoly    =   pDivide(tempPoly,pOne());
    991                 pSetCoeff(tempPoly,nOne);
     991                //poly tempPoly = pInit();
     992                //pLcm(pHead(gbPrev->m[0]),pHead(gbPrev->m[1]),tempPoly);
     993                //tempPoly    =   pDivide(tempPoly,pOne());
     994                //pSetCoeff(tempPoly,nOne);
    992995                rules    =   new RList();
    993996                rTag     =   new RTagList(rules->getFirst());
     
    995998                //Print("%p\n",rules->getFirst());
    996999                //pWrite(tempPoly);
    997                 rules->insert(2,tempPoly);
    998                 rTag->insert(rules->getFirst());
     1000                //rules->insert(2,tempPoly);
     1001                //rTag->insert(rules->getFirst());
    9991002                //Print("%p\n",rules->getFirst());
    10001003                //Print("%p\n",rTag->getFirst());
     
    10091012                        //pWrite(gbPrev->m[k]);
    10101013                        //pWrite(gbPrev->m[l]);
    1011                         poly tempPoly2  =   pOne();
    1012                         pLcm(pHead(gbPrev->m[k]),pHead(gbPrev->m[l]),tempPoly2);
    1013                         tempPoly2   =   pDivide(tempPoly2,pOne());
    1014                         pSetCoeff(tempPoly2,nOne);
     1014                        //poly tempPoly2  =   pOne();
     1015                        //pLcm(pHead(gbPrev->m[k]),pHead(gbPrev->m[l]),tempPoly2);
     1016                        //tempPoly2   =   pDivide(tempPoly2,pOne());
     1017                        //pSetCoeff(tempPoly2,nOne);
    10151018                        //pWrite(tempPoly2);
    1016                         rules->insert(k+1,tempPoly2);
     1019                        //rules->insert(k+1,tempPoly2);
    10171020                    }
    1018                     rTag->insert(rules->getFirst());
     1021                    //rTag->insert(rules->getFirst());
    10191022                }
    10201023            }
     
    10401043        //idShow(gbPrev);
    10411044    Print("\n\nNumber of zero-reductions:  %d\n",reductionsToZero);
     1045    timer   =   getTimer();
     1046    Print("Time for computations: %d/1000 seconds\n",timer);
    10421047    //LNode* temp    =   gPrev->getFirst();
    10431048    //while(NULL != temp) {
Note: See TracChangeset for help on using the changeset viewer.