Changeset 6094f4 in git
- Timestamp:
- Oct 9, 2009, 2:21:25 PM (14 years ago)
- Branches:
- (u'spielwiese', '828514cf6e480e4bafc26df99217bf2a1ed1ef45')
- Children:
- aa742138e4cde2d55a4f1744cc4b0e6446fb67d4
- Parents:
- d8db707919fa01c928aa53fa8778bf431ce018ef
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/longalg.cc
rd8db70 r6094f4 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: longalg.cc,v 1.5 6 2009-10-08 17:02:33Singular Exp $ */4 /* $Id: longalg.cc,v 1.57 2009-10-09 12:21:25 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT: algebraic numbers … … 164 164 #endif 165 165 166 #define napInit(i) (napoly)p_ISet(i,nacRing)167 166 #define napSetCoeff(p,n) {n_Delete(&pGetCoeff(p),nacRing);pGetCoeff(p)=n;} 168 167 #define napDelete1(p) p_LmDelete((poly *)p, nacRing) … … 326 325 { 327 326 q = p_Mult_q(q, qa,nacRing); 328 q = napAdd(q, napInit(1));327 q = napAdd(q, p_ISet(1,nacRing)); 329 328 nacNormalize(pGetCoeff(r)); 330 329 t = nacInvers(pGetCoeff(r)); … … 339 338 } 340 339 q = p_Mult_q(q, napCopy(qa),nacRing); 341 q = napAdd(q, napInit(1));340 q = napAdd(q, p_ISet(1,nacRing)); 342 341 qa = napNeg(qa); 343 342 loop … … 670 669 number x, y; 671 670 if (!naIsChar0) 672 return napInit(1);671 return p_ISet(1,nacRing); 673 672 x = nacCopy(pGetCoeff(a)); 674 673 if (nacIsOne(x)) … … 709 708 || ((pNext(b)==NULL)&&(nacIsZero(pGetCoeff(b))))) 710 709 { 711 return napInit(1);710 return p_ISet(1,nacRing); 712 711 } 713 712 return napCopy(b); … … 731 730 y = a; 732 731 } 733 if (!naIsChar0) g = napInit(1);732 if (!naIsChar0) g = p_ISet(1,nacRing); 734 733 else g = napGcd0(x, y); 735 734 if (pNext(y)==NULL) … … 767 766 g=a; 768 767 h=b; 769 if (!naIsChar0) x = napInit(1);768 if (!naIsChar0) x = p_ISet(1,nacRing); 770 769 else x = napGcd0(g,h); 771 770 for (i=(naNumbOfPar-1); i>=0; i--) … … 913 912 lnumber l = (lnumber)omAllocBin(rnumber_bin); 914 913 l->s = 2; 915 l->z = napInit(1);914 l->z = p_ISet(1,nacRing); 916 915 napSetExp(l->z,i,1); 917 916 p_Setm(l->z,nacRing); … … 1339 1338 lo->z = napCopy(b->n); 1340 1339 else 1341 lo->z = napInit(1);1340 lo->z = p_ISet(1,nacRing); 1342 1341 x = b->z; 1343 1342 if ((!napIsConstant(x)) || !nacIsOne(pGetCoeff(x))) … … 1598 1597 p_Delete(&a->z,nacRing); 1599 1598 p_Delete(&a->n,nacRing); 1600 a->z = napInit(1);1599 a->z = p_ISet(1,nacRing); 1601 1600 a->n = NULL; 1602 1601 a->s = 2; … … 2003 2002 //if (((naMinimalPoly==NULL) && (naI==NULL)) || !naIsChar0) 2004 2003 //{ 2005 // result->z = napInit(1);2004 // result->z = p_ISet(1,nacRing); 2006 2005 // return (number)result; 2007 2006 //}
Note: See TracChangeset
for help on using the changeset viewer.