Changeset a0350e9 in git for kernel/f5gb.cc
- Timestamp:
- Jan 15, 2009, 6:44:24 PM (15 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 2653d3c2809f0b96d6da699e54066fdb07623b92
- Parents:
- 5ec4a8bc9c1e8e23c5f5f95235a4e4839659b682
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/f5gb.cc
r5ec4a8 ra0350e9 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: f5gb.cc,v 1.1 5 2008-12-27 13:50:05ederc Exp $ */4 /* $Id: f5gb.cc,v 1.16 2009-01-15 17:44:23 ederc Exp $ */ 5 5 /* 6 6 * ABSTRACT: f5gb interface … … 24 24 #include "lpolynomial.h" 25 25 #include "lists.h" 26 27 28 /*29 ================================================30 computation of ONE polynomial as global variable31 ================================================32 */33 poly one_poly() {34 poly one = pInit();35 pSetCoeff(one,nInit(1));36 return one;37 }38 39 26 40 27 … … 80 67 ================================================== 81 68 */ 82 LList* F5inc(long* i, poly* f_i, LList* g_prev) { 83 poly one = pInit(); 84 pSetCoeff(one, nInit(1)); 85 static poly ONE = one; 86 //poly ONE = pOne(); 69 LList* F5inc(int* i, poly* f_i, LList* g_prev, poly* ONE) { 87 70 LList* g_curr = g_prev; 88 LNode* prev_last = g_prev->getLast(); //tags the end of g_prev->only 1 list for g_prev & g_curr 89 g_curr->append(&ONE,i,f_i); 71 g_curr->insert(ONE,i,f_i); 90 72 91 73 return g_curr; … … 98 80 */ 99 81 ideal F5main(ideal id, ring r) { 100 101 static poly ONE = pOne(); 102 long i,j; 103 82 int i,j; 83 const int idElems = IDELEMS(id); 84 // 1 polynomial for defining initial labels & further tests 85 static poly ONE = pOne(); 86 87 104 88 // definition of one-polynomial as global constant ONE 105 89 //poly one = pInit(); … … 122 106 idShow(id); 123 107 i = 1; 124 //lp->setIndex(&i);125 //lp->setTerm(&ONE);126 //lp->setPoly(&id->m[0]);127 LPoly* lp = new LPoly(&ONE,&i,&ONE);128 129 108 // only for debugging 130 long k = 2; 131 LList* g_prev = new LList(&ONE,&k,&id->m[2]); 132 LNode* current; 133 LPoly* lp2 = new LPoly(&ONE,&k,&ONE); 109 //LNode* current; 134 110 //LList* g_curr = new LList(lp); 135 k = 134; 136 g_prev->append(&ONE,&k,&id->m[3]); 137 g_prev->append(&ONE,&k,&id->m[3]); 138 g_prev->append(lp2); 139 g_prev->append(lp2); 140 g_prev->append(lp); 141 g_prev->append(&ONE,&k,&id->m[1]); 142 g_prev->append(&ONE,&k,&id->m[3]); 143 g_prev->append(lp2); 144 i = g_prev->getLength(); 145 Print("%ld\n\n",i); 146 current = g_prev->getFirst(); 147 while(NULL != current) { 148 Print("Index: %ld\n",*(current->getLPoly()->getIndex())); 149 Print("Pointer comparison: %p , %p\n\n",g_prev->getFirst(),current); 150 current = current->getNext(); 151 } 111 //} 152 112 //for(i=2; i<IDELEMS(id); i++) { 153 113 //g_curr = F5inc(&i,&id->m[i],g_prev);
Note: See TracChangeset
for help on using the changeset viewer.