Changeset e4b56d in git
- Timestamp:
- Aug 6, 2009, 4:36:51 PM (14 years ago)
- Branches:
- (u'spielwiese', '91fdef05f09f54b8d58d92a472e9c4a43aa4656f')
- Children:
- 89314848080dc671e1de502749b8f7a2b6fc2302
- Parents:
- 424cd646332e11e54c822873f6276ff3ddcb8751
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/longalg.cc
r424cd64 re4b56d 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: longalg.cc,v 1.4 0 2009-08-06 10:18:33Singular Exp $ */4 /* $Id: longalg.cc,v 1.41 2009-08-06 14:36:51 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT: algebraic numbers … … 209 209 210 210 /*3 211 * returns ph * z212 */213 static void napMultN(napoly p, number z)214 {215 number t;216 217 while (p!=NULL)218 {219 t = nacMult(napGetCoeff(p), z);220 nacNormalize(t);221 nacDelete(&napGetCoeff(p),nacRing);222 napGetCoeff(p) = t;223 napIter(p);224 }225 }226 227 /*3228 211 * division with rest; f = g*q + r, returns r and destroy f 229 212 */ … … 325 308 nacNormalize(t); 326 309 t = nacNeg(t); 327 napMultN(qa, t);310 qa=p_Mult_nn(qa,t,nacRing); 328 311 nacDelete(&t,nacRing); 329 312 napNormalize(qa); … … 345 328 nacNormalize(napGetCoeff(r)); 346 329 t = nacInvers(napGetCoeff(r)); 347 napMultN(q, t);330 q=p_Mult_nn(q,t,nacRing); 348 331 napNormalize(q); 349 332 nacDelete(&t,nacRing); … … 374 357 t = nacInvers(napGetCoeff(r)); 375 358 //nacNormalize(t); 376 napMultN(q, t);359 q=p_Mult_nn(q,t,nacRing); 377 360 napNormalize(q); 378 361 nacDelete(&t,nacRing); … … 763 746 napCleardenom(y); 764 747 if (!nacIsOne(napGetCoeff(g))) 765 napMultN(y, napGetCoeff(g));748 y=p_Mult_nn(y,napGetCoeff(g),nacRing); 766 749 napDelete1(&g); 767 750 return y; … … 1780 1763 /* collect all denoms from y and multiply x and y by it */ 1781 1764 number n=napLcm(y); 1782 napMultN(x,n);1783 napMultN(y,n);1765 x=p_Mult_nn(x,n,nacRing); 1766 y=p_Mult_nn(y,n,nacRing); 1784 1767 nacDelete(&n,nacRing); 1785 1768 while(x!=NULL) … … 1830 1813 number h1 = nacInvers(napGetCoeff(y)); 1831 1814 nacNormalize(h1); 1832 napMultN(x, h1);1815 x=p_Mult_nn(x,h1,nacRing); 1833 1816 nacDelete(&h1,nacRing); 1834 1817 napDelete1(&y); … … 1875 1858 { 1876 1859 number n=napLcm(y); 1877 napMultN(x,n);1878 napMultN(y,n);1860 x=p_Mult_nn(x,n,nacRing); 1861 y=p_Mult_nn(y,n,nacRing); 1879 1862 nacDelete(&n,nacRing); 1880 1863 while(x!=NULL) … … 2064 2047 number aa=nacInvers(a); 2065 2048 nacDelete(&a,nacRing); 2066 napMultN(x,aa);2049 x=p_Mult_nn(x,aa,nacRing); 2067 2050 nacDelete(&aa,nacRing); 2068 2051 }
Note: See TracChangeset
for help on using the changeset viewer.